De Novo Transcriptome Characterization and Growth-Related Gene Expression Profiling of Diploid and Triploid Bighead Catfish (Clarias macrocephalus Günther, 1864)

Impact of gonadotropin on certain testis-related genes identified through testicular transcriptome analysis in the Asian and the African catfish

Induction of spawning in catfish is well demonstrated in females while stimulation of spermiation in males seems difficult. This has led to least understanding of gene profile changes during testicular maturation. Expression of the factors after hormonal induction that control reproductive processes has become an intense research area in fish endocrinology. In view of this, de novo assemblies of testicular transcriptomes of the Asian catfish, Clarias batrachus and the African catfish, C. gariepinus, were performed to identify genes/factors involved in testicular maturation. For testicular induction, human chorionic gonadotropin (hCG) was administered in vivo to both the catfish species during active phase of reproduction. 1,68,071 and 1,26,232 assembled transcripts were obtained from C. batrachus and C. gariepinus testicular transcriptomes and were annotated using different databases. Further, in silico analysis suggested the presence of several transcripts that were involved in steroidogenesis and male reproduction. Comparison of transcriptomes of both species revealed the presence of certain unique genes related to reproduction differentially. Expression profiles after hCG induction in testis showed higher expression of certain steroidogenesis related genes such as star, cyp11b1, cyp17a and cyp21a in both the species. Further, expression levels of crucial factors related to testis, such as dmrt1/2/3, were also found to be significantly higher after gonadotropin induction. However, amh, tgfα and dmrt4/5 expression levels were significantly low. Factors related to male reproduction i.e., tekt1, tekt2, sox9, spag1, spata2 and spata7 were found to be differentially upregulated in hCG treated testicular tissues of both the species. Histology analysis indicated that the gonadotropin induction either short or long term is a better strategy to highlight expression profile changes during testicular maturation in teleost.

Muscle Transcriptome Sequencing Revealed Thermal Stress-Responsive Regulatory Genes in Farmed Rohu, Labeo rohita (Hamilton, 1822)

Rohu, Labeo rohita, is one of the most important aquaculture species in the Indian subcontinent. Understanding the molecular-level physiological responses to thermal stress or climate change is essential. In the present work, transcriptome sequencing was carried out in the muscle tissue of the rohu in response to heat stress (35 °C) in comparison with the control (28 °C). A total of 125 Gb of sequence data was generated, and the raw-reads were filtered and trimmed, which resulted in 484 million quality reads. Reference-based assembly of reads was performed using L. rohita genome, and a total of 90.17% of reads were successfully mapped. A total of 37,462 contigs were assembled with an N50 value of 1854. The differential expression analysis revealed a total of 107 differentially expressed genes (DEGs) (15 up-, 37 down-, and 55 neutrally regulated) as compared to the control group (Log2FC > 2, P < 0.05). Gene enrichment analysis of DEGs indicates that transcripts were associated with molecular, biological, and cellular activities. The randomly selected differentially expressed transcripts were validated by RT-qPCR and found consistent expression patterns in line with the RNA-seq data. Several transcripts such as SERPINE1(HSP47), HSP70, HSP90alpha, Rano class II histocompatibility A beta, PGC-1 and ERR-induced regulator, proto-oncogene c-Fos, myozenin2, alpha-crystallin B chain-like protein, angiopoietin-like protein 8, and acetyl-CoA carboxylases have been identified in muscle tissue of rohu that are associated with stress/immunity. This study identified the key biomarker SERPINE1 (HSP47), which showed significant upregulation (~ 2- to threefold) in muscle tissue of rohu exposed to high temperature. This study can pave a path for the identification of stress-responsive biomarkers linked with thermal adaptations in the farmed carps.

Telomerase Activity in Somatic Tissues and Ovaries of Diploid and Triploid Rainbow Trout (Oncorhynchus mykiss) Females

