Comparative analysis of a teleost skeleton transcriptome provides insight into its regulation

Estradiol-17β and bisphenol A affect growth and mineralization in early life stages of seabass

Natural and synthetic estrogens are contaminants present in aquatic ecosystems. They can have significant consequences on the estrogen-sensitive functions of organisms, including skeletal development and growth of vertebrate larvae. Synthetic polyphenols represent a group of environmental xenoestrogens capable of binding the receptors for the natural hormone estradiol-17β (E2). To better understand how (xeno-)estrogens can affect the skeleton in fish species with high ecological and commercial interest, 16 days post-hatch larvae of the seabass were experimentally exposed for 7 days to E2 and Bisphenol A (BPA), both used at the regulatory concentration of surface water quality (E2: 0.4 ng.L-1, BPA: 1.6 μg.L-1) or at a concentration 100 times higher. Skeletal mineralization levels were evaluated using Alizarin red staining, and expression of several genes playing key roles in growth, skeletogenesis and estrogen signaling pathways was assessed by qPCR. Our results show that E2 exerts an overall negative effect on skeletal mineralization at the environmental concentration of 0.4 ng.L-1, correlated with an increase in the expression of genes associated only with osteoblast bone cells. Both BPA exposures inhibited mineralization with less severe effects and modified bone homeostasis by regulating the expression of gene encoding osteoblasts and osteoclasts markers. Our results demonstrate that environmental E2 exposure inhibits larval growth and has an additional inhibitory effect on skeleton mineralization while both BPA exposures have marginal inhibitory effect on skeletal mineralization. All exposures have significant effects on transcriptional levels of genes involved in the skeletal development of seabass larvae.

A comprehensive coding and microRNA transcriptome of vertebral bone in postlarvae and juveniles of Senegalese sole (Solea senegalensis)

Understanding vertebral bone development is essential to prevent skeletal malformations in farmed fish related to genetic and environmental factors. This is an important issue in Solea senegalensis, with special impact of spinal anomalies in postlarval and juvenile stages. Vertebral bone transcriptomics in farmed fish mainly comes from coding genes, and barely on miRNA expression. Here, we used RNA-seq of spinal samples to obtain the first comprehensive coding and miRNA transcriptomic repertoire for postlarval and juvenile vertebral bone, covering different vertebral phenotypes and egg-incubation temperatures related to skeleton health in S. senegalensis. Coding genes, miRNA and pathways regulating bone development and growth were identified. Differential transcriptomic profiles and suggestive mRNA-miRNA interactions were found between postlarvae and juveniles. Bone-related genes and functions were associated with the extracellular matrix, development and regulatory processes, calcium binding, retinol and lipid metabolism or response to stimulus, including those revealed by the miRNA targets related to signaling, cellular and metabolic processes, growth, cell proliferation and biological adhesion. Pathway enrichment associated with fish skeleton were identified when comparing postlarvae and juveniles: growth and bone development functions in postlarvae, while actin cytoskeleton, focal adhesion and proteasome related to bone remodeling in juveniles. The transcriptome data disclosed candidate coding and miRNA gene markers related to bone cell processes, references for functional studies of the anosteocytic bone of S. senegalensis. This study establishes a broad transcriptomic foundation to study healthy and anomalous spines under early thermal conditions across life-stages in S. senegalensis, and for comparative analysis of skeleton homeostasis and pathology in fish and vertebrates.

Exploring Omega-3's Impact on the Expression of Bone-Related Genes in Meagre (Argyrosomus regius)

Dietary supplementation with Omega-3 fatty acids seems to promote skeletal health. Therefore, their consumption at imbalanced or excessive levels has offered less beneficial or even prejudicial effects. Fish produced in aquaculture regimes are prone to develop abnormal skeletons. Although larval cultures are usually fed with diets supplemented with Omega-3 Long Chain Polyunsaturated fatty acids (LC-PUFAs), the lack of knowledge about the optimal requirements for fatty acids or about their impact on mechanisms that regulate skeletal development has impeded the design of diets that could improve bone formation during larval stages when the majority of skeletal anomalies appear. In this study, Argyrosomus regius larvae were fed different levels of Omega-3s (2.6% and 3.6% DW on diet) compared to a commercial diet. At 28 days after hatching (DAH), their transcriptomes were analyzed to study the modulation exerted in gene expression dynamics during larval development and identify impacted genes that can contribute to skeletal formation. Mainly, both levels of supplementation modulated bone-cell proliferation, the synthesis of bone components such as the extracellular matrix, and molecules involved in the interaction and signaling between bone components or in important cellular processes. The 2.6% level impacted several genes related to cartilage development, denoting a special impact on endochondral ossification, delaying this process. However, the 3.6% level seemed to accelerate this process by enhancing skeletal development. These results offered important insights into the impact of dietary Omega-3 LC-PUFAs on genes involved in the main molecular mechanism and cellular processes involved in skeletal development.

Skeletal Morphogenesis and Anomalies in Gilthead Seabream: A Comprehensive Review

The gilthead seabream, one of the most important species in Mediterranean aquaculture, with an increasing status of exploitation in terms of production volume and aquafarming technologies, has become an important research topic over the years. The accumulation of knowledge from several studies conducted during recent decades on their functional and biological characteristics has significantly improved their aquacultural aspects, namely their reproductive success, survival, and growth. Despite the remarkable progress in the aquaculture industry, hatchery conditions are still far from ideal, resulting in frequent abnormalities at the beginning of intensive culture, entailing significant economic losses. Those deformities are induced during the embryonic and post-embryonic periods of life, and their development is still poorly understood. In the present review, we created a comprehensive synthesis that covers the various aspects of skeletal morphogenesis and anomalies in the gilthead seabream, highlighting the genetic, environmental, and nutritional factors contributing to bone deformities and emphasized the potential of the gilthead seabream as a model organism for understanding bone morphogenesis in both aquaculture and translational biological research. This review article addresses the existing lack in the literature regarding gilthead seabream bone deformities, as there are currently no comprehensive reviews on this subject.

