Characterization and Expression Dynamics of Key Genes Involved in the Gilthead Sea Bream (Sparus aurata) Cortisol Stress Response during Early Ontogeny

Temperature-Biased miRNA Expression Patterns during European Sea Bass (Dicentrarchus labrax) Development

Environmental effects and, particularly, temperature changes have been demonstrated to influence the activity, function, and well-being of teleosts. Temperature may change seasonally in the wild, and in captivity under aquaculture operations. Moreover, climate change is expected to shift temperature profiles worldwide. MicroRNAs (miRNA) are important temperature-sensitive gene-expression regulators acting at the post-transcriptional level. They are known to be key regulators in development, reproduction, and immune responses. Therefore, early larval development of the European sea bass (Dicentrarchus labrax), one of the most extensively cultured species in Mediterranean aquaculture, was investigated at early rearing temperatures, i.e., 15, 17.5, and 20 °C, in regard to the impact of temperatures on miRNAs through sncRNA high-throughput sequencing but also at the phenotypic level in terms of growth, sex, vision, and skeletal deformities. Expression profiling revealed stage- and temperature-specific miRNA expression targeting genes with roles in reproduction and immune response mainly at the flexion and all-fins stages. Similar stage- and temperature-specific results were also observed concerning the number of rod cells and lower jaw elongation. The present work presents for the first time highly promising results on the influence of early rearing temperature at the post-transcriptional level during European sea bass development, with a putative impact on reproduction and immune response, as well as regarding teleost vision and larval development.

Differential Metabolic and Transcriptional Responses of Gilthead Seabream (Sparus aurata) Administered with Cortisol or Cortisol-BSA

Simple Summary

Cortisol is a key stress hormone in teleosts. Cortisol exerts its effects through genomic—and membrane-initiated mechanisms, however, the role of the latter in long-term stress responses is unknown. Here, we treated Sparus aurata with cortisol or cortisol-BSA (exclusive inductor to membrane-initiated effects) to emulate a long-term stress situation. We found that cortisol, but not cortisol-BSA, promotes energy substrate mobilization in the liver, together with the regulation of metabolism-related genes. We suggest that genomic cortisol actions exclusively participate in metabolic responses during prolonged treatment using cortisol in S. aurata. This study contributes to the current knowledge on cortisol’s involvement in stress responses in fish.

Abstract

Cortisol is the main glucocorticoid hormone promoting compensatory metabolic responses of stress in teleosts. This hormone acts through genomic and membrane-initiated actions to exert its functions inside the cell. Experimental approaches, using exogenous cortisol administration, confirm the role of this hormone during short (minutes to hours)- and long-term (days to weeks) responses to stress. The role of membrane-initiated cortisol signaling during long-term responses has been recently explored. In this study, Sparus aurata were intraperitoneally injected with coconut oil alone or coconut oil containing cortisol, cortisol-BSA, or BSA. After 3 days of treatment, plasma, liver, and skeletal muscle were extracted. Plasma cortisol, as well as metabolic indicators in the plasma and tissues collected, and metabolism-related gene expression, were measured. Our results showed that artificially increased plasma cortisol levels in S. aurata enhanced plasma glucose and triacylglycerols values as well as hepatic substrate energy mobilization. Additionally, cortisol stimulated hepatic carbohydrates metabolism, as seen by the increased expression of metabolism-related genes. All of these responses, observed in cortisol-administered fish, were not detected by replicating the same protocol and instead using cortisol-BSA, which exclusively induces membrane-initiated effects. Therefore, we suggest that after three days of cortisol administration, only genomic actions are involved in the metabolic responses in S. aurata.

Insights into early ontogenesis: characterization of stress and development key genes of pikeperch (Sander lucioperca) in vivo and in vitro