Telomerase activity has been found in the somatic tissues of rainbow trout. The enzyme is essential for maintaining telomere length but also assures homeostasis of the fish organs, playing an important role during tissue regeneration. The unique morphological and physiological characteristics of triploid rainbow trout, when compared to diploid specimens, make them a promising model for studies concerning telomerase activity. Thus, in this study, we examined the expression of the Tert gene in various organs of subadult and adult diploid and triploid rainbow trout females. Upregulated Tert mRNA transcription was observed in all the examined somatic tissues sampled from the triploid fish when compared to diploid individuals. Contrastingly, Tert expression in the ovaries was significantly decreased in the triploid specimens. Within the diploids, the highest expression of Tert was observed in the liver and in the ovaries of the subadult individuals. In the triploids, Tert expression was increased in the somatic tissues, while the ovaries exhibited lower activity of telomerase compared to other organs and decreased compared to the ovaries in the diploids. The ovaries of triploid individuals were underdeveloped, consisting of only a few oocytes. The lack of germ cells, which are usually characterized by high Tert expression, might be responsible for the decrease in telomerase activity in the triploid ovaries. The increase in Tert expression in triploid somatic tissues suggests that they require higher telomerase activity to cope with environmental stress and maintain internal homeostasis.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Comparative Transcriptome Analysis Revealed Genes Involved in Sexual and Polyploid Growth Dimorphisms in Loach (Misgurnus anguillicaudatus)

Simple Summary

Misgurnus anguillicaudatus not only exhibits sexual size dimorphism, but also shows polyploid size dimorphism. Here, we performed comparative transcriptome integration analysis of multiple tissues of diploid and tetraploid M. anguillicaudatus of both sexes. We found that differences in energy metabolism and steroid hormone synthesis levels may be the main causes of sexual and polyploidy growth dimorphisms of M. anguillicaudatus. Fast-growing M. anguillicaudatus (tetraploids, females) have higher levels of energy metabolism and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing M. anguillicaudatus (diploids, males).

Abstract

Sexual and polyploidy size dimorphisms are widespread phenomena in fish, but the molecular mechanisms remain unclear. Loach (Misgurnus anguillicaudatus) displays both sexual and polyploid growth dimorphism phenomena, and are therefore ideal models to study these two phenomena. In this study, RNA-seq was used for the first time to explore the differentially expressed genes (DEGs) between both sexes of diploid and tetraploid loaches in four tissues (brain, gonad, liver, and muscle). Results showed that 21,003, 17, and 1 DEGs were identified in gonad, liver, and muscle tissues, respectively, between females and males in both diploids and tetraploids. Regarding the ploidy levels, 4956, 1496, 2187, and 1726 DEGs were identified in the brain, gonad, liver, and muscle tissues, respectively, between tetraploids and diploids of the same sex. When both sexual and polyploid size dimorphisms were considered simultaneously in the four tissues, only 424 DEGs were found in the gonads, indicating that these gonadal DEGs may play an important regulatory role in regulating sexual and polyploid size dimorphisms. Regardless of the sex or ploidy comparison, the significant DEGs involved in glycolysis/gluconeogenesis and oxidative phosphorylation pathways were upregulated in faster-growing individuals, while steroid hormone biosynthesis-related genes and fatty acid degradation and elongation-related genes were downregulated. This suggests that fast-growing loaches (tetraploids, females) have higher energy metabolism levels and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing loaches (diploids, males). Our findings provide an archive for future systematic research on fish sexual and polyploid dimorphisms.

Short- and long-term probiotic effects of Enterococcus hirae isolated from fermented vegetable wastes on the growth, immune responses, and disease resistance of hybrid catfish (Clarias gariepinus × Clarias macrocephalus)