Interaction between Dietary Vitamin D3 and Vitamin K3 in Gilthead Seabream Larvae (Sparus aurata) in Relation to Growth and Expression of Bone Development-Related Genes

Vitamins D and K are essential fat-soluble nutrients that intervene in bone development processes among other biological functions. The present study is aimed at investigating the potential combined effect of dietary supplementation with vitamin D₃ (cholecalciferol) and vitamin K₃ (menadione) in gilthead seabream (Sparus aurata) larvae. For that purpose, seabream diets were supplemented with different combinations of vitamin D₃/vitamin K₃ (mg/kg diet) as follows: 0.00/0, 0.06/70, 0.06/170, 0.13/70, 0.13/170, 0.40/70, and 0.40/170. Feeding gilthead seabream larvae (22 days post hatch) for 21 days with the diets supplemented with 0.06-0.13 mg/kg vitamin D₃ and 70 mg/kg vitamin K₃ (diets 0.06/70 and 0.13/70) led to the highest larval growth and survival and the highest expression of important biomarkers of both bone development and health, such as bmp2, osx, and mgp, and calcium homeostasis, such as pthrp and casr. However, the increased supplementation with both vitamins at 0.40 mg/kg vitamin D₃ and 170 mg/kg vitamin K₃ (diet 0.40/170) reduced larval growth and survival, downregulated bmp2 and pthrp expressions, and upregulated osx and mgp, causing an unbalance in the relative expression of these genes. The results of the present study have shown the interaction between vitamin D₃ supplementation and vitamin K₃ supplementation in larval performance and gene expression related to bone development and calcium homeostasis, denoting the significance of a correct balance between both vitamins in larval diets.

Muscle regeneration in gilthead sea bream: Implications of endocrine and local regulatory factors and the crosstalk with bone

Fish muscle regeneration is still a poorly known process. In the present study, an injury was done into the left anterior epaxial skeletal muscle of seventy 15 g gilthead sea bream (Sparus aurata) juveniles to evaluate at days 0, 1, 2, 4, 8, 16 and 30 post-wound, the expression of several muscle genes. Moreover, transcripts' expression in the bone (uninjured tissue) was also analyzed. Histology of the muscle showed the presence of dead tissue the first day after injury and how the damaged fibers were removed and replaced by new muscle fibers by day 16 that kept growing up to day 30. Gene expression results showed in muscle an early upregulation of igf-2 and a downregulation of ghr-1 and igf-1. Proteolytic systems expression increased with capn2 and ctsl peaking at 1 and 2 days post-injury, respectively and mafbx at day 8. A pattern of expression that fitted well with active myogenesis progression 16 days after the injury was then observed, with the recovery of igf-1, pax7, cmet, and cav1 expression; and later on, that of cav3 as well. Furthermore, the first days post-injury, the cytokines il-6 and il-15 were also upregulated confirming the tissue inflammation, while tnfα was only upregulated at days 16 and 30 to induce satellite cells recruitment; overall suggesting a possible role for these molecules as myokines. The results of the bone transcripts showed an upregulation first, of bmp2 and ctsk at days 1 and 2, respectively; then, ogn1 and ocn peaked at day 4 in parallel to mstn2 downregulation, and runx2 and ogn2 increased after 8 days of muscle injury, suggesting a possible tissue crosstalk during the regenerative process. Overall, the present model allows studying the sequential involvement of different regulatory molecules during muscle regeneration, as well as the potential relationship between muscle and other tissues such as bone to control musculoskeletal development and growth, pointing out an interesting new line of research in this group of vertebrates.

Hydroxytyrosol-rich extract from olive juice as an additive in gilthead sea bream juveniles fed a high-fat diet: Regulation of somatic growth

The dietary inclusion of plant-based products in fish feeds formulation is required for the sustainable development of aquaculture. Moreover, considering functional diets, hydroxytyrosol, one of the major phenolic compounds found in olives (Olea europaea), has been identified as a potential candidate to be used in the aquafeeds industry due to its health promoting abilities. The aim of this study was to evaluate the effects of the inclusion of an olive juice extract rich in hydroxytyrosol as an additive (0.52 g HT/kg feed) in a high-fat (24% lipids) diet in gilthead sea bream (Sparus aurata) juveniles. Moreover, the experimental diets, with or without the extract, were administered daily at a standard (3% of total biomass in the tank) or restricted ration (40% reduction) for 8–9 weeks. Growth and biometric parameters, insulin-like growth factor 1 (IGF-1) plasma levels and growth hormone/IGF axis-, myogenic- and osteogenic-related genes expression in liver, white muscle and/or bone were analyzed. Moreover, in vitro cultures of vertebra bone-derived cells from fish fed the diets at a standard ration were performed at weeks 3 and 9 to explore the effects of hydroxytyrosol on osteoblasts development. Although neither body weight or any other biometric parameter were affected by diet composition after 4 or 8 weeks, the addition of the hydroxytyrosol-rich extract to the diet increased IGF-1 plasma levels, regardless of the ration regime, suggesting an anabolic condition. In muscle, the higher mRNA levels of the binding protein igfbp-5b and the myoblast fusion marker dock5 in fish fed with the hydroxytyrosol-rich diet suggested that this compound may have a role in muscle, inducing development and a better muscular condition. Furthermore in bone, increased osteogenic potential while delayed matrix mineralization after addition to the diet of the olive juice extract was supported by the upregulated expression of igf-1 and bmp4 and reduced transcript levels of osteopontin. Overall, this study provides new insights into the beneficial use of hydroxytyrosol as a dietary additive in gilthead sea bream functional diets to improve muscle-skeletal condition and, the aquaculture industry.

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple Summary

