Molecular characterization and gene expression of thyrotropin receptor (TSHR) and a truncated TSHR-like in Senegalese sole

Transcriptome characterization of BPG axis and expression profiles of ovarian steroidogenesis-related genes in the Japanese sardine

Background

The clupeoid fishes are ecologically and commercially important fish species worldwide that exhibit a high level of population fluctuation, accompanied by alteration of reproductive traits. However, knowledge about their reproductive physiology in order to understand mechanisms underlying such population dynamics is limited. The endocrine system along with the brain–pituitary–gonadal (BPG) axis is critical for regulating reproduction. The aims of this study were to provide transcript data and genes related to the BPG axis, and to characterize the expression profiles of ovarian steroidogenesis-related genes in the Japanese sardine (Sardinops melanostictus, Clupeidae).

Results

RNA sequencing was performed using the sardine brain, pituitary, and gonad in both sexes. A total of 290,119 contigs were obtained and 115,173 non-redundant ORFs were annotated. The genes differentially expressed between ovary and testis were strongly associated with GO terms related to gamete production. The tissue-specific profile of the abundance of transcripts was characterized for the major regulators in the BPG axis, such as gonadotropin-releasing hormone, gonadotropin, and steroidogenic enzyme. By comparing between ovary and testis, out of eight different 17β-hydroxysteroid dehydrogenase (Hsd17b) genes identified, higher hsd17b7 expression was found in testis, whereas higher expression of hsd17b8, hsd17b10, hsd17b12a, and hsd17b12b was found in ovary. The cDNAs encoding key endocrine factors in the ovarian steroidogenic pathway were cloned, sequenced, and quantitatively assayed. In the pituitary, follicle-stimulating hormone beta peaked during vitellogenesis, while luteinizing hormone beta peaked at the completion of vitellogenesis. In the ovary, follicle-stimulating hormone receptor and luteinizing hormone receptor were upregulated from mid- to late phase of vitellogenesis. Furthermore, three steroidogenic enzyme genes (cyp11a1, cyp17a1, and cyp19a1a) gradually increased their expression during ovarian development, accompanying a rise in serum estradiol-17β, while 3β-hydroxysteroid dehydrogenase and steroidogenic acute regulatory protein did not change significantly.

Conclusions

This is the first report of deep RNA sequencing analysis of Japanese sardine, in which many key genes involved in the BPG axis were identified. Expression profiles of ovarian steroidogenesis-related genes provide a molecular basis of the physiological processes underlying ovarian development in the sardine. Our study will be a valuable resource for clarifying the molecular biology of clupeoid fishes.

Teleost Metamorphosis: The Role of Thyroid Hormone

In most teleosts, metamorphosis encompasses a dramatic post-natal developmental process where the free-swimming larvae undergo a series of morphological, cellular and physiological changes that enable the larvae to become a fully formed, albeit sexually immature, juvenile fish. In all teleosts studied to date thyroid hormones (TH) drive metamorphosis, being the necessary and sufficient factors behind this developmental transition. During metamorphosis, negative regulation of thyrotropin by thyroxine (T4) is relaxed allowing higher whole-body levels of T4 that enable specific responses at the tissue/cellular level. Higher local thyroid cellular signaling leads to cell-specific responses that bring about localized developmental events. TH orchestrate in a spatial-temporal manner all local developmental changes so that in the end a fully functional organism arises. In bilateral teleost species, the most evident metamorphic morphological change underlies a transition to a more streamlined body. In the pleuronectiform lineage (flatfishes), these metamorphic morphological changes are more dramatic. The most evident is the migration of one eye to the opposite side of the head and the symmetric pelagic larva development into an asymmetric benthic juvenile. This transition encompasses a dramatic loss of the embryonic derived dorsal-ventral and left-right axis. The embryonic dorsal-ventral axis becomes the left-right axis, whereas the embryonic left-right axis becomes, irrespectively, the dorsal-ventral axis of the juvenile animal. This event is an unparalleled morphological change in vertebrate development and a remarkable display of the capacity of TH-signaling in shaping adaptation and evolution in teleosts. Notwithstanding all this knowledge, there are still fundamental questions in teleost metamorphosis left unanswered: how the central regulation of metamorphosis is achieved and the neuroendocrine network involved is unclear; the detailed cellular and molecular events that give rise to the developmental processes occurring during teleost metamorphosis are still mostly unknown. Also in flatfish, comparatively little is still known about the developmental processes behind asymmetric development. This review summarizes the current knowledge on teleost metamorphosis and explores the gaps that still need to be challenged.