This study evaluated the short- and long-term effects of dietary supplementation with Enterococcus hirae strain UPM02 on the growth performance, immunity, and disease resistance of hybrid catfish (Clarias gariepinus × Clarias macrocephalus) against Aeromonas hydrophila infection. In the long-term trial, fingerling fish were fed diets containing 0 (control), 2 × 10⁵, or 2 × 10⁷ CFU/g E. hirae UPM02 for 120 days. Administration of E. hirae UPM02 had significant effects on the specific growth rate (SGR), feed utilization efficiency, body indices (P < 0.05), and gut villus physiology of the catfish. E. hirae UPM02 application also significantly increased the complete blood cell counts, phagocytic activity, respiratory burst, lysozyme activity, and alternative complement pathway hemolytic (ACH₅₀) activity in tested catfish throughout the experimental periods (P < 0.05). Dietary E. hirae UPM02 at both concentrations significantly increased the expression levels of the alpha-2-macroglobulin (α2M), CC chemokines, CXC chemokines, lysozyme c (LYZC), myeloperoxidase (MYE), NF-kappa-B1 p105 subunit (NF-K), and bactericidal permeability-increasing protein (BPIP) genes in the head kidney, liver, and spleen (P < 0.05) at days 80, 100 and 120 after application. However, heat shock protein 70 (HSP70) gene expression was slightly downregulated in these organs. Interestingly, fish fed the diets containing 2 × 10⁵ and 2 × 10⁷ CFU/g E. hirae UPM02 exhibited a significantly lower (P < 0.05) postchallenge mortality rates (32% and 30%, respectively) after 14 days of A. hydrophila challenge than the control fish (58%). In short-term (28 days) application to juvenile catfish, the two concentrations of E. hirae did not affect all growth parameters. Nevertheless, these concentrations markedly elevated all tested immune parameters, similarly to long-term application. Immune-related gene expression was significantly upregulated at day 28 in the head kidney, at day 14 in the liver, and at day 7 in the spleen in fish treated with the two concentrations of the probiotics (P < 0.05). Mortality at 14 days after challenge with A. hydrophila in the groups receiving the two concentrations of the probiotic was significantly lower than that in the control group, at 28, 24, and 48%, respectively (P < 0.05). These results collectively suggest that dietary supplementation with E. hirae UPM02 at 2 × 10⁵ and 2 × 10⁷ CFU/g effectively influenced immune responses, enhanced disease protection, and stimulated immunity-related gene expression in hybrid catfish under both short- and long-term application. However, growth enhancement was significantly evidenced with long-term application only.

Comparative Physiological and Transcriptomic Profiling Offers Insight into the Sexual Dimorphism of Hepatic Metabolism in Size-Dimorphic Spotted Scat (Scatophagus argus)

The spotted scat (Scatophagus argus) is an economically important cultured marine fish that exhibits a typical sexual size dimorphism (SSD). SSD has captivated considerable curiosity for farmed fish production; however, up till now the exact underlying mechanism remains largely unclear. As an important digestive and metabolic organ, the liver plays key roles in the regulation of fish growth. It is necessary to elucidate its significance as a downstream component of the hypothalamic-pituitary-liver axis in the formation of SSD. In this study, the liver physiological differences between the sexes were evaluated in S. argus, and the activity of several digestive and metabolic enzymes were affected by sex. Females had higher amylase, protease, and glucose-6-phosphate dehydrogenase activities, while males exhibited markedly higher hepatic lipase and antioxidant enzymes activities. A comparative transcriptomics was then performed to characterize the responsive genes. Illumina sequencing generated 272.6 million clean reads, which were assembled into 79,115 unigenes. A total of 259 differentially expressed genes were identified and a few growth-controlling genes such as igf1 and igfbp1 exhibited female-biased expression. Further analyses showed that several GO terms and pathways associated with metabolic process, particularly lipid and energy metabolisms, were significantly enriched. The male liver showed a more active mitochondrial energy metabolism, implicating an increased energy expenditure associated with reproduction. Collectively, the female-biased growth dimorphism of S. argus may be partially attributed to sexually dimorphic metabolism in the liver. These findings would facilitate further understanding of the nature of SSD in teleost fish.

Analysis of the gene transcription patterns and DNA methylation characteristics of triploid sea cucumbers (Apostichopus japonicus)

Breeding of polyploid aquatic animals is still an important approach and research hotspot for realizing the economic benefits afforded by the improvement of aquatic animal germplasm. To better understand the molecular mechanisms of the growth of triploid sea cucumbers, we performed gene expression and genome-wide comparisons of DNA methylation using the body wall tissue of triploid sea cucumbers using RNA-seq and MethylRAD-seq technologies. We clarified the expression pattern of triploid sea cucumbers and found no dosage effect. DEGs were significantly enriched in the pathways of nucleic acid and protein synthesis, cell growth, cell division, and other pathways. Moreover, we characterized the methylation pattern changes and found 615 differentially methylated genes at CCGG sites and 447 differentially methylated genes at CCWGG sites. Integrative analysis identified 23 genes (such as Guf1, SGT, Col5a1, HAL, HPS1, etc.) that exhibited correlations between promoter methylation and expression. Altered DNA methylation and expression of various genes suggested their roles and potential functional interactions in the growth of triploid sea cucumbers. Our data provide new insights into the epigenetic and transcriptomic alterations of the body wall tissue of triploid sea cucumbers and preliminarily elucidate the molecular mechanism of their growth, which is of great significance for the breeding of fine varieties of sea cucumbers.