The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.

Abstract

The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F₁ juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

Comparative transcriptome analyses and identification of candidate genes involved in vertebral abnormality of bighead carp Hypophthalmichthys nobilis

Body deformity occurs both in wild and farmed fishes, which is one of the most challenging problems for aquaculture industry. In most cases, such body deformities are linked to skeletal deformities. Currently, very limited information is available on skeletal deformities of farmed fish species which may be caused by genetic factor. In this study, we performed muscle and vertebra transcriptome analyses in body deformity and normality of bighead carp Hypophthalmichthys nobilis (from one meiotic gynogenesis family) using RNA-Seq. A total of 43,923 and 44,416 unigenes were predicted in muscles and vertebrae, respectively. Based on these data, we further explored the gene expression profiles in gynogenetic normal and abnormal bighead carp. No differentially expressed gene (DEG) was found in transcriptome data of muscles. Totally, 20 key DEGs were identified in transcriptome data of vertebrae, such as low density lipoprotein-related protein 2 (lrp2), bone morphogenetic protein 2B (bmp2b) and collagen alpha-1(IV) (col4a1). 12 potential pathways were also identified in vertebra transcriptome data, which were mainly involved in development, growth, cytoskeleton and energy metabolism, such as MAPK signaling pathway, regulation of actin cytoskeleton and TGF-beta signaling pathway. Results of this study will be informative for the understanding of genetic mechanisms for body shape formation and also provide potential candidate genes for selection program involved in body shape and skeletal development in H. nobilis.

Spatiotemporal expression of activin receptor-like kinase-5 and bone morphogenetic protein receptor type II in the ovary of shortfinned eel, Anguilla australis

In the eel ovary, the expression of growth differentiation factor-9 (Gdf9) appears to be largely confined to the germ cell in early stages of oogenesis. However, both the target tissue and the function of Gdf9 in fish remain unknown. This study aimed to describe the abundance and localization of activin receptor-like kinase-5 (Alk5) and bone morphogenetic protein receptor type II (Bmpr2), which together mediate the Gdf9 signal, in the ovary of a basal teleost, the shortfinned eel, Anguilla australis, during early folliculogenesis. The cDNA encoding eel alk5 and bmpr2 genes were cloned, characterized and the transcript abundances of these receptors quantified by quantitative real-time PCR. Ovarian transcript abundance for both receptors, along with that of gdf9 and of its paralogue bmp15, increased from the previtellogenic to early vitellogenic stage. Localization of receptor mRNAs by in situ hybridization revealed that these receptors are located in the somatic cells surrounding the oocyte. Furthermore, tissue distribution analysis showed that the expression of alk5 and bmpr2 were highest in ovary and thyroid, respectively. Unexpectedly, however, bmpr2 mRNA levels were lower in the ovary than in any of the other 17 tissues examined, and indeed, lower than ovarian gdf9 transcript abundance. These findings, together with the ovarian expression pattern of Gdf9, suggest that Gdf9, and conceivably, Bmp15, from the oocyte can signal through receptors that are located on the somatic cells surrounding the oocyte; this, in turn, facilitates elucidation of the function of these growth factors during oogenesis in teleost fish.

Regulatory mechanisms involved in muscle and bone remodeling during refeeding in gilthead sea bream

The tolerance of fish to fasting offers a model to study the regulatory mechanisms and changes produced when feeding is restored. Gilthead sea bream juveniles were exposed to a 21-days fasting period followed by 2 h to 7-days refeeding. Fasting provoked a decrease in body weight, somatic indexes, and muscle gene expression of members of the Gh/Igf system, signaling molecules (akt, tor and downstream effectors), proliferation marker pcna, myogenic regulatory factors, myostatin, and proteolytic molecules such as cathepsins or calpains, while most ubiquitin-proteasome system members increased or remained stable. In bone, downregulated expression of Gh/Igf members and osteogenic factors was observed, whereas expression of the osteoclastic marker ctsk was increased. Refeeding recovered the expression of Gh/Igf system, myogenic and osteogenic factors in a sequence similar to that of development. Akt and Tor phosphorylation raised at 2 and 5 h post-refeeding, much faster than its gene expression increased, which occurred at day 7. The expression in bone and muscle of the inhibitor myostatin (mstn2) showed an inverse profile suggesting an inter-organ coordination that needs to be further explored in fish. Overall, this study provides new information on the molecules involved in the musculoskeletal system remodeling during the early stages of refeeding in fish.

Multiomics analyses reveal that NOD-like signaling pathway plays an important role against Streptococcus agalactiae in the spleen of tilapia

