Zebrafish pituitary gene expression before and after sexual maturation

Enhancement of zebrafish sperm production via a large body-sized surrogate with germ cell transplantation

Zebrafish (Danio rerio) is a commonly-used vertebrate model species for many research areas. However, its low milt volume limits effective cryopreservation of sperm from a single individual and often precludes dividing a single semen sample to conduct multiple downstream procedures such as genomic DNA/RNA extraction and in-vitro fertilization. Here, we apply germ stem cell transplantation to increase zebrafish sperm production in a closely related larger species from the same subfamily, giant danio Devario aequipinnatus. The endogenous germ cell of the host is depleted by dead-end morpholino antisense oligonucleotide. Histology of the sterile gonad and quantitative PCR of gonadal tissue reveals all sterile giant danio develop the male phenotype. Spermatogonial cells of Tg(ddx4:egfp) transgenic zebrafish are transplanted into sterile giant danio larvae, and 22% of recipients (germline chimera) produce donor-derived sperm at sexual maturation. The germline chimera produce approximately three-fold the volume of sperm and 10-fold the spermatozoon concentration of the donor. The donor-derived sperm is functional and gives rise to viable progeny upon fertilization of donor oocytes. We show that the issue of low milt volume can be effectively addressed by employing a larger surrogate parent.

Transcriptome analyses reveal differentially expressed genes associated with development of the palatal organ in bighead carp (Hypophthalmichthys nobilis)

The palatal organ is a filter-feeding related organ and occupies a considerable proportion of the head of bighead carp (Hypophthalmichthys nobilis), a large cyprinid fish intensive aquaculture in Asia. In this study, we performed RNA-seq of the palatal organ during growth periods of two (M2), six (M6) and 15 (M15) months of age after hatching. The numbers of differentially expressed genes (DEGs) were 1384, 481 and 1837 for M2 VS M6, M6 VS M15 and M2 VS M15 respectively. The following signaling pathways of energy metabolism and cytoskeleton function were enriched, including ECM-receptor interaction, Cardiac muscle contraction, Steroid biosynthesis and PPAR signaling pathway. Several members of collagen family (col1a1, col2a1, col6a2, col6a3, col9a2), Laminin gamma 1 (lamc1), integrin alpha 1 (itga1), Fatty acid binding protein 2 (fads2) and lipoprotein lipase (lpl), and Protein tyrosine kinase 7 (Ptk7) are candidate genes for growth and development of basic tissues of the palatal organ. Furthermore, taste-related genes such as fgfrl1, fgf8a, fsta and notch1a were also identified, which may be involved in the development of taste buds of the palatal organ. The transcriptome data obtained in this study provide insights into the understanding functions and development mechanisms of palatal organ, and potential candidate genes that may be related to the genetic modulation of head size of bighead carp.

Age-associated changes in gene expression in the anterior pituitary glands of female Japanese black cattle

Proper functioning of the anterior pituitary (AP) gland is imperative, however, is suppressed by aging via unclear mechanisms. Therefore, we identified differentially expressed genes (DEGs) in the AP glands of Japanese Black young heifers (approximately 22 months old) compared to old cows (approximately 120 months old) via deep sequencing of the transcriptome (RNA-seq) to characterize potentially important pathways. The young and old AP glands expressed 20,171 annotated genes. Of the total transcripts per million, approximately 41.6% and 35.5% were the sum of seven AP hormone genes in young and old AP glands, respectively, with difference observed in the sum between the young and old AP glands (P < 0.05). Moreover, we identified 48 downregulated genes and 218 upregulated genes in old compared to young AP glands (P < 0.01, fold change > 120%). The DEGs included 1 cytokine (AIMP1), 3 growth factors (NRG2, PTN, and TGFB1), 1 receptor-associated protein gene (AGTRAP), and 10 receptor genes, including PRLHR and two orphan G-protein-coupled receptors (GPR156 and GPR176). Metascape analysis of the DEGs revealed "Peptide metabolic process," "Regulation of hormone levels," and "Peptide hormone processing" as enriched pathways. Furthermore, Ingenuity Pathway analysis of the DEGs revealed (1) a network of 24 genes (including GPR156 and PRLHR) named "Neurological disease, organismal injury and abnormalities, and psychological disorders", and (2) two canonical pathways (P < 0.01), namely "Huntington's disease signaling", and "AMPK signaling". Thus, the findings of the current study revealed relevant DEGs, while identifying important pathways that occur during aging in AP glands of female cattle.

Pituitary Gonadotropin Gene Expression During Induced Onset of Postsmolt Maturation in Male Atlantic Salmon: In Vivo and Tissue Culture Studies