Comparative Genomic and Transcriptomic Analyses Revealed Twenty-Six Candidate Genes Involved in the Air-Breathing Development and Function of the Bighead Catfish Clarias macrocephalus

The bighead catfish (Clarias macrocephalus) and channel catfish (Ictalurus punctatus) are freshwater species in the Siluriformes order. C. macrocephalus has both gills and modified gill structures serving as an air-breathing organ (ABO), while I. punctatus does not possess such an organ, and cannot breathe in air, providing an excellent model for studying the molecular basis of ABO development in teleost fish. To investigate the critical time window for the development of air-breathing function, seven development stages were selected based on hypoxia challenge results, and RNA-seq was performed upon C. macrocephalus to compare with the non-air-breathing I. punctatus. Five-hundred million reads were generated and 25,239 expressed genes were annotated in C. macrocephalus. Among those, 8675 genes were differentially expressed across developmental stages. Comparative genomic analysis identified 1458 C. macrocephalus specific genes, which were absent in I. punctatus. Gene network and protein-protein interaction analyses identified 26 key hub genes involved in the air-breathing function. Three top candidate genes, mb, ngb, hbae, are mainly associated with oxygen carrying, oxygen binding, and heme binding activities. Our study provides a rich data set for exploring the genomic basis of air-breathing function in C. macrocephalus and offers insights into the adaption to hypoxic environments.

Comparative transcriptomic analysis for identification of candidate sex-related genes and pathways in Crimson seabream (Parargyrops edita)

Teleost fishes display the largest array of sex-determining systems among animals, resulting in various reproductive strategies. Research on sex-related genes in teleosts will broaden our understanding of the process, and provide important insight into the plasticity of the sex determination process in vertebrates in general. Crimson seabream (Parargyrops edita Tanaka, 1916) is one of the most valuable and abundant fish resources throughout Asia. However, little genomic information on P. edita is available. In the present study, the transcriptomes of male and female P. edita were sequenced with RNA-seq technology. A total of 388,683,472 reads were generated from the libraries. After filtering and assembling, a total of 79,775 non redundant unigenes were obtained with an N50 of 2,921 bp. The unigenes were annotated with multiple public databases, including NT (53,556, 67.13%), NR (54,092, 67.81%), Swiss-Prot (45,265, 56.74%), KOG (41,274, 51.74%), KEGG (46,302, 58.04%), and GO (11,056, 13.86%) databases. Comparison of the unigenes of different sexes of P. edita revealed that 11,676 unigenes (9,335 in females, 2,341 in males) were differentially expressed between males and females. Of these, 5,463 were specifically expressed in females, and 1,134 were specifically expressed in males. In addition, the expression levels of ten unigenes were confirmed to validate the transcriptomic data by qRT-PCR. Moreover, 34,473 simple sequence repeats (SSRs) were identified in SSR-containing sequences, and 50 loci were randomly selected for primer development. Of these, 36 loci were successfully amplified, and 19 loci were polymorphic. Finally, our comparative analysis identified many sex-related genes (zps, amh, gsdf, sox4, cyp19a, etc.) and pathways (MAPK signaling pathway, p53 signaling pathway, etc.) of P. edita. This informative transcriptomic analysis provides valuable data to increase genomic resources of P. edita. The results will be useful for clarifying the molecular mechanism of sex determination and for future functional analyses of sex-associated genes.

Transcriptome sequencing and histology reveal dosage compensation in the liver of triploid pre-smolt Atlantic salmon