Streptococcus aglactiae(GBS) infection in tilapia is a serious global disease that causes significant production loss. Here, we studied the role of GBS in the spleen and the spleen's response against the pathogen through dual RNA-seq and proteome technology. Animals were divided into three groups: control, virulent treated (HN016), and attenuated treated (YM001). Spleen samples were collected and analysis when a disease outbreak. Dual RNA-seq result showed the virulence factor genes of GBS, included CAMP factor, PGK, OCT, enolase, scpB, Sip, bca, were upregulation. downregulation of GapA, cylE, OCT, scpB, C5AP, rlmB, hly, FBP, in HN016 and YM001. But for proteomic, OCT and bca were downregulation, the others were upregulation. For host transcriptome KEGG analysis showed, the NOD-like receptor signaling pathway (NLRs) and TOLL-like receptor signaling pathway (TLRs) were upreguoation in HN016 infected fish than the control fish; But for proteome KEGG, only the NLRS was up, the TLRS was not change. Compared with YM001 infected fishes, for transcriptome, NLRs and TLRs in infected HN016 fishes were significance rise (p < 0.01); for proteome, the NLRs was up (p < 0.05), but TLRs was no change.Analysis of pathogen-host interaction showed that the peptidoglycan (PNG), CD2, LCK, and host's Zap70 were involved in the regulation of NLRs; PNG, LCK, and ZAP70 were involved in the regulation of TRLs. Conclusion: the virulent strain HN016 and attenuated strainYM001 differed in the quantity of virulence factors. In tilapia's innate immune system, NLRs was the main defense factors, but bacteria avoided the host defense through TLRs.

Development of Teleost Intermuscular Bones Undergoing Intramembranous Ossification Based on Histological-Transcriptomic-Proteomic Data

Intermuscular bones (IBs) specially exist in lower teleost fish and the molecular mechanism for its development remains to be clarified. In this study, different staining methods and comparative proteomics were conducted to investigate the histological structure and proteome of IB development in Megalobrama amblycephala, including four key IB developmental stages (S1-IBs have not emerged in the tail part; S2-several small IBs started to ossify in the tail part; S3-IBs appeared rapidly; S4-all the IBs appeared with mature morphology). Alcian blue and alizarin red S stained results indicated that IBs were gradually formed from S2 to S4, undergoing intramembranous ossification without a cartilaginous phase. A total of 3368 proteins were identified by using the isobaric tags for relative and absolute quantitation (iTRAQ) approach. Functional annotation showed that proteins which were differentially expressed among stages were involved in calcium, MAPK, Wnt, TGF-β, and osteoclast pathways which played a critical role in bone formation and differentiation. Three proteins (collagen9α1, stat1, tnc) associated with chondrocytes did not exhibit significant changes through S2 to S4; however, proteins (entpd5, casq1a, pvalb, anxa2a, anxa5) which associated with osteoblasts and bone formation and differentiation showed significantly a higher expression level from S1 to S2, as well as to S3 and S4. These further demonstrated that development of IBs did not go through a cartilaginous phase. The inhibitors of TGF-β and Wnt pathways were tested on zebrafish (sp7/eGFP) and the results indicated that both inhibitors significantly delayed IB development. This study provides a comprehensive understanding of the IB ossification pattern, which will help further elucidate the molecular mechanisms for IB development in teleosts.

Gene expression analyses in malformed skeletal structures of gilthead sea bream (Sparus aurata)

The incidence of skeletal anomalies in reared fish has been translated for years in important economic losses for the aquaculture industry. In the present study, we have analysed the gene expression of extracellular matrix components and transcription factors involved in bone development in gilthead sea bream presenting different skeletal anomalies: lordosis (LD), lordosis-scoliosis-kyphosis (LSK) or opercular, dental or jaw malformations in comparison with control (CT) specimens. Results showed a possible link between the presence of LD and LSK and the significant downregulation of genes involved in osteoblasts' maturation and matrix mineralization (collagen type 1-alpha, osteopontin, osteocalcin, matrix Gla protein and tissue non-specific alkaline phosphatase), as well as in bone resorption (cathepsin K and matrix metalloproteinase 9) compared to CT animals. Contrarily, the key osteogenic transcription factor runx2 was upregulated in the malformed vertebra suggesting impaired determination of mesenchymal stem cells towards the osteoblastic lineage. Despite the gene expression patterns of the other malformed structures were not affected in comparison with CT fish, the results of the present study may contribute in the long term to identify potential candidate gene profiles associated with column deformities that may help reducing the incidence of appearance of skeletal anomalies in this important aquaculture species.

Fatty acids from fish or vegetable oils promote the adipogenic fate of mesenchymal stem cells derived from gilthead sea bream bone potentially through different pathways

Fish are rich in n-3 long-chain polyunsaturated fatty acids (LC-PUFA), such as eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids, thus they have a great nutritional value for human health. In this study, the adipogenic potential of fatty acids commonly found in fish oil (EPA and DHA) and vegetable oils (linoleic (LA) and alpha-linolenic (ALA) acids), was evaluated in bone-derived mesenchymal stem cells (MSCs) from gilthead sea bream. At a morphological level, cells adopted a round shape upon all treatments, losing their fibroblastic form and increasing lipid accumulation, especially in the presence of the n-6 PUFA, LA. The mRNA levels of the key transcription factor of osteogenesis, runx2 significantly diminished and those of relevant osteogenic genes remained stable after incubation with all fatty acids, suggesting that the osteogenic process might be compromised. On the other hand, transcript levels of the main adipogenesis-inducer factor, pparg increased in response to EPA. Nevertheless, the specific PPARγ antagonist T0070907 appeared to suppress the effects being caused by EPA over adipogenesis. Moreover, LA, ALA and their combinations, significantly up-regulated the fatty acid transporter and binding protein, fatp1 and fabp11, supporting the elevated lipid content found in the cells treated with those fatty acids. Overall, this study has demonstrated that fatty acids favor lipid storage in gilthead sea bream bone-derived MSCs inducing their fate into the adipogenic versus the osteogenic lineage. This process seems to be promoted via different pathways depending on the fatty acid source, being vegetable oils-derived fatty acids more prone to induce unhealthier metabolic phenotypes.

