The pituitary gland of the European eel reveals massive expression of genes involved in the melanocortin system

Pituitary multi-hormone cells in mammals and fish: history, origin, and roles

The vertebrate pituitary is a dynamic organ, capable of adapting its hormone secretion to different physiological demands. In this context, endocrinologists have debated for the past 40 years if endocrine cells are mono- or multi-hormonal. Since its establishment, the dominant "one cell, one hormone" model has been continuously challenged. In mammals, the use of advanced multi-staining approaches, sensitive gene expression techniques, and the analysis of tumor tissues have helped to quickly demonstrate the existence of pituitary multi-hormone cells. In fishes however, only recent advances in imaging and transcriptomics have enabled the identification of such cells. In this review, we first describe the history of the discovery of cells producing multiple hormones in mammals and fishes. We discuss the technical limitations that have led to uncertainties and debates. Then, we present the current knowledge and hypotheses regarding their origin and biological role, which provides a comprehensive review of pituitary plasticity.

Intracerebroventricular administration of α-melanocyte-stimulating hormone (α-MSH) enhances thigmotaxis and induces anxiety-like behavior in the goldfish Carassius auratus

α-Melanocyte-stimulating hormone (α-MSH) is a body pigmentation-regulating hormone secreted from the intermediate lobe of the pituitary in vertebrates. It is also produced in the brain, and acts as an anorexigenic neuropeptide involved in feeding regulation. In rodents, intracerebroventricular (ICV) administration of α-MSH has been shown to affect not only feeding behavior, but also psychomotor activity. However, there is still no information regarding the psychophysiological effects of α-MSH on behavior in fish. Therefore, we examined the effect of synthetic α-MSH on psychomotor activity in goldfish. Since this species prefers the edge to the central area of a tank, we used this as a preference test for assessing psychomotor activity. When α-MSH was administered ICV at 1 and 10 pmol g^-1 body weight (BW), the time spent in the edge area of a tank was prolonged at 10 pmol g^-1 BW. However, α-MSH at these doses did not affect locomotor activity. The action of α-MSH mimicked those of FG-7142 (a central-type benzodiazepine receptor (CBR) inverse agonist with an anxiogenic effect) at 10 pmol g^-1 BW and melanotan II (a melanocortin 4 receptor (MC4R) agonist) at 50 pmol g^-1 BW, whereas ICV administration of tofisopam (a CBR agonist with an anxiolytic effect) at 10 pmol g^-1 BW prolonged the time spent in the central area. The anxiogenic-like effect of α-MSH was abolished by treatment with the MC4R antagonist HS024 at 50 pmol g^-1 BW. These data indicate that α-MSH affects psychomotor activity in goldfish, and exerts an anxiogenic-like effect via the MC4R-signaling pathway.

Characterization of hormone-producing cell types in the teleost pituitary gland using single-cell RNA-seq

The pituitary is the vertebrate endocrine gland responsible for the production and secretion of several essential peptide hormones. These, in turn, control many aspects of an animal’s physiology and development, including growth, reproduction, homeostasis, metabolism, and stress responses. In teleost fish, each hormone is presumably produced by a specific cell type. However, key details on the regulation of, and communication between these cell types remain to be resolved. We have therefore used single-cell sequencing to generate gene expression profiles for 2592 and 3804 individual cells from the pituitaries of female and male adult medaka (Oryzias latipes), respectively. Based on expression profile clustering, we define 15 and 16 distinct cell types in the female and male pituitary, respectively, of which ten are involved in the production of a single peptide hormone. Collectively, our data provide a high-quality reference for studies on pituitary biology and the regulation of hormone production, both in fish and in vertebrates in general.

Measurement(s)	RNA-seq gene expression profiling assay
Technology Type(s)	tag based single cell RNA sequencing
Factor Type(s)	sex
Sample Characteristic - Organism	Oryzias latipes
Sample Characteristic - Environment	fresh water aquarium

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.16592621

Unravelling the changes during induced vitellogenesis in female European eel through RNA-Seq: What happens to the liver?