Triploid Atlantic salmon (Salmo salar L.) is seen as one of the best solutions to solve key issues in the salmon farming industry, such as the impact of escapees on wild stocks and pre-harvest sexual maturation. However, the effects of triploidy on salmon smoltification are poorly understood at the molecular level, even though smoltification is a very sensitive period that has a major influence on survival rate and performance of farmed salmon. In this study, we have compared the liver transcriptomes of diploid and triploid Atlantic salmon at three ontogeny stages: fry, parr and smolt. In diploid fish, a total of 2,655 genes were differentially expressed between fry and parr, whereas 506 genes had significantly different transcript levels between parr and smolts. In triploids, 1,507 and 974 genes were differentially expressed between fry and parr, and between parr and smolts, respectively. Most of these genes were down-regulated and 34 genes were differentially expressed between ploidies at the same stage. In both ploidy groups, the top differentially expressed genes with ontogeny stage belonged to common functional categories that can be related to smoltification. Nucleotide and energy metabolism were significantly down-regulated in fry when compared to parr, while immune system processes were significantly down-regulated in parr when compared to smolts. The close resemblance of enriched biological processes and pathways between ploidy groups suggests that triploidy is regulated by genome dosage compensation in Atlantic salmon. Histological analysis revealed that areas of vacuolization (steatosis) were present only in fry and parr stages, in contrast to a compact cellular histology with glycogen granules after smoltification. There was no significant difference in vacuolization between ploidy groups at the fry stage but the liver of diploid parr had a 33.5% higher vacuolization area compared to their triploid counterparts. Taken together, our data provide novel insights into the changes that occur at the molecular and histological level in the liver of both diploid and triploid Atlantic salmon prior to and during smoltification.

A transcriptomic analysis of diploid and triploid Atlantic salmon lenses with and without cataracts

To avoid negative environmental impacts of escapees and potential inter-breeding with wild populations, the Atlantic salmon farming industry has and continues to extensively test triploid fish that are sterile. However, they often show differences in performance, physiology, behavior and morphology compared to diploid fish, with increased prevalence of vertebral deformities and ocular cataracts as two of the most severe disorders. Here, we investigated the mechanisms behind the higher prevalence of cataracts in triploid salmon, by comparing the transcriptional patterns in lenses of diploid and triploid Atlantic salmon, with and without cataracts. We assembled and characterized the Atlantic salmon lens transcriptome and used RNA-seq to search for the molecular basis for cataract development in triploid fish. Transcriptional screening showed only modest differences in lens mRNA levels in diploid and triploid fish, with few uniquely expressed genes. In total, there were 165 differentially expressed genes (DEGs) between the cataractous diploid and triploid lens. Of these, most were expressed at lower levels in triploid fish. Differential expression was observed for genes encoding proteins with known function in the retina (phototransduction) and proteins associated with repair and compensation mechanisms. The results suggest a higher susceptibility to oxidative stress in triploid lenses, and that mechanisms connected to the ability to handle damaged proteins are differentially affected in cataractous lenses from diploid and triploid salmon.

Transcriptome analysis of Clarias magur brain and gonads suggests neuro-endocrine inhibition of milt release from captive GnRH-induced males

Transcriptome analysis of Clarias magur brain and gonads at preparatory, mature, 6 and 16 h post-GnRH injection (hpi) stages yielded 9.5 GB data with 39,738 contigs. Sequences of 45 reproductive genes were identified for the first time in C. magur along with unique and differentially expressed genes. The expression of 20 genes was validated by qRT-PCR. Upregulation of Cyp11A1, Cyp17A1 and FTZF1 genes in the 16hpi testis accompanied by the 17β-HSD3 expression indicates testosterone (T) synthesis in response to LH surge, while reduced expression of CYP11B1 suggests a high T: 11-KT ratio. It is evident by the gene expression analysis that the inhibitory neurotransmitter GABA, altered T: 11-KT, increased testicular bile acids, and oxytocin-like neuropeptide in the male brain, appear to be involved in arresting the pulsatile motion of testicular smooth muscles. The work generates important leads for an effective induced breeding strategy for silurid catfish.

Development of a Monoclonal Antibody Specific to the IgM Heavy Chain of Bighead Catfish (Clarias macrocephalus): A Biomolecular Tool for the Detection and Quantification of IgM Molecules and IgM+ Cells in Clarias Catfish