Unravelling paralogous gene expression dynamics during three-spined stickleback embryogenesis

Development requires the implementation of a plethora of molecular mechanisms, involving a large set of genes to ensure proper cell differentiation, morphogenesis of tissues and organs as well as the growth of the organism. Genome duplication and resulting paralogs are considered to provide the raw genetic materials important for new adaptation opportunities and boosting evolutionary innovation. The present study investigated paralogous genes, involved in three-spined stickleback (Gasterosteus aculeatus) development. Therefore, the transcriptomes of five early stages comprising developmental leaps were explored. Obtained expression profiles reflected the embryo's needs at different stages. Early stages, such as the morula stage comprised transcripts mainly involved in energy requirements while later stages were mostly associated with GO terms relevant to organ development and morphogenesis. The generated transcriptome profiles were further explored for differential expression of known and new paralogous genes. Special attention was given to hox genes, with hoxa13a being of particular interest and to pigmentation genes where itgb1, involved in the melanophore development, displayed a complementary expression pattern throughout studied stages. Knowledge obtained by untangling specific paralogous gene functions during development might not only significantly contribute to the understanding of teleost ontogenesis but might also shed light on paralogous gene evolution.

Effects of the soya isoflavone genistein in early life stages of the Senegalese sole, Solea senegalensis: Thyroid, estrogenic and metabolic biomarkers

This study examines the effects induced by environmentally relevant concentrations of the isoflavone genistein (3mg/L and 10mg/L) during early life stages of the Senegalese sole. Throughout the hypothalamus-pituitary-thyroid (HPT) axis, several neurohormonal regulatory thyroid signalling patterns (thyroglobulin/Tg, thyroid peroxidase/TPO, transthyretin/TTR, thyroid receptors/TRβ, and iodothrynonine deiodinases, Dio2 and Dio3) were analysed. Furthermore, the expression patterns of estrogen receptor ERβ and haemoprotein Cyp1a were also evaluated. In the control larvae, progressive increases of constitutive hormonal signalling pathways have been evidenced from the pre-metamorphosis phase onwards, reaching the highest expression basal levels at the metamorphosis (Tg, TPO, Dio2) and/or during post-metamorphosis (TTR, TRβ, ERβ). When the early larvae were exposed to both genistein concentrations (3mg/L and 10mg/L), a statistically significant down-regulation of TPO, TTR and Tg mRNA levels was clearly detected at the metamorphic stages. In addition, the Dio2 and Dio3 transcript expression levels were also down and up-regulated when exposed to both genistein concentrations. In the larvae exposed to genistein, no statistically significant responses were recorded for the TRβ expression patterns. Nevertheless, the ERβ and Cyp1a transcript levels were up-regulated at the middle metamorphic stage (S2, at 16 dph) in the larvae exposed to high genistein concentrations and, only the ERβ was down-regulated (S1, at 12dph) at the lower doses. Finally, all these pointed out imbalances were only temporarily disrupted by exposure to genistein, since most of the modulated transcriptional signals (i.e. up or down-regulation) were quickly restored to the baseline levels. Additionally, the control and genistein-exposed Senegalese sole specimens showed characteristic ontogenetic patterns and completely suitable for an optimal development, metamorphosis, and growth.