Precocious male maturation causes reduced welfare and increased production costs in Atlantic salmon (Salmo salar) aquaculture. The pituitary produces and releases follicle-stimulating hormone (Fsh), the gonadotropin triggering puberty in male salmonids. However, little is known about how Fsh production is regulated in Atlantic salmon. We examined, in vivo and ex vivo, transcriptional changes of gonadotropin-related genes accompanying the initial steps of testis maturation, in pituitaries of males exposed to photoperiod and temperature conditions promoting maturation (constant light and 16°C). Pituitary fshb, lhb and gnrhr2bba transcripts increased in vivo in maturing males (gonado-somatic index > 0.1%). RNA sequencing (RNAseq) analysis using pituitaries from genetically similar males carrying the same genetic predisposition to mature, but differing by responding or not responding to stimulatory environmental conditions, revealed 144 differentially expressed genes, ~2/3rds being up-regulated in responders, including fshb and other pituitary hormones, steroid-related and other puberty-associated transcripts. Functional enrichment analyses confirmed gene involvement in hormone/steroid production and gonad development. In ex vivo studies, whole pituitaries were exposed to a selection of hormones and growth factors. Gonadotropin-releasing hormone (Gnrh), 17β-estradiol (E₂) and 11-ketotestosterone (11-KT) up-regulated gnrhr2bba and lhb, while fshb was up-regulated by Gnrh but down-regulated by 11-KT in pituitaries from immature males. Also pituitaries from maturing males responded to Gnrh and sex steroids by increased gnrhr2bba and lhb transcript levels, but fshb expression remained unchanged. Growth factors (inhibin A, activin A and insulin-like growth factor 1) did not change gnrhr2bba, lhb or fshb transcript levels in pituitaries either from immature or maturing males. Additional pituitary ex vivo studies on candidates identified by RNAseq showed that these transcripts were preferentially regulated by Gnrh and sex steroids, but not by growth factors, and that Gnrh/sex steroids were less effective when incubating pituitaries from maturing males. Our results suggest that a yet to be characterized mechanism up-regulating fshb expression in the salmon pituitary is activated in response to stimulatory environmental conditions prior to morphological signs of testis maturation, and that the transcriptional program associated with this mechanism becomes unresponsive or less responsive to most stimulators ex vivo once males had entered pubertal developmental in vivo.

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

3D Atlas of the Pituitary Gland of the Model Fish Medaka (Oryzias latipes)

In vertebrates, the anterior pituitary plays a crucial role in regulating several essential physiological processes via the secretion of at least seven peptide hormones by different endocrine cell types. Comparative and comprehensive knowledge of the spatial distribution of those endocrine cell types is required to better understand their physiological functions. Using medaka as a model and several combinations of multi-color fluorescence in situ hybridization, we present the first 3D atlas revealing the gland-wide distribution of seven endocrine cell populations: lactotropes, thyrotropes, Lh and Fsh gonadotropes, somatotropes, and pomca-expressing cells (corticotropes and melanotropes) in the anterior pituitary of a teleost fish. By combining in situ hybridization and immunofluorescence techniques, we deciphered the location of corticotropes and melanotropes within the pomca-expressing cell population. The 3D localization approach reveals sexual dimorphism of tshba-, pomca-, and lhb-expressing cells in the adult medaka pituitary. Finally, we show the existence of bi-hormonal cells co-expressing lhb-fshb, fshb-tshba and lhb-sl using single-cell transcriptomics analysis and in situ hybridization. This study offers a solid basis for future comparative studies of the teleost pituitary and its functional plasticity.

Screening and evaluating of long non-coding RNAs in prenatal and postnatal pituitary gland of sheep

Long non-coding RNAs are transcribed into RNA molecules that are >200 nucleotides in length. However, the expression and function analysis of lncRNAs in the sheep pituitary gland are still lacking. In this study, we identified 1755 lncRNAs (545 annotated lncRNAs and 1210 novel lncRNAs) from RNA-seq data in the pituitary gland of embryonic and adult sheep. A total of 235 lncRNAs were differentially expressed between embryonic and adult group. We verified the presence of some lncRNAs using RT-PCR and DNA sequencing, and identified some differentially expressed lncRNAs using qPCR. We also investigated the role of cis-acting lncRNAs on target genes. GO and KEGG enrichment analysis revealed that the target genes of lncRNAs were involved in the regulation of hormones secretion and some signaling pathways in the sheep pituitary gland. Our study provides comprehensive expression profiles of lncRNAs and valuable resource for understanding their function in the pituitary gland.