Vertebrate SLRP family evolution and the subfunctionalization of osteoglycin gene duplicates in teleost fish

Background

Osteoglycin (OGN, a.k.a. mimecan) belongs to cluster III of the small leucine-rich proteoglycans (SLRP) of the extracellular matrix (ECM). In vertebrates OGN is a characteristic ECM protein of bone. In the present study we explore the evolution of SLRP III and OGN in teleosts that have a skeleton adapted to an aquatic environment.

Results

The SLRP gene family has been conserved since the separation of chondrichthyes and osteichthyes. Few gene duplicates of the SLRP III family exist even in the teleosts that experienced a specific whole genome duplication. One exception is ogn for which duplicate copies were identified in fish genomes. The ogn promoter sequence and in vitro mesenchymal stem cell (MSC) cultures suggest the duplicate ogn genes acquired divergent functions. In gilthead sea bream (Sparus aurata) ogn1 was up-regulated during osteoblast and myocyte differentiation in vitro, while ogn2 was severely down-regulated during bone-derived MSCs differentiation into adipocytes in vitro.

Conclusions

Overall, the phylogenetic analysis indicates that the SLRP III family in vertebrates has been under conservative evolutionary pressure. The retention of the ogn gene duplicates in teleosts was linked with the acquisition of different functions. The acquisition by OGN of functions other than that of a bone ECM protein occurred early in the vertebrate lineage.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1310-2) contains supplementary material, which is available to authorized users.

Temperature responsiveness of gilthead sea bream bone; an in vitro and in vivo approach

This study aimed to characterize the molecules involved in osteogenesis in seabream and establish using in vitro/in vivo approaches the responsiveness of selected key genes to temperature. The impact of a temperature drop from 23 to 13 °C was evaluated in juvenile fish thermally imprinted during embryogenesis. Both, in vitro/in vivo, Fib1a, appeared important in the first stages of bone formation, and Col1A1, ON and OP, in regulating matrix production and mineralization. OCN mRNA levels were up-regulated in the final larval stages when mineralization was more intense. Moreover, temperature-dependent differential gene expression was observed, with lower transcript levels in the larvae at 18 °C relative to those at 22 °C, suggesting bone formation was enhanced in the latter group. Results revealed that thermal imprinting affected the long-term regulation of osteogenesis. Specifically, juveniles under the low and low-to-high-temperature regimes had reduced levels of OCN when challenged, indicative of impaired bone development. In contrast, gene expression in fish from the high and high-to-low-temperature treatments was unchanged, suggesting imprinting may have a protective effect. Overall, the present study revealed that thermal imprinting modulates bone development in seabream larvae, and demonstrated the utility of the in vitro MSC culture as a reliable tool to investigate fish osteogenesis.

Matricellular Proteins: Functional Insights From Non-mammalian Animal Models

The extracellular matrix (ECM) has central roles in tissue integrity and remodeling throughout the life span of animals. While collagens are the most abundant structural components of ECM in most tissues, tissue-specific molecular complexity is contributed by ECM glycoproteins. The matricellular glycoproteins are categorized primarily according to functional criteria and represented predominantly by the thrombospondin, tenascin, SPARC/osteonectin, and CCN families. These proteins do not self-assemble into ECM fibrils; nevertheless, they shape ECM properties through interactions with structural ECM proteins, growth factors, and cells. Matricellular proteins also promote cell migration or morphological changes through adhesion-modulating or counter-adhesive actions on cell-ECM adhesions, intracellular signaling, and the actin cytoskeleton. Typically, matricellular proteins are most highly expressed during embryonic development. In adult tissues, expression is more limited unless activated by cues for dynamic tissue remodeling and cell motility, such as occur during inflammatory response and wound repair. Many insights in the complex roles of matricellular proteins have been obtained from studies of gene knockout mice. However, with the exception of chordate-specific tenascins, these are highly conserved proteins that are encoded in many animal phyla. This review will consider the increasing body of research on matricellular proteins in nonmammalian animal models. These models provide better access to the very earliest stages of embryonic development and opportunities to study biological processes such as limb and organ regeneration. In aggregate, this research is expanding concepts of the functions and mechanisms of action of matricellular proteins.

A transcriptomic approach to study the effect of long-term starvation and diet composition on the expression of mitochondrial oxidative phosphorylation genes in gilthead sea bream (Sparus aurata)

Background

The impact of nutritional status and diet composition on mitochondrial oxidative phosphorylation (OXPHOS) in fish remains largely unknown. To identify biomarkers of interest in nutritional studies, herein we obtained a deep-coverage transcriptome by 454 pyrosequencing of liver and skeletal muscle cDNA normalised libraries from long-term starved gilthead sea bream (Sparus aurata) and fish fed different diets.

Results

After clean-up of high-throughput deep sequencing reads, 699,991 and 555,031 high-quality reads allowed de novo assembly of liver and skeletal muscle sequences, respectively (average length: 374 and 441 bp; total megabases: 262 and 245 Mbp). An additional incremental assembly was completed by integrating data from both tissues (hybrid assembly). Assembly of hybrid, liver and skeletal muscle transcriptomes yielded, respectively, 19,530, 11,545 and 10,599 isotigs (average length: 1330, 1208 and 1390 bp, respectively) that were grouped into 15,954, 10,033 and 9189 isogroups. Following annotation, hybrid transcriptomic data were used to construct an oligonucleotide microarray to analyse nutritional regulation of the expression of 129 genes involved in OXPHOS in S. aurata. Starvation upregulated cytochrome c oxidase components and other key OXPHOS genes in the liver, which exhibited higher sensitive to food deprivation than the skeletal muscle. However, diet composition affected OXPHOS in the skeletal muscle to a greater extent than in the liver: most of genes upregulated under starvation presented higher expression among fish fed a high carbohydrate/low protein diet.