There are still numerous difficulties in the successful farming of pikeperch in the anthropogenic environment of various aquaculture systems, especially during early developmental steps in the hatchery. To investigate the physiological processes involved on the molecular level, we determined the basal expression patterns of 21 genes involved in stress and immune responses and early ontogenesis of pikeperch between 0 and 175 days post hatch (dph). Their transcription patterns most likely reflect the challenges of growth and feed conversion. The gene coding for apolipoprotein A (APOE) was strongly expressed at 0 dph, indicating its importance for yolk sac utilization. Genes encoding bone morphogenetic proteins 4 and 7 (BMP4, BMP7), creatine kinase M (CKM), and SRY-box transcription factor 9 (SOX9) were highly abundant during the peak phases of morphological changes and acclimatization processes at 4-18 dph. The high expression of genes coding for peroxisome proliferator-activated receptors alpha and delta (PPARA, PPARD) at 121 and 175 dph, respectively, suggests their importance during this strong growth phase of juvenile stages. As an alternative experimental model to replace further in vivo investigations of ontogenetically important processes, we initiated the first approach towards a long-lasting primary cell culture from whole pikeperch embryos. The present study provides a set of possible biomarkers to support the monitoring of pikeperch farming and provides a first basis for the establishment of a suitable cell model of this emerging aquaculture species.

Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism for cells to generate transcript variability and proteome diversity. No systematic investigation of AS events among different tissues in response to stressors is available for tilapia currently. In this study, AS among different tissues was identified and the cold stress-related AS events were explored in a Nile tilapia (Oreochromis niloticus) line based on 42 RNA-seq datasets using a bioinformatics pipeline. 14,796 (82.76%; SD = 2,840) of the expression genes showed AS events. The two most abundant AS types were alternative transcription start site (TSS) and terminal site (TTS) in tilapia. Testis, brain and kidney possess the most abundant AS gene number, while the blood, muscle and liver possess the least number in each tissue. Furthermore, 208 differentially alternative splicing (DAS) genes in heart and 483 DAS in brain in response to cold stress. The number of AS types for alternative exon end, exon skipping and retention of single intron increased significantly under cold stress. GO enrichment and pathway overrepresentation analysis indicated that many DAS genes, e.g., genes in circadian clock pathway, may influence expression of downstream genes under cold stress. Our study revealed that AS exists extensively in tilapia and plays an important role in cold adaption.

Early-life social environment alters juvenile behavior and neuroendocrine function in a highly social cichlid fish

Early-life experiences can shape adult behavior, with consequences for fitness and health, yet fundamental questions remain unanswered about how early-life social experiences are translated into variation in brain and behavior. The African cichlid fish Astatotilapia burtoni, a model system in social neuroscience, is well known for its highly plastic social phenotypes in adulthood. Here, we rear juveniles in either social groups or pairs to investigate the effects of early-life social environments on behavior and neuroendocrine gene expression. We find that both juvenile behavior and neuroendocrine function are sensitive to early-life effects. Behavior robustly co-varies across multiple contexts (open field, social cue investigation, and dominance behavior assays) to form a behavioral syndrome, with pair-reared juveniles towards the end of syndrome that is less active and socially interactive. Pair-reared juveniles also submit more readily as subordinates. In a separate cohort, we measured whole brain expression of stress and sex hormone genes. Expression of glucocorticoid receptor 1a was elevated in group-reared juveniles, supporting a highly-conserved role for the stress axis mediating early-life effects. The effect of rearing environment on androgen receptor α and estrogen receptor α expression was mediated by treatment duration (1 vs. 5 weeks). Finally, expression of corticotropin-releasing factor and glucocorticoid receptor 2 decreased significantly over time. Rearing environment also caused striking differences in gene co-expression, such that expression was tightly integrated in pair-reared juveniles but not group-reared or isolates. Together, this research demonstrates the important developmental origins of behavioral phenotypes and identifies potential behavioral and neuroendocrine mechanisms.

Transport and Recovery of Gilthead Seabream (Sparus aurata L.) Sedated With Clove Oil and MS-222: Effects on Stress Axis Regulation and Intermediary Metabolism