The life cycle of European eel (Anguilla anguilla), a catadromous species, is complex and enigmatic. In nature, during the silvering process prior to their long spawning migration, reproductive development is arrested, and they cease feeding. In studies of reproduction using hormonal induction, eels are equivalently not feed. Therefore, in female eels that undergo vitellogenesis, the liver plays different, essential roles being involved both in vitellogenins synthesis and in reallocating resources for the maintenance of vital functions, performing the transoceanic reproductive migration and completing reproductive development. The present work aimed at unravelling the major transcriptomic changes that occur in the liver during induced vitellogenesis in female eels. mRNA-Seq data from 16 animals (eight before induced vitellogenesis and eight after nine weeks of hormonal treatment) were generated and differential expression analysis was performed comparing the two groups. This analysis detected 1,328 upregulated and 1,490 downregulated transcripts. Overrepresentation analysis of the upregulated genes included biological processes related to biosynthesis, response to estrogens, mitochondrial activity and localization, while downregulated genes were enriched in processes related to morphogenesis and development of several organs and tissues, including liver and immune system. Among key genes, the upregulated ones included vitellogenin genes (VTG1 and VTG2) that are central in vitellogenesis, together with ESR1 and two novel genes not previously investigated in European eel (LMAN1 and NUPR1), which have been linked with reproduction in other species. Moreover, several upregulated genes, such as CYC1, ELOVL5, KARS and ACSS1, are involved in the management of the effect of fasting and NOTCH, VEGFA and NCOR are linked with development, autophagy and liver maintenance in other species. These results increase the understanding of the molecular changes that occur in the liver during vitellogenesis in this complex and distinctive fish species, bringing new insights on European eel reproduction and broodstock management.

De novo European eel transcriptome provides insights into the evolutionary history of duplicated genes in teleost lineages

Paralogues pairs are more frequently observed in eels (Anguilla sp.) than in other teleosts. The paralogues often show low phylogenetic distances; however, they have been assigned to the third round of whole genome duplication (WGD), shared by all teleosts (3R), due to their conserved synteny. The apparent contradiction of low phylogenetic difference and 3R conserved synteny led us to study the duplicated gene complement of the freshwater eels. With this aim, we assembled de novo transcriptomes of two highly relevant freshwater eel species: The European (Anguilla anguilla) and the Japanese eel (Anguilla japonica). The duplicated gene complement was analysed in these transcriptomes, and in the genomes and transcriptomes of other Actinopterygii species. The study included an assessment of neutral genetic divergence (4dTv), synteny, and the phylogenetic origins and relationships of the duplicated gene complements. The analyses indicated a high accumulation of duplications (1217 paralogue pairs) among freshwater eel genes, which may have originated in a WGD event after the Elopomorpha lineage diverged from the remaining teleosts, and thus not at the 3R. However, very similar results were observed in the basal Osteoglossomorpha and Clupeocephala branches, indicating that the specific genomic regions of these paralogues may still have been under tetrasomic inheritance at the split of the teleost lineages. Therefore, two potential hypotheses may explain the results: i) The freshwater eel lineage experienced an additional WGD to 3R, and ii) Some duplicated genomic regions experienced lineage specific rediploidization after 3R in the ancestor to freshwater eels. The supporting/opposing evidence for both hypotheses is discussed.

Polygenic and sex specific architecture for two maturation traits in farmed Atlantic salmon

BACKGROUND

A key developmental transformation in the life of all vertebrates is the transition to sexual maturity, whereby individuals are capable of reproducing for the first time. In the farming of Atlantic salmon, early maturation prior to harvest size has serious negative production impacts.

RESULTS

We report genome wide association studies (GWAS) using fish measured for sexual maturation in freshwater or the marine environment. Genotypic data from a custom 50 K single nucleotide polymorphism (SNP) array was used to identify 13 significantly associated SNP for freshwater maturation with the most strongly associated on chromosomes 10 and 11. A higher number of associations (48) were detected for marine maturation, and the two peak loci were found to be the same for both traits. The number and broad distribution of GWAS hits confirmed a highly polygenetic nature, and GWAS performed separately within males and females revealed sex specific genetic behaviour for loci co-located with positional candidate genes phosphatidylinositol-binding clathrin assembly protein-like (picalm) and membrane-associated guanylate kinase, WW and PDZ domain-containing protein 2 (magi2).