Catfish is a commonly-cultivated freshwater fish in Thailand and many Southeast Asian countries. The molecular data obtained for the IgM heavy chain (IgMH) of catfish have been useful for distinguishing monoclonal antibodies (mAbs). A mAb specific to Cμ₁ of the IgMH of catfish (IgMHCμ₁ mAb) was developed in a rabbit model using sequence information from bighead catfish (Clarias macrocephalus). The IgMHCμ₁ mAb strongly recognized the IgM heavy chain of the tested catfish, namely, bighead catfish, African catfish (Clarias gariepinus) and their hybrid (C. macrocephalus × C. gariepinus), in immunological Western blot analysis and competitive ELISAs. Additionally, the IgMHCμ₁ mAb successfully recognized IgM⁺ cells by detecting IgM molecules in both secreted and membrane-bound forms in peripheral blood leukocytes (PBLs). The IgMHCμ₁ mAb was further used to quantify the percentage of IgM⁺ cells among PBLs through flow cytophotometry. The IgM⁺ cell percentages of healthy bighead catfish, African catfish and their hybrid were 38.0–39.9%, 45.6–53.2%, and 58.7–60.0%, respectively. Furthermore, the IgMHCμ₁ mAb showed no cross-reactivity with the IgM of zebrafish. These findings suggest that this mAb can be used as an immunological tool for monitoring the health, immune status, and immune development of cultivated Clarias catfish.

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.

Probiotic Effects of a Novel Strain, Acinetobacter KU011TH, on the Growth Performance, Immune Responses, and Resistance against Aeromonas hydrophila of Bighead Catfish (Clarias macrocephalus Günther, 1864)

In the present study, the novel probiotic strain Acinetobacter KU011TH with an evident lack of pathogenicity in catfish was experimented. Three practical administration routes, namely, feed additive (FD), water-soluble additive (SOL), and a combination route (FD+SOL), were applied in two sizes of catfish. After 120 days of FD+SOL administration, catfish fingerlings (15 g) exhibited a significant improvement in all tested growth performance parameters. For 15- and 30-day applications at the juvenile stage (150 g), phagocytic activity, phagocytic index, lysozyme activity, respiratory burst activity, alternative complement pathway, and bactericidal activity were significantly increased. Furthermore, probiotic-administered bighead catfish exhibited an upregulated expression of several immune-related genes in tested organs. Significant colonization by Acinetobacter KU011TH in rearing water and on skin and gills was observed among experimental groups. Histological analysis clearly indicated enhanced physical characteristics of skin mucosal immunity in the treated groups. No histopathological changes in the gills, skin, intestine or liver were observed among the fish groups. Interestingly, after challenge with Aeromonas hydrophila, the survival rates of the treated groups were significantly higher than those of the controls. In conclusion, the novel probiont Acinetobacter KU011TH provides a potent strategy for improvement in growth and disease resistance, which is an important steppingstone for sustaining catfish aquaculture.

An SNP-Based Genetic Map and QTL Mapping for Growth Traits in the Red-Spotted Grouper (Epinephelus akaara)

The red-spotted grouper (Epinephelus akaara) is one of the most commercially important aquatic species in China. However, its seedstock has low larval survival rates, and its stability is confronted with the danger of overexploitation. In this study, a high-density genetic map was constructed using 3435 single nucleotide polymorphisms (SNPs) from 142 first generation (F₁) full-sib offspring and two parents of a red-spotted grouper population. The total genetic length of the map was 2300.12 cM with an average intermarker distance of 0.67 cM. Seventeen genome-wide significant quantitative trait loci (QTLs) for growth-related traits were detected on 24 linkage groups, including 5 QTLs for full length, 7 QTLs for body length, and 5 QTLs for body weight. The contribution values of explained phenotypic variance ranged from 10.7% to 12.9%. Moreover, 13 potential candidate genes for growth-related traits were identified. Collectively, these findings will be useful for conducting marker-assisted selection of the red-spotted grouper in future studies.

Differential Transcriptomic and Metabolomic Responses in the Liver of Nile Tilapia (Oreochromis niloticus) Exposed to Acute Ammonia

Ammonia is toxic to aquatic animal. Currently, only limited works were reported on the responses of aquatic animals after ammonia exposure using "omics" technologies. Tilapia suffers from the stress of ammonia-nitrogen during intensive recirculating aquaculture. Optimizing ammonia stress tolerance has become an important issue in tilapia breeding. The molecular and biochemical mechanisms of ammonia-nitrogen toxicity have not been understood comprehensively in tilapia yet. In this study, using RNA-seq and gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF-MS) techniques, we investigated differential expressed genes (DEGs) and metabolomes in the liver at 6 h post-challenges (6 hpc) and 24 h post-challenges (24 hpc) under high concentration of ammonia-nitrogen treatment. We detected 2258 DEGs at 6 hpc and 315 DEGs at 24 hpc. Functional enrichment analysis indicated that DEGs were significantly associated with cholesterol biosynthesis, steroid and lipid metabolism, energy conservation, and mitochondrial tissue organization. Metabolomic analysis detected 31 and 36 metabolites showing significant responses to ammonia-nitrogen stress at 6 and 24 hpc, respectively. D-(Glycerol 1-phosphate), fumaric acid, and L-malic acid were found significantly down-regulated at both 6 and 24 hpc. The integrative analysis of transcriptomics and metabolomics suggested considerable alterations and precise control of gene expression at both physiological and molecular levels in response to the stress of ammonia-nitrogen in tilapia.