Transport processes between aquaculture facilities activate the stress response in fish. To deal with these situations, the hypothalamic-pituitary-interrenal (HPI) axis releases cortisol, leading to an increase in circulating energy resources to restore homeostasis. However, if the allostatic load generated exceeds fish tolerance limits, stress-related responses will compromise health and welfare of the animals. In this context, anesthetics have arisen as potential agents aiming to reduce negative effects of stress response. Here we assessed the effects of a sedative dose of clove oil (CO) and MS-222 on hallmarks involved in HPI axis regulation and energy management after simulated transport, and further recovery, in gilthead seabream (Sparus aurata L.) juveniles. Fish were placed in a mobile setup of water tanks where transport conditions were simulated for 6 h. Sedation doses of either CO (2.5 mg L⁻¹) or MS-222 (5 mg L⁻¹) were added in the water tanks. A control group without anesthetics was also included in the setup. Half of the animals (n = 12 per group) were sampled immediately after transport, while remaining animals were allowed to recover for 18 h in clean water tanks and then sampled. Our results showed that the HPI axis response was modified at peripheral level, with differences depending on the anesthetic employed. Head kidney gene-expressions related to cortisol production (star and cyp11b1) matched concomitantly with increased plasma cortisol levels immediately after transport in CO-sedated fish, but these levels remained constant in MS-222-sedated fish. Differential changes in the energy management of carbohydrates, lipids and amino acids, depending on the anesthetic employed, were also observed. The use of CO stimulated amino acids catabolism, while MS-222-sedated fish tended to consume liver glycogen and mobilize triglycerides. Further studies, including alternative doses of both anestethics, as well as the assessment of time-course HPI activation and longer recovery periods, are necessary to better understand if the use of clove oil and MS-222 is beneficial for S. aurata under these circumstances.

Identifying a Long QTL Cluster Across chrLG18 Associated with Salt Tolerance in Tilapia Using GWAS and QTL-seq

Understanding the genetic mechanism of osmoregulation is important for the improvement of salt tolerance in tilapia. In our previous study, we have identified a major quantitative trait locus (QTL) region located at 23.0 Mb of chrLG18 in a Nile tilapia line by QTL-seq. However, the conservation of these QTLs in other tilapia populations or species is not clear. In this study, we successfully investigated the QTLs associated with salt tolerance in a mass cross population from the GIFT line of Nile tilapia (Oreochromis niloticus) using a ddRAD-seq-based genome-wide association study (GWAS) and in a full-sib family from the Malaysia red tilapia strain (Oreochromis spp) using QTL-seq. Our study confirmed the major QTL interval that is located at nearly 23.0 Mb of chrLG18 in Nile tilapia and revealed a long QTL cluster across chrLG18 controlling for the salt-tolerant trait in both red tilapia and Nile tilapia. This is the first GWAS analysis on salt tolerance in tilapia. Our finding provides important insights into the genetic architecture of salinity tolerance in tilapia and supplies a basis for fine mapping QTLs, marker-assisted selection, and further detailed functional analysis of the underlying genes for salt tolerance in tilapia.

Expressions of Shell Matrix Protein Genes in the Pearl Sac and Its Correlation with Pearl Weight in the First 6 Months of Pearl Formation in Hyriopsis cumingii

Matrix proteins regulate crystal nucleation, morphology, and polymorphism during pearl biomineralization and have significant correlations with pearl quality traits in nucleated pearls. However, there is little information about the connection between pearl quality traits and matrix proteins in non-nucleated pearls. In this study, we analyzed CaCO₃ deposition during the first month of non-nucleated pearl formation and examined the expression patterns of ten shell matrix protein genes (Hcperlucin, hic31, silkmapin, hic22, hic74, hic52, HcTyr, HcCA3, hic24 and Hc-upsalin) in the pearl sac of Hyriopsis cumingii. During pearl formation, CaCO₃ crystals were initially deposited in a disorderly manner during days 12 and 15 of pearl formation. On days 18 and 21, CaCO₃ crystals gradually nucleated on an organic membrane, and the pattern of crystal deposition changed markedly. Between days 24 and 30, crystals similar to nacre tablets were deposited; they then grew and formed connections in a more orderly fashion, eventually forming the nacreous layer. We observed high expression levels of shell matrix proteins during the phases of disordered or ordered CaCO₃ deposition, suggesting they were involved in non-nucleated pearl formation. Furthermore, the expressions of nine matrix proteins were significantly correlated with pearl weight during the first 6 months after grafting. The prismatic-layer matrix protein hic31 and nacreous-layer matrix protein hic22 showed negative correlations with pearl weight, but the other seven nacreous-layer matrix proteins had significantly positive correlations with pearl weight. These results show the involvement of matrix proteins in pearl formation and in determination of quality traits.