Transcriptional regulation of genes involved in retinoic acid metabolism in Senegalese sole larvae

The aim of this study was the characterization of transcriptional regulatory pathways mediated by retinoic acid (RA) in Senegalese sole larvae. For this purpose, pre-metamorphic larvae were treated with a low concentration of DEAB, an inhibitor of RALDH enzyme, until the end of metamorphosis. No differences in growth, eye migration or survival were observed. Nevertheless, gene expression analysis revealed a total of 20 transcripts differentially expressed during larval development and only six related with DEAB treatments directly involved in RA metabolism and actions (rdh10a, aldh1a2, crbp1, igf2r, rarg and cyp26a1) to adapt to a low-RA environment. In a second experiment, post-metamorphic larvae were exposed to the all-trans RA (atRA) observing an opposite regulation for those genes involved in RA synthesis and degradation (rdh10a, aldh1a2, crbp1 and cyp26a1) as well as other related with thyroid- (dio2) and IGF-axes (igfbp1, igf2r and igfbp5) to balance RA levels. In a third experiment, DEAB-pretreated post-metamorphic larvae were exposed to atRA and TTNPB (a specific RAR agonist). Both drugs down-regulated rdh10a and aldh1a2 and up-regulated cyp26a1 expression demonstrating their important role in RA homeostasis. Moreover, five retinoic receptors that mediate RA actions, the thyroid receptor thrb, and five IGF binding proteins changed differentially their expression. Overall, this study demonstrates that exogenous RA modulates the expression of some genes involved in the RA synthesis, degradation and cellular transport through RAR-mediated regulatory pathways establishing a negative feedback regulatory mechanism necessary to balance endogenous RA levels and gradients.

Flatfish metamorphosis: a hypothalamic independent process?

Anuran and flatfish metamorphosis are tightly regulated by thyroid hormones that are the necessary and sufficient factors that drive this developmental event. In the present study whole mount in situ hybridization (WISH) and quantitative PCR in sole are used to explore the central regulation of flatfish metamorphosis. Central regulation of the thyroid in vertebrates is mediated by the hypothalamus-pituitary-thyroid (HPT) axis. Teleosts diverge from other vertebrates as hypothalamic regulation in the HPT axis is proposed to be through hypothalamic inhibition although the regulatory factor remains enigmatic. The dynamics of the HPT axis during sole metamorphosis revealed integration between the activity of the thyrotrophes in the pituitary and the thyroid follicles. No evidence was found supporting a role for thyroid releasing hormone (trh) or corticotrophin releasing hormone (crh) in hypothalamic control of TH production during sole metamorphosis. Intriguingly the results of the present study suggest that neither hypothalamic trh nor crh expression changes during sole metamorphosis and raises questions about the role of these factors and the hypothalamus in regulation of thyrotrophs.

Multiple thyrotropin β-subunit and thyrotropin receptor-related genes arose during vertebrate evolution

Thyroid-stimulating hormone (TSH) is composed of a specific β subunit and an α subunit that is shared with the two pituitary gonadotropins. The three β subunits derive from a common ancestral gene through two genome duplications (1R and 2R) that took place before the radiation of vertebrates. Analysis of genomic data from phylogenetically relevant species allowed us to identify an additional Tshβ subunit-related gene that was generated through 2R. This gene, named Tshβ2, present in cartilaginous fish, little skate and elephant shark, and in early lobe-finned fish, coelacanth and lungfish, was lost in ray-finned fish and tetrapods. The absence of a second type of TSH receptor (Tshr) gene in these species suggests that both TSHs act through the same receptor. A novel Tshβ sister gene, named Tshβ3, was generated through the third genomic duplication (3R) that occurred early in the teleost lineage. Tshβ3 is present in most teleost groups but was lostin tedraodontiforms. The 3R also generated a second Tshr, named Tshrb. Interestingly, the new Tshrb was translocated from its original chromosomic position after the emergence of eels and was then maintained in its new position. Tshrb was lost in tetraodontiforms and in ostariophysians including zebrafish although the latter species have two TSHs, suggesting that TSHRb may be dispensable. The tissue distribution of duplicated Tshβs and Tshrs was studied in the European eel. The endocrine thyrotropic function in the eel would be essentially mediated by the classical Tshβ and Tshra, which are mainly expressed in the pituitary and thyroid, respectively. Tshβ3 and Tshrb showed a similar distribution pattern in the brain, pituitary, ovary and adipose tissue, suggesting a possible paracrine/autocrine mode of action in these non-thyroidal tissues. Further studies will be needed to determine the binding specificity of the two receptors and how these two TSH systems are interrelated.

