Identification and expression of Smads associated with TGF-β/activin/nodal signaling pathways in the rainbow trout (Oncorhynchus mykiss)

Transcriptomic assessment of dietary fishmeal partial replacement by soybean meal and prebiotics inclusion in the liver of juvenile Pacific yellowtail (Seriola lalandi)

BACKGROUND

Seriola lalandi is an important species for aquaculture, due to its rapid growth, adaptation to captivity and formulated diets, and high commercial value. Due to the rise in fishmeal (FM) price, efforts have been and still are made to replace it partially or entirely with vegetable meals in diets for carnivorous fish. The use of prebiotics when feeding vegetable meals has improved fish health.

METHODS

Four experimental diets were assessed in juveniles, the control diet consisted of FM as the main protein source, the second diet included 2% of GroBiotic®-A (FM-P), in the third diet FM was partially replaced (25%) by soybean meal (SM25), and the fourth consisted of SM25 with 2% of GroBiotic®-A (SM25-P). Growth was evaluated and RNA-seq of the liver tissue was performed, including differential expression analysis and functional annotation to identify genes affected by the diets.

RESULTS

Growth was not affected by this level of FM replacement, but it was improved by prebiotics. Annotation was achieved for 59,027 transcripts. Gene expression was affected by the factors: 225 transcripts due to FM replacement, 242 due to prebiotics inclusion, and 62 due to the interaction of factors. The SM25-P diet showed the least amount of differentially expressed genes against the control diet.

CONCLUSION

The replacement of FM (25%) by soybean meal combined with prebiotics (2%) represents a good cost-benefit balance for S. lalandi juveniles since the fish growth increased and important metabolic and immune system genes in the liver were upregulated with this diet.

Detection of selection signatures in the genome of a farmed population of anadromous rainbow trout (Oncorhynchus mykiss)

Domestication processes and artificial selection are likely to leave signatures that can be detected at a molecular level in farmed rainbow trout (Oncorhynchus mykiss). These signatures of selection are genomic regions that contain functional genetic variants conferring a higher fitness to their bearers. We genotyped 749 rainbow trout from a commercial population using a rainbow trout Axiom 57 K SNP array panel and identified putative genomic regions under selection using the pcadapt, Composite Likelihood Ratio (CLR) and Integrated Haplotype Score (iHS) methods. After applying quality-control pipelines and statistical analyses, we detected 12, 96 and 16 SNPs putatively under selection, associated with 96, 781 and 115 candidate genes, respectively. Several of these candidate genes were associated with growth, early development, reproduction, behavior and immune system traits. In addition, some of the SNPs were found in interesting regions located in autosomal inversions on Omy05 and Omy20. These findings could represent a genome-wide map of selection signatures in farmed rainbow trout and could be important in explaining domestication and selection for genetic traits of commercial interest.

Identification, molecular evolution, and expression analysis of the transcription factor Smad gene family in lamprey

Transcription factor small mothers against decapentaplegic (Smad) family SMAD proteins are the essential intracellular signal mediators and transcription factors for transforming growth factor β (TGF-β) signal transduction pathway, which usually exert pleiotropic actions on cell physiology, including immune response, cell migration and differentiation. In this study, the Smad family was identified in the most primitive vertebrates through the investigation of the transcriptome data of lampreys. The topology of phylogenetic tree showed that the four Smads (Smad1, Smad3, Smad4 and Smad6) in lampreys were subdivided into four different groups. Meanwhile, homology analysis indicated that most Smads were conserved with typical Mad Homology (MH) 1 and MH2 domains. In addition, Lethenteron reissneri Smads (Lr-Smads) adopted general Smads folding structure and had high tertiary structural similarity with human Smads (H-Smads). Genomic synteny analysis revealed that the large-scale duplication blocks were not found in lamprey genome and neighbor genes of lamprey Smads presented dramatic differences compared with jawed vertebrates. Importantly, quantitative real-time PCR analysis demonstrated that Smads were widely expressed in lamprey, and the expression level of Lr-Smads mRNA was up-regulated with different pathogenic stimulations. Moreover, depending on the weighted gene co-expression network analysis (WGCNA), four Lr-Smads were identified as two meaningful modules (green and gray). The functional analysis of these two modules showed that they might have a correlation with ployI:C. And these genes presented strong positive correlation during the immune response from the results of Pearson's correlation analysis. In conclusion, our results would not only enrich the information of Smad family in jawless vertebrates, but also lay the foundation for immunity in further study.