Hepatic transcriptome profile of sheep (Ovis aries) in response to overgrazing: novel genes and pathways revealed

BACKGROUND

Overgrazing is a major factor that causes steppe degradation in Inner Mongolian, resulting in extensive ecosystem damage. Scarcity of grass means sheep are smaller and therefore mutton and cashmere production is greatly reduced, which has resulted in massive annual economic losses. Liver is the primary metabolic organ in mammals. It is also the key source of energy supply and detoxification of metabolites in animals, has a close relationship with animal growth. However, investigations on the responses of sheep induced by consequence of overgrazing, particularly those relating to liver-related molecular mechanisms and related metabolic pathways, remain elusive.

RESULTS

The body weight daily gain of sheep, immune organ indices (liver and spleen), and serum parameters related to immune response, protein synthesis and energy supply (IgG, albumin, glucose and non-esterified fatty acid) were significantly lower in the overgrazing group. Other serum parameters including alanine aminotransferase, aspartate aminotransferase, alkaline phosphatase, total bilirubin, blood urea nitrogen and interleukin-6 were significantly higher in the overgrazing group. For the RNA-Seq results, we identified approximately 50 differentially expressed genes, of which half of were up-regulated and the other half were down-regulated (overgrazing group versus light grazing group). Bioinformatics analysis identified two enriched KEGG pathways including peroxisome proliferator-activated receptor (PPAR) signaling pathway (related to lipolysis) and ECM-receptor interaction (related to liver injury and apoptosis). Additionally, several of the down-regulated genes were related to detoxification and immune response.

CONCLUSIONS

Overall, based on the high-throughput RNA sequencing profile integrated with the results of serum biochemical analyses, consequences of lower forage availability and quality under overgrazing condition induced altered expression levels of genes participating in energy metabolism (particularly lipid metabolism) and detoxification and immune responses, causing lipolysis and impaired health status, which might be key reasons for the reduced growth performance of sheep. This investigation provides a novel foundation for the development of sheep hepatic gene interactive networks that are a response to the degraded forage availability under overgrazing condition.

Effect of triploidy on liver gene expression in coho salmon (Oncorhynchus kisutch) under different metabolic states

BACKGROUND

Triploid coho salmon are excellent models for studying gene dosage and the effects of increased cell volume on gene expression. Triploids have an additional haploid genome in each cell and have fewer but larger cells than diploid coho salmon to accommodate the increased genome size. Studying gene expression in triploid coho salmon provides insight into how gene expression may have been affected after the salmonid-specific genome duplication which occurred some 90 MYA. Triploid coho salmon are sterile and consequently can live longer and grow larger than diploid congeners in many semelparous species (spawning only once) because they never reach maturity and post-spawning mortality is averted. Triploid fishes are also of interest to the commercial sector (larger fish are more valuable) and to fisheries management since sterile fish can potentially minimize negative impacts of escaped fish in the wild.

RESULTS

The vast majority of genes in liver tissue had similar expression levels between diploid and triploid coho salmon, indicating that the same amount of mRNA transcripts were being produced per gene copy (positive gene dosage effects) within a larger volume cell. Several genes related to nutrition and compensatory growth were differentially expressed between diploid and triploid salmon, indicating that some loci are sensitive to cell size and/or DNA content per cell. To examine how robust expression between ploidies is under different conditions, a genetic/metabolic modifier in the form of different doses of a growth hormone transgene was used to assess gene expression under conditions that the genome has not naturally experienced or adapted to. While many (up to 1400) genes were differentially expressed between non-transgenic and transgenic fish, relatively few genes were differentially expressed between diploids and triploids with similar doses of the transgene. These observations indicate that the small effect of ploidy on gene expression is robust to large changes in physiological state.

CONCLUSIONS