Conclusions

Our findings indicate that the expression of coenzyme Q-binding protein (COQ10), cytochrome c oxidase subunit 6A2 (COX6A2) and ADP/ATP translocase 3 (SLC25A6) in the liver, and cytochrome c oxidase subunit 5B isoform 1 (COX5B1) in the liver and the skeletal muscle, are sensitive markers of the nutritional condition that may be relevant to assess the effect of changes in the feeding regime and diet composition on fish farming.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4148-x) contains supplementary material, which is available to authorized users.

Comparative proteomics analysis of teleost intermuscular bones and ribs provides insight into their development

BACKGROUND

Intermuscular bones (IBs) and ribs both are a part of skeletal system in teleosts, but with different developing process. The chemical composition of fish IBs and ribs as well as the underlying mechanism about their development have not been investigated. In the present study, histological structures showed that one bone cavity containing osteoclasts were existed in ribs, but not in IBs of Megalobrama amblycephala. We constructed the first proteomics map for fish bones including IBs and ribs, and identified the differentially expressed proteins between IBs and ribs through iTRAQ LC-MS/MS proteomic analysis.

RESULTS

The proteins extracted from IBs and ribs at 1- to 2-year old M. amblycephala were quantified 2,342 proteins, with 1,451 proteins annotated with GO annotation in biological processes, molecular function and cellular component. A number of bone related proteins as well as pathways were identified in the study. A total of 93 and 154 differently expressed proteins were identified in comparison groups of 1-IB-vs-1-Rib and 2-IB-vs-2-Rib, which indicated the obvious differences of chemical composition between these two bone tissues. The two proteins (vitronectin b precursor and matrix metalloproteinase-2) related to osteoclasts differentiation were significantly up-regulated in ribs compared with IBs (P < 0.05), which was in accordance with the results from histological structures. In comparison groups of 1-IB-vs-2-IB and 1-Rib-vs-2-Rib, 33 and 51 differently expressed proteins were identified and the function annotation results showed that these proteins were involved in regulating bone development and differentiation. Subsequently, 11 and 13 candidate proteins in comparison group of 1-IB-vs-1-Rib and 1-IB-vs-2-IB related to bone development were validated by MRM assays.

CONCLUSIONS

Our present study suggested the different key proteins involved in the composition of fish ribs and IBs as well as their growth development. These findings could provide important clues towards further understanding of fish skeletal system and the roles of proteins playing in regulating diverse biological processes in fish.

Thermal imprinting modifies bone homeostasis in cold challenged sea bream (Sparus aurata, L.)

Fish are ectotherms and temperature plays a determinant role in their physiology, biology and ecology and is a driver of seasonal responses. The present study assessed how thermal imprinting during embryonic and larval stages modified the response of adult fish to low water temperature. We targeted the gilthead sea bream that develops a condition known as winter syndrome when it is exposed to low water temperatures. Eggs and larvae of sea bream were exposed to four different thermal regimes and then the response of the resulting adults to a low temperature challenge was assessed. Sea bream exposed to a high-low thermal regime as eggs and larvae (HLT, 22°C until hatch and then 18°C until larvae-juvenile transition) had increased plasma cortisol and lower sodium and potassium in response to a cold challenge compared to the other thermal history groups. Plasma glucose and osmolality were increased in cold challenge HLT fish relative to the unchallenged HLT fish. Cold challenge modified bone homeostasis/responsiveness in the low-high thermal regime group (LHT) relative to other groups and ocn, ogn1/2, igf1, gr and trα/β transcripts were all down-regulated. In the low temperature group (LT) and HLT group challenged with a low temperature, ALP/TRAP activities were decreased relative to unchallenged groups and bone calcium content also decreased in the LT group. Overall, the results indicate that thermal imprinting during early development of sea bream causes a change in the physiological response of adults to a cold challenge.

Having a BLAST: Searchable transcriptome resources for the gilthead sea bream and the European sea bass

The gilthead sea bream (Sparus aurata) and the European sea bass (Dicentrarchus labrax) are the most important aquaculture species in the Mediterranean Sea and since the last decade it has been seen an exponential increase in their available molecular resources. In order to improve accessibility to transcriptome resources, Expressed Sequence Tags (ESTs), mRNA sequences and raw read sequences were assembled and deposited in BLAST queryable databases. The publicly available sea bream and sea bass sequences (6.4 and 247.5 million) generated 45,094 and 68,117 assembled sequences, with, respectively, arithmetic mean size of 998 and 2125bp and N50 of 1302 and 2966bp. The assemblies will be regularly updated and new analytical tools added to the web server at http://sea.ccmar.ualg.pt.

Dynamic mRNA and miRNA expression analysis in response to intermuscular bone development of blunt snout bream (Megalobrama amblycephala)