Expression and regulation of Smad2 by gonadotropins in the protogynous hermaphroditic ricefield eel (Monopterus albus)

Smad2, a receptor-activated Smad, plays a critical role in regulating gametogenesis. In this study, a smad2 homologue was identified and sequenced from ricefield eel ovary cDNA, and its mRNA and protein expression levels were analysed during oocyte development. The cDNA sequence of ricefield eel smad2 consisted of 1863 bp encoding a 467-amino acid protein that had high sequence homology with Smad proteins in other teleosts, especially in Poeciliopsis prolifica. The results of real-time quantitative PCR (RT-qPCR) analysis revealed that smad2 is expressed in the ovary during gonad development, increased continuously until the early vitellogenic stage in the ovaries, and then decreased with ovary maturation. Smad2 protein immunoreactivity was localized in the cytoplasm of follicular cells, oogonia, and primary growth stage oocytes. In vitro experiments revealed that follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG) promoted smad2 expression in ovary tissue in a time- and dose-dependent manner, respectively. In summary, Smad2 plays a potentially vital role in ricefield eel ovary development.

The functional differentiation of four smad4 paralogs in TGF-β signaling pathway of Japanese flounder (Paralichthys olivaceus)

As a classical signaling pathway, transforming growth factor β (TGF-β) has been studied in various animals for more than decade years. However, the members of TGF-β were markedly expanded in teleost specific third and fourth rounds of whole genome duplication (WGD). Here, four smad4s named Posmad4a, Posmad4b, Posmad4c and Posmad4d were identified in Japanese flounder. Our study showed that four flounder smad4s had distinct properties in terms of their protein structure, expression pattern, protein interaction and subcellular localization. PoSMAD4a/b were mainly located in the cytoplasm, and could co-localize in the nucleus with PoSMAD3a after TGF-β activator stimulation. PoSMAD4c was mainly located in nucleus, whereas PoSMAD4d distributed in the whole cell. Both PoSMAD4c and PoSMAD4d could co-localize in the nucleus with PoSMAD3b after TGF-β activator stimulation. Furthermore, Posmad4c responded most strongly to TGF-β signal stimulation. Dual-luciferase reporter assay also showed that Posmad4c could specifically up-regulate the TGF-β signal luciferase reporter gene, Posmad4b could enhance Wnt signal luciferase reporter gene, while both Posmad4b and Posmad4d could markedly up-regulate Notch signal reporter gene. All results indicated that Posmad4a/b/c/d had significantly functional differences among TGF-β, Notch and Wnt signaling pathways. Our study provided important understanding to the biology of smad4s and its pathway crosstalk in teleost.

Effects of sex steroids on expression of genes regulating growth-related mechanisms in rainbow trout (Oncorhynchus mykiss)

Effects of a single injection of 17β-estradiol (E2), testosterone (T), or 5β-dihydrotestosterone (DHT) on expression of genes central to the growth hormone (GH)/insulin-like growth factor (IGF) axis, muscle-regulatory factors, transforming growth factor-beta (TGFβ) superfamily signaling cascade, and estrogen receptors were determined in rainbow trout (Oncorhynchus mykiss) liver and white muscle tissue. In liver in addition to regulating GH sensitivity and IGF production, sex steroids also affected expression of IGF binding proteins, as E2, T, and DHT increased expression of igfbp2b and E2 also increased expression of igfbp2 and igfbp4. Regulation of this system also occurred in white muscle in which E2 increased expression of igf1, igf2, and igfbp5b1, suggesting anabolic capacity may be maintained in white muscle in the presence of E2. In contrast, DHT decreased expression of igfbp5b1. DHT and T decreased expression of myogenin, while other muscle regulatory factors were either not affected or responded similarly for all steroid treatments. Genes within the TGFβ superfamily signaling cascade responded to steroid treatment in both liver and muscle, suggesting a regulatory role for sex steroids in the ability to transmit signals initiated by TGFβ superfamily ligands, with a greater number of genes responding in liver than in muscle. Estrogen receptors were also regulated by sex steroids, with era1 expression increasing for all treatments in muscle, but only E2- and T-treatment in liver. E2 reduced expression of erb2 in liver. Collectively, these data identify how physiological mechanisms are regulated by sex steroids in a manner that promotes the disparate effects of androgens and estrogens on growth in salmonids.

Regulation of skeletal muscle growth in fish by the growth hormone--insulin-like growth factor system