These findings are of interest from a gene regulatory perspective, but also valuable for understanding phenotypic effects in triploids, transgenics, and triploid transgenics that could affect their utility in culture conditions and their fitness and potential consequences of release into nature.

Comparative transcriptome analysis of diploid and triploid hybrid groupers (Epinephelus coioides♀ × E. lanceolatus♂) reveals the mechanism of abnormal gonadal development in triploid hybrids

In our previous studies, diploid and triploid hybrids have been detected from the hybridization of Epinephelus coioides♀ × E. lanceolatus♂. The triploid groupers have been found to be delayed in gonadal development, but the mechanism remains poorly understood. In this study, we examined the gonadal development, assayed the serum steroid hormone levels, and compared the BP (brain and pituitary) and G (gonad) transcriptomes of 18-month-old diploid and triploid hybrids. The results showed that levels of serum estradiol-17β and testosterone were significantly higher in triploid groupers. The RNA-seq data revealed that 1518 and 14,963 differentially expressed genes were identified in the BP and G transcriptome, respectively. Further analysis revealed that the expression levels of genes involved in the sexual differentiation pathway and sex steroid synthesis pathway are significantly higher in triploid hybrids. Our findings provided a comprehensive insight into a better understanding of the regulatory mechanisms of sterility in triploid hybrid fish.

Gonadal Transcriptome Analysis of Pacific Abalone Haliotis discus discus: Identification of Genes Involved in Germ Cell Development

Little is known about the molecular mechanisms governing gonadal developmental processes in abalones. Here, we conducted transcriptome analysis of Pacific abalone Haliotis discus discus for gene discovery in the brain, ovary, testis, and unfertilized eggs. Among the annotated unigenes, 48.6% of unigenes were identified by Venn diagram analysis as having universal or tissue-specific expression. Twenty-three genes with gonad-biased gene ontology (GO) terms were first obtained. Secondly, 36 genes were found by screening known gene names related to germ cell development. Finally, 17 genes were obtained by querying the annotated unigene database for zygotically expressed gonadal genes (ovary and testis) and maternally expressed gonadal genes (ovary, testis, and unfertilized eggs) using keywords related to reproduction. To further verify tissue distribution pattern and subcellular localization of these genes, RT-PCR and in situ hybridization were performed using a unigene encoding a germ cell marker, vasa, as control. The results showed that vasa was expressed mainly in the early developmental stages of germ cells in both sexes. One of the candidate genes, vitelline envelope zona pellucida domain protein 12 (ZP12), was expressed in the primordial germ cells of immature gonad and early developmental stages of germ cells of the adult female. The results obtained from the present study suggest that vasa and ZP12 are involved in germ cell development of Pacific abalone and that ZP12 is an especially useful germ cell-specific marker in immature adults. The current gonadal transcriptome profile is an extensive resource for future reproductive molecular biology studies of this species.

Comparative Transcriptome Analysis Reveals Growth-Related Genes in Juvenile Chinese Sea Cucumber, Russian Sea Cucumber, and Their Hybrids

Heterosis is important for sea cucumber breeding, but its molecular mechanism remains largely unexplored. In this study, parental lines of Apostichopus japonicus from Russia (R) and China (C) were used to construct hybrids (CR and RC) by reciprocal crossing. We examined the transcriptional profiles of the hybrids (CR and RC) and the purebreds (CC and RR) at different developmental times. A total of 60.27 Gb of clean data was obtained, and 176,649 unigenes were identified, of which 50,312 unigenes were annotated. A total of 414,536 SNPs were identified. A total of 7011 differentially expressed genes (DEGs) were obtained between the purebreds and hybrids at 45 days after fertilization (DAF), and a total of 8218 DEGs were obtained between the purebreds and hybrids at 75 DAF. In addition, a total of 7652 DEGs were obtained between 45 DAF and 75 DAF. The significant DEGs were mainly involved in the MAPK and FOXO signaling pathways, especially in the Ras-Raf-MEK1/2-ERK module, which may be a key regulator of development and growth in juvenile A. japonicus. In addition, we also identified key growth-related genes, such as fgfs, igfs, megfs and hgfs, which were upregulated in the hybrids (RC and CR); these genes may play important roles in heterosis in A. japonicus. Our study provides fundamental information on the molecular mechanisms underlying heterosis in sea cucumber and might suggest strategies for the selection of rapidly growing strains of sea cucumber in aquaculture.