Evolution of proopiomelanocortin

Proopiomelanocortin (POMC) belongs to the opioid/orphanin gene family whose peptide precursors include either opioid (YGGF) or the orphanin/nociceptin core sequences (FGGF). In addition to POMC the family includes the proenkephalin (PENK), prodynorphin (PDYN), and nociceptin/proorphanin (PNOC) precursors. The opioid core sequence in POMC is incorporated by the β-endorphin that occupies the C-terminal region but this propeptide also exhibits at least two "alien" melanocortin core sequences (HFRW). An ACTH/MSH fragment merged into the opioid fragment not earlier than the two tetraploidizations of the vertebrate genome. Therefore, POMC exhibit a complex "evolutionary life" since the gene has coevolved together with two different receptor systems, i.e., opioid and melanocortin following a horse trading system. In this article, we summarize the different evolutionary hypotheses proposed for POMC evolution.

Combined transcriptome and proteome profiling reveals specific molecular brain signatures for sex, maturation and circalunar clock phase

Many marine animals, ranging from corals to fishes, synchronise reproduction to lunar cycles. In the annelid Platynereis dumerilii, this timing is orchestrated by an endogenous monthly (circalunar) clock entrained by moonlight. Whereas daily (circadian) clocks cause extensive transcriptomic and proteomic changes, the quality and quantity of regulations by circalunar clocks have remained largely elusive. By establishing a combined transcriptomic and proteomic profiling approach, we provide first systematic insight into the molecular changes in Platynereis heads between circalunar phases, and across sexual differentiation and maturation. Whereas maturation elicits large transcriptomic and proteomic changes, the circalunar clock exhibits only minor transcriptomic, but strong proteomic regulation. Our study provides a versatile extraction technique and comprehensive resources. It corroborates that circadian and circalunar clock effects are likely distinct and identifies key molecular brain signatures for reproduction, sex and circalunar clock phase. Examples include prepro-whitnin/proctolin and ependymin-related proteins as circalunar clock targets.

Forever young: Endocrinology of paedomorphosis in the Mexican axolotl (Ambystoma mexicanum)

The Mexican axolotl (Ambystoma mexicanum) is a salamander species that does not undergo metamorphosis, resulting in the retention of juvenile characteristics in the mature breeding stage (paedomorphosis). Here we review the endocrinological studies investigating the proximate cause of axolotl paedomorphosis with a focus on the hypothalamo-pituitary-thyroid (HPT) axis. It is well established that axolotl paedomorphosis is a consequence of low activity of the HPT axis. The pituitary hormone thyrotropin (TSH) is capable of inducing metamorphosis in the axolotl, which indicates that all processes and interactions in the HPT axis below the pituitary level are functional, but that TSH release is impaired. In metamorphosing species, TSH secretion is largely controlled by the hypothalamic neuropeptide corticotropin-releasing hormone (CRH), which seems to have lost its thyrotropic activity in the axolotl. However, preliminary experiments have not yet confirmed a role for faulty CRH signalling in axolotl paedomorphosis. Other hypothalamic factors and potential pituitary inhibitors need to be investigated to identify their roles in amphibian metamorphosis and axolotl paedomorphosis.

Deep sequencing of the transcriptome in the anterior pituitary of heifers before and after ovulation

We aimed to determine gene expression patterns in the anterior pituitary (AP) of heifers before and after ovulation via deep sequencing of the transcriptome (RNA-seq) to identify new genes and clarify important pathways. Heifers were slaughtered on the estrus day (pre-ovulation; n=5) or 3 days after ovulation (post-ovulation; n=5) for AP collection. We randomly selected 4 pre-ovulation and 4 post-ovulation APs, and the ribosomal RNA-depleted poly (A)+RNA were prepared to assemble next-generation sequencing libraries. The bovine APs expressed 12,769 annotated genes at pre- or post-ovulation. The sum of the reads per kilobase of exon model per million mapped reads (RPKM) values of all transcriptomes were 599,676 ± 38,913 and 668,209 ± 23,690, and 32.2 ± 2.6% and 44.0 ± 4.4% of these corresponded to the AP hormones in the APs of pre- and post-ovulation heifers, respectively. The bovine AP showed differential expression of 396 genes (P<0.05) in the pre- and post-ovulation APs. The 396 genes included two G-protein-coupled receptor (GPCR) genes (GPR61 and GPR153) and those encoding 13 binding proteins. The AP also expressed 259 receptor and other 364 binding proteins. Moreover, ingenuity pathway analysis for the 396 genes revealed (P=2.4 × 10⁻³) a canonical pathway linking GPCR to cytoskeleton reorganization, actin polymerization, microtubule growth, and gene expression. Thus, the present study clarified the novel genes found to be differentially expressed before and after ovulation and clarified an important pathway in the AP.