Contribution of Non-canonical Cortisol Actions in the Early Modulation of Glucose Metabolism of Gilthead Sea Bream (Sparus aurata)

Teleost fish are exposed to diverse stressors in farming and wildlife conditions during their lifespan. Cortisol is the main glucocorticoid hormone involved in the regulation of their metabolic acclimation under physiological stressful conditions. In this context, increased plasma cortisol is associated with energy substrate mobilization from metabolic tissues, such as liver and skeletal muscle, to rapidly obtain energy and cope with stress. The metabolic actions of cortisol have primarily been attributed to its genomic/classic action mechanism involving the interaction with intracellular receptors, and regulation of stress-responsive genes. However, cortisol can also interact with membrane components to activate rapid signaling pathways. In this work, using the teleost fish gilthead sea bream (Sparus aurata) as a model, we evaluated the effects of membrane-initiated cortisol actions on the early modulation of glucose metabolism. For this purpose, S. aurata juveniles were intraperitoneally administrated with cortisol and with its membrane impermeable analog, cortisol-BSA. After 1 and 6 h of each treatment, plasma cortisol levels were measured, together with glucose, glycogen and lactate in plasma, liver and skeletal muscle. Transcript levels of corticosteroids receptors (gr1, gr2, and mr) and key gluconeogenesis (g6pc and pepck)- and glycolysis (pgam1 and aldo) related genes in the liver were also measured. Cortisol and cortisol-BSA administration increased plasma cortisol levels in S. aurata 1 h after administration. Plasma glucose levels enhanced 6 h after each treatment. Hepatic glycogen content decreased in the liver at 1 h of both cortisol and cortisol-BSA administration, while increased at 6 h due to cortisol but not in response to cortisol-BSA. Expression of gr1, g6pc, pgam1, and aldo were preferentially increased by cortisol-BSA in the liver. Taking all these results in consideration, we suggest that non-canonical cortisol mechanisms contribute to the regulation of the early glucose metabolism responses to stress in S. aurata.

Allostatic Load and Stress Physiology in European Seabass (Dicentrarchus labrax L.) and Gilthead Seabream (Sparus aurata L.)

The present study aimed to compare effects of increasing chronic stress load on the stress response of European seabass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata) to identify neuroendocrine functions that regulate this response. Fish were left undisturbed (controls) or exposed to three levels of chronic stress for 3 weeks and then subjected to an acute stress test (ACT). Chronic stress impeded growth and decreased feed consumption in seabass, not in seabream. In seabass basal cortisol levels are high and increase with stress load; the response to a subsequent ACT decreases with increasing (earlier) load. Basal cortisol levels in seabream increase with the stress load, whereas the ACT induced a similar response in all groups. In seabass and seabream plasma α-MSH levels and brain stem serotonergic activity and turnover were similar and not affected by chronic stress. Species-specific molecular neuro-regional differences were seen. In-situ hybridization analysis of the early immediate gene cfos in the preoptic area showed ACT-activation in seabream; in seabass the expression level was not affected by ACT and seems constitutively high. In seabream, expression levels of telencephalic crf, crfbp, gr1, and mr were downregulated; the seabass hypothalamic preoptic area showed increased expression of crf and gr1, and decreased expression of mr, and this increased the gr1/mr ratio considerably. We substantiate species-specific physiological differences to stress coping between seabream and seabass at an endocrine and neuroendocrine molecular level. Seabass appear less resilient to stress, which we conclude from high basal activities of stress-related parameters and poor, or absent, responses to ACT. This comparative study reveals important aquaculture, husbandry, and welfare implications for the rearing of these species.