The growth hormone (GH)-insulin-like growth factor (IGF) system is the key promoter of growth in vertebrates; however, how this system modulates muscle mass in fish is just recently becoming elucidated. In fish, the GH induces muscle growth by modulating the expression of several genes belonging to the myostatin (MSTN), atrophy, GH, and IGF systems as well as myogenic regulatory factors (MRFs). The GH controls the expression of igf1 via Janus kinase 2 (JAK2)/signal transducers and activators of the transcription 5 (STAT5) signaling pathway, but it seems that it is not the major regulator. These mild effects of the GH on igf1 expression in fish muscle seem to be related with the presence of higher contents of truncated GH receptor1 (tGHR1) than full length GHR (flGHR1). IGFs in fish stimulate myogenic cell proliferation, differentiation, and protein synthesis through the MAPK/ERK and PI3K/AKT/TOR signaling pathways, concomitant with abolishing protein degradation and atrophy via the PI3K/AKT/FOXO signaling pathway. Besides these signaling pathways control the expression of several genes belonging to the atrophy and IGF systems. Particularly, IGFs and amino acid control the expression of igf1, thus, suggesting other of alternative signaling pathways regulating the transcription of this growth factor. The possible role of IGF binding proteins (IGFBPs) and the contribution of muscle-derived versus hepatic-produced IGF1 on fish muscle growth is also addressed. Thus, a comprehensive overview on the GH-IGF system regulating fish skeletal muscle growth is presented, as well as perspectives for future research in this field.

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

An articulated endoskeleton that is calcified is a unifying innovation of the vertebrates, however the molecular basis of the structural divergence between terrestrial and aquatic vertebrates, such as teleost fish, has not been determined. In the present study long-read next generation sequencing (NGS, Roche 454 platform) was used to characterize acellular perichondral bone (vertebrae) and chondroid bone (gill arch) in the gilthead sea bream (Sparus auratus). A total of 15.97 and 14.53Mb were produced, respectively from vertebrae and gill arch cDNA libraries and yielded 32,374 and 28,371 contigs (consensus sequences) respectively. 10,455 contigs from vertebrae and 10,625 contigs from gill arches were annotated with gene ontology terms. Comparative analysis of the global transcriptome revealed 4249 unique transcripts in vertebrae, 4201 unique transcripts in the gill arches and 3700 common transcripts. Several core gene networks were conserved between the gilthead sea bream and mammalian skeleton. Transcripts for putative endocrine factors were identified in acellular gilthead sea bream bone suggesting that in common with mammalian bone it can act as an endocrine tissue. The acellular bone of the vertebra, in contrast to current opinion based on histological analysis, was responsive to a short fast and significant (p<0.05) down-regulation of several transcripts identified by NGS, osteonectin, osteocalcin, cathepsin K and IGFI occurred. In gill arches fasting caused a significant (p<0.05) down-regulation of osteocalcin and up-regulation of MMP9.

Genomic patterns of introgression in rainbow and westslope cutthroat trout illuminated by overlapping paired-end RAD sequencing

Rapid and inexpensive methods for genomewide single nucleotide polymorphism (SNP) discovery and genotyping are urgently needed for population management and conservation. In hybridized populations, genomic techniques that can identify and genotype thousands of species-diagnostic markers would allow precise estimates of population- and individual-level admixture as well as identification of 'super invasive' alleles, which show elevated rates of introgression above the genomewide background (likely due to natural selection). Techniques like restriction-site-associated DNA (RAD) sequencing can discover and genotype large numbers of SNPs, but they have been limited by the length of continuous sequence data they produce with Illumina short-read sequencing. We present a novel approach, overlapping paired-end RAD sequencing, to generate RAD contigs of >300-400 bp. These contigs provide sufficient flanking sequence for design of high-throughput SNP genotyping arrays and strict filtering to identify duplicate paralogous loci. We applied this approach in five populations of native westslope cutthroat trout that previously showed varying (low) levels of admixture from introduced rainbow trout (RBT). We produced 77 141 RAD contigs and used these data to filter and genotype 3180 previously identified species-diagnostic SNP loci. Our population-level and individual-level estimates of admixture were generally consistent with previous microsatellite-based estimates from the same individuals. However, we observed slightly lower admixture estimates from genomewide markers, which might result from natural selection against certain genome regions, different genomic locations for microsatellites vs. RAD-derived SNPs and/or sampling error from the small number of microsatellite loci (n = 7). We also identified candidate adaptive super invasive alleles from RBT that had excessively high admixture proportions in hybridized cutthroat trout populations.