Novel hypophysiotropic AgRP2 neurons and pineal cells revealed by BAC transgenesis in zebrafish

The neuropeptide agouti-related protein (AgRP) is expressed in the arcuate nucleus of the mammalian hypothalamus and plays a key role in regulating food consumption and energy homeostasis. Fish express two agrp genes in the brain: agrp1, considered functionally homologous with the mammalian AgRP, and agrp2. The role of agrp2 and its relationship to agrp1 are not fully understood. Utilizing BAC transgenesis, we generated transgenic zebrafish in which agrp1- and agrp2-expressing cells can be visualized and manipulated. By characterizing these transgenic lines, we showed that agrp1-expressing neurons are located in the ventral periventricular hypothalamus (the equivalent of the mammalian arcuate nucleus), projecting throughout the hypothalamus and towards the preoptic area. The agrp2 gene was expressed in the pineal gland in a previously uncharacterized subgroup of cells. Additionally, agrp2 was expressed in a small group of neurons in the preoptic area that project directly towards the pituitary and form an interface with the pituitary vasculature, suggesting that preoptic AgRP2 neurons are hypophysiotropic. We showed that direct synaptic connection can exist between AgRP1 and AgRP2 neurons in the hypothalamus, suggesting communication and coordination between AgRP1 and AgRP2 neurons and, therefore, probably also between the processes they regulate.

Anatomical and functional gonadotrope networks in the teleost pituitary

Mammalian pituitaries exhibit a high degree of intercellular coordination; this enables them to mount large-scale coordinated responses to various physiological stimuli. This type of communication has not been adequately demonstrated in teleost pituitaries, which exhibit direct hypothalamic innervation and expression of luteinizing hormone (LH) and follicle-stimulating hormone (FSH) in distinct cell types. We found that in two fish species, namely tilapia and zebrafish, LH cells exhibit close cell-cell contacts and form a continuous network throughout the gland. FSH cells were more loosely distributed but maintained some degree of cell-cell contact by virtue of cytoplasmic processes. These anatomical differences also manifest themselves at the functional level as evidenced by the effect of gap-junction uncouplers on gonadotropin release. These substances abolished the LH response to gonadotropin-releasing hormone stimulation but did not affect the FSH response to the same stimuli. Dye transfer between neighboring LH cells provides further evidence for functional coupling. The two gonadotropins were also found to be differently packaged within their corresponding cell types. Our findings highlight the evolutionary origin of pituitary cell networks and demonstrate how the different levels of cell-cell coordination within the LH and FSH cell populations are reflected in their distinct secretion patterns.

Depletion of suppressor of cytokine signaling-1a causes hepatic steatosis and insulin resistance in zebrafish

Suppressor of cytokine signaling-1a (SOCS1a) is a member of the suppressor of cytokine signaling family, a group of related molecules that mediate the negative regulation of the JAK-STAT pathway. Here, we depleted SOCS1a using the transcription activator-like (TAL) effector nuclease (TALEN) technique to understand its physiological roles in zebrafish. Although elevated levels of JAK-STAT5 activation and erythropoiesis have been observed in socs1a-deficient zebrafish, these animals exhibited normal growth during the early stages. Socs1a-deficient zebrafish began to grow slowly with certain mortalities after 20 days postfertilization (dpf), whereas the heterozygous socs1a-deficient zebrafish exhibited enhanced somatic growth. Decreased adiposity, hepatic steatosis, and insulin resistance were observed in our socs1a-deficient adult zebrafish, which is similar to the lipodystrophy phenotypes observed in mammals. Comparative transcriptomic analyses revealed elevated levels of gluconeogenesis, lipolysis, and hypoxia-inducible response and decreased activities of lipogenesis and glycolysis in the hepatocytes of socs1a-deflicient adult zebrafish. Evident mitochondrial dysfunction has also been observed in hepatocytes from socs1a-deficient zebrafish. Taken together, our results suggest that the negative regulatory roles of SOCS1a on JAK-STAT5 signaling may be involved in the suppression of the erythropoiesis and growth hormone activities, which was also reflected in the enhanced somatic growth performance observed in the heterozygous socs1a-deficient fish. The differences in the effects caused by SOCS1a depletion on insulin sensitivity, lipid metabolism, and inflammatory responses between zebrafish and mammalian models observed here may reflect differences between the functional mechanisms of SOCS members in terrestrial mammals and aquatic teleosts.