Subfunctionalization of POMC paralogues in Senegalese sole (Solea senegalensis)

The precursor protein proopiomelanocortin (POMC) gives rise to a variety of biologically active peptides through cell-specific posttranslational processing. Two transcripts of pomc were found in the flatfish Solea senegalensis (ssePOMC-A and ssePOMC-B), that most likely represent subfunctionalized paralogues: ssePOMC-A lacks the N-terminal cleavage site for β-MSH, whereas ssePOMC-B cannot yield ACTH and completely lacks the opioid consensus sequence in the β-END region. An analysis of nucleotide substitution rates shows that the POMC-derived peptides possess well-conserved regions under purifying selection, except the β-END derived from POMC-B, which has undergone positive selection. The calculated K(s) values for ssePOMC-A versus ssePOMC-B and zebrafish POMCαversus zebrafish POMCβ are 0.40 and 0.72, respectively, indicating that the zebrafish POMC paralogues started to evolve almost twice as early in evolution, and that the Solea POMC paralogues arose independently from the whole genome duplication event that gave rise to the zebrafish paralogues. This makes ssePOMC-B the first identified POMCα orthologue that lacks the opioid consensus. Furthermore, pomc-a expression is down-regulated in chronic stressed S. senegalensis juveniles, whereas pomc-b expression levels remain unaffected, indicating different physiological roles for both POMC paralogues. The distribution of functional POMC-derived peptide hormones over two pomc genes in S. senegalensis suggests subfunctionalization of the paralogues, a relevant notion when studying POMC function in endocrine responses.

TSH receptor function is required for normal thyroid differentiation in zebrafish

TSH is the primary physiological regulator of thyroid gland function. The effects of TSH on thyroid cells are mediated via activation of its membrane receptor [TSH receptor (TSHR)]. In this study, we examined functional thyroid differentiation in zebrafish and characterized the role of TSHR signaling during thyroid organogenesis. Cloning of a cDNA encoding zebrafish Tshr showed conservation of primary structure and functional properties between zebrafish and mammalian TSHR. In situ hybridization confirmed that the thyroid is the major site of tshr expression during zebrafish development. In addition, we identified tpo, iyd, duox, and duoxa as novel thyroid differentiation markers in zebrafish. Temporal analyses of differentiation marker expression demonstrated the induction of an early thyroid differentiation program along with thyroid budding, followed by a delayed onset of duox and duoxa expression coincident with thyroid hormone synthesis. Furthermore, comparative analyses in mouse and zebrafish revealed for the first time a thyroid-enriched expression of cell death regulators of the B-cell lymphoma 2 family during early thyroid morphogenesis. Knockdown of tshr function by morpholino microinjection into embryos did not affect early thyroid morphogenesis but caused defects in later functional differentiation. The thyroid phenotype observed in tshr morphants at later stages comprised a reduction in number and size of functional follicles, down-regulation of differentiation markers, as well as reduced thyroid transcription factor expression. A comparison of our results with phenotypes observed in mouse models of defective TSHR and cAMP signaling highlights the value of zebrafish as a model to enhance the understanding of functional differentiation in the vertebrate thyroid.