Intermuscular bone (IB), which occurs only in the myosepta of lower teleosts, is attracting more attention because they are difficult to remove and make the fish unpleasant to eat. By gaining a better understanding of the genetic regulation of IB development, an integrated analysis of miRNAs and mRNAs expression profiling was performed on Megalobrama amblycephala. Four key development stages were selected for transcriptome and small RNA sequencing. A number of significantly differentially expressed miRNAs/genes associated with bone formation and differentiation were identified and the functional characteristics of these miRNAs/genes were revealed by GO function and KEGG pathway analysis. These were involved in TGF-β, ERK and osteoclast differentiation pathways known in the literature to affect bone formation and differentiation. MiRNA-mRNA interaction pairs were detected from comparison of expression between different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating bone development and differentiation. A negative regulation effect of two miRNAs was verified through dual luciferase reporter assay. As a unique public resource for gene expression and regulation during the IB development, this study is expected to provide forwards ideas and resources for further biological researches to understand the IBs’ development.

Effects of high temperatures on threatened estuarine fishes during periods of extreme drought

Climate change and associated increases in water temperatures may impact physiological performance in ectotherms and exacerbate endangered species declines. We used an integrative approach to assess the impact of elevated water temperature on two fishes of immediate conservation concern in a large estuary system, the threatened longfin smelt (Spirinchus thaleichthys) and endangered delta smelt (Hypomesus transpacificus). Abundances have reached record lows in California, USA, and these populations are at imminent risk of extirpation. California is currently impacted by a severe drought, resulting in high water temperatures, conditions that will become more common as a result of climate change. We exposed fish to environmentally relevant temperatures (14°C and 20°C) and used RNA sequencing to examine the transcriptome-wide responses to elevated water temperature in both species. Consistent with having a lower temperature tolerance, longfin smelt exhibited a pronounced cellular stress response, with an upregulation of heat shock proteins, after exposure to 20°C that was not observed in delta smelt. We detected an increase in metabolic rate in delta smelt at 20°C and increased expression of genes involved in metabolic processes and protein synthesis, patterns not observed in longfin smelt. Through examination of responses across multiple levels of biological organization, and by linking these responses to habitat distributions in the wild, we demonstrate that longfin smelt may be more susceptible than delta smelt to increases in temperatures, and they have little room to tolerate future warming in California. Understanding the species-specific physiological responses of sensitive species to environmental stressors is crucial for conservation efforts and managing aquatic systems globally.

Chronic exposures to low and high concentrations of ibuprofen elicit different gene response patterns in a euryhaline fish

Ibuprofen is one of the most commonly detected pharmaceuticals in wastewater effluent; however, the effects of ibuprofen on aquatic organisms are poorly understood. This study presents the transcriptome-wide response of the inland silverside, Menidia beryllina, to chronic exposure to ibuprofen. At the lowest exposure concentration (0.0115 mg/L), we detected a downregulation of many genes involved in skeletal development, aerobic respiration, and immune function. At the highest exposure concentration (1.15 mg/L), we detected increased expression of regulatory genes in the arachidonic acid metabolism pathway and several immune genes involved in an inflammatory response. Additionally, there was differential expression of genes involved in oxidative stress responses and a downregulation of genes involved in osmoregulation. This study provides useful information for monitoring the effects of this common wastewater effluent contaminant in the environment and for the generation of biomarkers of exposure to ibuprofen that may be transferable to other fish species.

On the evolutionary relationship between chondrocytes and osteoblasts

Vertebrates are the only animals that produce bone, but the molecular genetic basis for this evolutionary novelty remains obscure. Here, we synthesize information from traditional evolutionary and modern molecular genetic studies in order to generate a working hypothesis on the evolution of the gene regulatory network (GRN) underlying bone formation. Since transcription factors are often core components of GRNs (i.e., kernels), we focus our analyses on Sox9 and Runx2. Our argument centers on three skeletal tissues that comprise the majority of the vertebrate skeleton: immature cartilage, mature cartilage, and bone. Immature cartilage is produced during early stages of cartilage differentiation and can persist into adulthood, whereas mature cartilage undergoes additional stages of differentiation, including hypertrophy and mineralization. Functionally, histologically, and embryologically, these three skeletal tissues are very similar, yet unique, suggesting that one might have evolved from another. Traditional studies of the fossil record, comparative anatomy and embryology demonstrate clearly that immature cartilage evolved before mature cartilage or bone. Modern molecular approaches show that the GRNs regulating differentiation of these three skeletal cell fates are similar, yet unique, just like the functional and histological features of the tissues themselves. Intriguingly, the Sox9 GRN driving cartilage formation appears to be dominant to the Runx2 GRN of bone. Emphasizing an embryological and evolutionary transcriptomic view, we hypothesize that the Runx2 GRN underlying bone formation was co-opted from mature cartilage. We discuss how modern molecular genetic experiments, such as comparative transcriptomics, can test this hypothesis directly, meanwhile permitting levels of constraint and adaptation to be evaluated quantitatively. Therefore, comparative transcriptomics may revolutionize understanding of not only the clade-specific evolution of skeletal cells, but also the generation of evolutionary novelties, providing a modern paradigm for the evolutionary process.

Differential expression of the aryl hydrocarbon receptor pathway associates with craniofacial polymorphism in sympatric Arctic charr

Background

The developmental basis of craniofacial morphology hinges on interactions of numerous signalling systems. Extensive craniofacial variation in the polymorphic Arctic charr, a member of the salmonid family, from Lake Thingvallavatn (Iceland), offers opportunities to find and study such signalling pathways and their key regulators, thereby shedding light on the developmental pathways, and the genetics of trophic divergence.

Results