CONCLUSIONS

The results extend earlier work and have implications for future applied breeding strategies to delay maturation in this important aquaculture species.

Molecular biology and functional genomics of immune-endocrine interactions in the Japanese eel, Anguilla japonica

Immune-endocrine interactions are an important pathogen resistance mechanism in fish. We review the immune-endocrine interactions in the Japanese eel, Anguilla japonica, with special reference to high throughput gene sequencing. These data may be relevant to the significant decrease in the eel harvest in recent years and will aid in the selection of appropriate disease-resistant strains for aquaculture. More than 1000 sequences that whose expression in elvers responded to air exposure were identified through comprehensive gene expression analysis using next-generation sequencing. These included transcription factors within the MAPK pathway. Significant changes in expression after air exposure were detected by quantitative polymerase chain reaction analysis in many genes related to disease resistance. These factors include innate immune system factors and cytokines that interact with the endocrine system during the stress response. Other applications of immune-endocrine interactions in eel culture are discussed.

Rapid de novo assembly of the European eel genome from nanopore sequencing reads

We have sequenced the genome of the endangered European eel using the MinION by Oxford Nanopore, and assembled these data using a novel algorithm specifically designed for large eukaryotic genomes. For this 860 Mbp genome, the entire computational process takes two days on a single CPU. The resulting genome assembly significantly improves on a previous draft based on short reads only, both in terms of contiguity (N50 1.2 Mbp) and structural quality. This combination of affordable nanopore sequencing and light weight assembly promises to make high-quality genomic resources accessible for many non-model plants and animals.

A full-body transcriptome and proteome resource for the European common carp

Background

The common carp (Cyprinus carpio) is the oldest, most domesticated and one of the most cultured fish species for food consumption. Besides its economic importance, the common carp is also highly suitable for comparative physiological and disease studies in combination with the animal model zebrafish (Danio rerio). They are genetically closely related but offer complementary benefits for fundamental research, with the large body mass of common carp presenting possibilities for obtaining sufficient cell material for advanced transcriptome and proteome studies.

Results

Here we have used 19 different tissues from an F1 hybrid strain of the common carp to perform transcriptome analyses using RNA-Seq. For a subset of the tissues we also have performed deep proteomic studies. As a reference, we updated the European common carp genome assembly using low coverage Pacific Biosciences sequencing to permit high-quality gene annotation. These annotated gene lists were linked to zebrafish homologs, enabling direct comparisons with published datasets. Using clustering, we have identified sets of genes that are potential selective markers for various types of tissues. In addition, we provide a script for a schematic anatomical viewer for visualizing organ-specific expression data.

Conclusions

The identified transcriptome and proteome data for carp tissues represent a useful resource for further translational studies of tissue-specific markers for this economically important fish species that can lead to new markers for organ development. The similarity to zebrafish expression patterns confirms the value of common carp as a resource for studying tissue-specific expression in cyprinid fish. The availability of the annotated gene set of common carp will enable further research with both applied and fundamental purposes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3038-y) contains supplementary material, which is available to authorized users.

Using normalization to resolve RNA-Seq biases caused by amplification from minimal input

RNA-Seq has become a widely used method to study transcriptomes, and it is now possible to perform RNA-Seq on almost any sample. Nevertheless, samples obtained from small cell populations are particularly challenging, as biases associated with low amounts of input RNA can have strong and detrimental effects on downstream analyses. Here we compare different methods to normalize RNA-Seq data obtained from minimal input material. Using RNA from isolated medaka pituitary cells, we have amplified material from six samples before sequencing. Both synthetic and real data are used to evaluate different normalization methods to obtain a robust and reliable pipeline for analysis of RNA-Seq data from samples with very limited input material. The analysis outlined here shows that quantile normalization outperforms other more commonly used normalization procedures when using amplified RNA as input and will benefit researchers employing low amounts of RNA in similar experiments.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.