To identify genes involved in the craniofacial differences between benthic and limnetic Arctic charr, we used transcriptome data from different morphs, spanning early development, together with data on craniofacial expression patterns and skeletogenesis in model vertebrate species. Out of 20 genes identified, 7 showed lower gene expression in benthic than in limnetic charr morphs. We had previously identified a conserved gene network involved in extracellular matrix (ECM) organization and skeletogenesis, showing higher expression in developing craniofacial elements of benthic than in limnetic Arctic charr morphs. The present study adds a second set of genes constituting an expanded gene network with strong, benthic–limnetic differential expression. To identify putative upstream regulators, we performed knowledge-based motif enrichment analyses on the regulatory sequences of the identified genes which yielded potential binding sites for a set of known transcription factors (TFs). Of the 8 TFs that we examined using qPCR, two (Ahr2b and Ap2) were found to be differentially expressed between benthic and limnetic charr. Expression analysis of several known AhR targets indicated higher activity of the AhR pathway during craniofacial development in benthic charr morphotypes.

Conclusion

These results suggest a key role of the aryl hydrocarbon receptor (AhR) pathway in the observed craniofacial differences between distinct charr morphotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0022-6) contains supplementary material, which is available to authorized users.

Evolutionary History of Cathepsin L (L-like) Family Genes in Vertebrates

Cathepsin L family, an important cysteine protease found in lysosomes, is categorized into cathepsins B, F, H, K, L, S, and W in vertebrates. This categorization is based on their sequence alignment and traditional functional classification, but the evolutionary relationship of family members is unclear. This study determined the evolutionary relationship of cathepsin L family genes in vertebrates through phylogenetic construction. Results showed that cathepsins F, H, S and K, and L and V were chronologically diverged. Tandem-repeat duplication was found to occur in the evolutionary history of cathepsin L family. Cathepsin L in zebrafish, cathepsins S and K in xenopus, and cathepsin L in mice and rats underwent evident tandem-repeat events. Positive selection was detected in cathepsin L-like members in mice and rats, and amino acid sites under positive selection pressure were calculated. Most of these sites appeared at the connection of secondary structures, suggesting that the sites may slightly change spatial structure. Severe positive selection was also observed in cathepsin V (L2) of primates, indicating that this enzyme had some special functions. Our work provided a brief evolutionary history of cathepsin L family and differentiated cathepsins S and K from cathepsin L based on vertebrate appearance. Positive selection was the specific cause of differentiation of cathepsin L family genes, confirming that gene function variation after expansion events was related to interactions with the environment and adaptability.

The developmental transcriptome of contrasting Arctic charr (Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charrSalvelinus alpinuspopulations in Iceland.  To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, includinglysozymeandnatterin-likewhich was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

The developmental transcriptome of contrasting Arctic charr ( Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charr Salvelinus alpinus populations in Iceland.  

To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.

The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, including lysozyme and natterin-like which was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.

The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

Effects of estrogens and estrogenic disrupting compounds on fish mineralized tissues

Estrogens play well-recognized roles in reproduction across vertebrates, but also intervene in a wide range of other physiological processes, including mineral homeostasis. Classical actions are triggered when estrogens bind and activate intracellular estrogen receptors (ERs), regulating the transcription of responsive genes, but rapid non-genomic actions initiated by binding to plasma membrane receptors were recently described. A wide range of structurally diverse compounds from natural and anthropogenic sources have been shown to interact with and disrupt the normal functions of the estrogen system, and fish are particularly vulnerable to endocrine disruption, as these compounds are frequently discharged or run-off into waterways. The effect of estrogen disruptors in fish has mainly been assessed in relation to reproductive endpoints, and relatively little attention has been given to other disruptive actions. This review will overview the actions of estrogens in fish, including ER isoforms, their expression, structure and mechanisms of action. The estrogen functions will be considered in relation to mineral homeostasis and actions on mineralized tissues. The impact of estrogenic endocrine disrupting compounds on fish mineralized tissues will be reviewed, and the potential adverse outcomes of exposure to such compounds will be discussed. Current lacunae in knowledge are highlighted along with future research priorities.

Transcriptional assessment by microarray analysis and large-scale meta-analysis of the metabolic capacity of cardiac and skeletal muscle tissues to cope with reduced nutrient availability in Gilthead Sea Bream (Sparus aurata L.)

The effects of nutrient availability on the transcriptome of cardiac and skeletal muscle tissues were assessed in juvenile gilthead sea bream fed with a standard diet at two feeding levels: (1) full ration size and (2) 70 % satiation followed by a finishing phase at the maintenance ration. Microarray analysis evidenced a characteristic transcriptomic profile for each muscle tissue following changes in oxidative capacity (heart > red skeletal muscle > white skeletal muscle). The transcriptome of heart and secondly that of red skeletal muscle were highly responsive to nutritional changes, whereas that of glycolytic white skeletal muscle showed less ability to respond. The highly expressed and nutritionally regulated genes of heart were mainly related to signal transduction and transcriptional regulation. In contrast, those of white muscle were enriched in gene ontology (GO) terms related to proteolysis and protein ubiquitination. Microarray meta-analysis using the bioinformatic tool Fish and Chips ( http://fishandchips.genouest.org/index.php ) showed the close association of a representative cluster of white skeletal muscle with some of cardiac and red skeletal muscle, and many GO terms related to mitochondrial function appeared to be common links between them. A second round of cluster comparisons revealed that mitochondria-related GOs also linked differentially expressed genes of heart with those of liver from cortisol-treated gilthead sea bream. These results show that mitochondria are among the first responders to environmental and nutritional stress stimuli in gilthead sea bream, and functional phenotyping of this cellular organelle is highly promising to obtain reliable markers of growth performance and well-being in this fish species.