Transcriptome Analysis for Identification of Genes Related to Gonad Differentiation, Growth, Immune Response and Marker Discovery in The Turbot (Scophthalmus maximus)

Whole-Genome Sequencing Analyses Reveal the Evolution Mechanisms of Typical Biological Features of Decapterus maruadsi

Decapterus maruadsi is a typical representative of small pelagic fish characterized by fast growth rate, small body size, and high fecundity. It is a high-quality marine commercial fish with high nutritional value. However, the underlying genetics and genomics research focused on D. maruadsi is not comprehensive. Herein, a high-quality chromosome-level genome of a male D. maruadsi was assembled. The assembled genome length was 716.13 Mb with contig N50 of 19.70 Mb. Notably, we successfully anchored 95.73% contig sequences into 23 chromosomes with a total length of 685.54 Mb and a scaffold N50 of 30.77 Mb. A total of 22,716 protein-coding genes, 274.90 Mb repeat sequences, and 10,060 ncRNAs were predicted, among which 22,037 (97%) genes were successfully functionally annotated. The comparative genome analysis identified 459 unique, 73 expanded, and 52 contracted gene families. Moreover, 2804 genes were identified as candidates for positive selection, of which some that were related to the growth and development of bone, muscle, cardioid, and ovaries, such as some members of the TGF-β superfamily, were likely involved in the evolution of typical biological features in D. maruadsi. The study provides an accurate and complete chromosome-level reference genome for further genetic conservation, genomic-assisted breeding, and adaptive evolution research for D. maruadsi.

Sequencing-based network analysis provides a core set of genes for understanding hemolymph immune response mechanisms against Poly I:C stimulation in Amphioctopus fangsiao

Aquatic viruses can spread rapidly and widely in seawater for their high infective ability. Polyinosinic-polycytidylic acid (Poly I:C), a viral dsRNA analog, is an immunostimulant that has been proved to activate various immune responses of immune cells in invertebrate. Hemolymph is a critical site that host immune response in invertebrates, and its transcriptome information obtained from Amphioctopus fangsiao stimulated by Poly I:C is crucial for understanding the antiviral molecular mechanisms of this species. In this study, we analyzed gene expression data in A. fangsiao hemolymph tissue within 24 h under Poly I:C stimulation and found 1082 and 299 differentially expressed genes (DEGs) at 6 and 24 h, respectively. Union set (1,369) DEGs were selected for subsequent analyses. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses were carried out for identifying DEGs related to immunity. Several significant immune-related terms and pathways, such as toll-like receptor signaling pathways term, inflammatory response term, TNF signaling pathway, and chemokine signaling pathway were identified. A protein-protein interaction (PPI) network was constructed for examining the relationships among immune-related genes. Finally, 12 hub genes, including EGFR, ACTG1, MAP2K1, and other nine hub genes, were identified based on the KEGG enrichment analysis and PPI network. The quantitative RT-PCR (qRT-PCR) was used to verify the expression profile of 12 hub genes. This research provides a reference for solving the problem of high mortality of A. fangsiao and other mollusks and provides a reference for the future production of some disease-resistant A. fangsiao.

Comparative transcriptome sequencing analysis of female and male Decapterus macrosoma

Sexual growth dimorphism is a common phenomenon in teleost fish and has led to many reproductive strategies. Growth- and sex-related gene research in teleost fish would broaden our understanding of the process. In this study, transcriptome sequencing of shortfin scad Decapterus macrosoma was performed for the first time, and a high-quality reference transcriptome was constructed. After identification and assembly, a total of 58,475 nonredundant unigenes were obtained with an N50 length of 2,266 bp, and 28,174 unigenes were successfully annotated with multiple public databases. BUSCO analysis determined a level of 92.9% completeness for the assembled transcriptome. Gene expression analysis revealed 2,345 differentially expressed genes (DEGs) in the female and male D. macrosoma, 1,150 of which were female-biased DEGs, and 1,195 unigenes were male-biased DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the DEGs were mainly involved in biological processes including protein synthesis, growth, rhythmic processes, immune defense, and vitellogenesis. Then, we identified many growth- and sex-related genes, including Igf, Fabps, EF-hand family genes, Zp3, Zp4 and Vg. In addition, a total of 19,573 simple sequence repeats (SSRs) were screened and identified from the transcriptome sequences. The results of this study can provide valuable information on growth- and sex-related genes and facilitate further exploration of the molecular mechanism of sexual growth dimorphism.

Transcriptome profiling based on protein-protein networks provides a core set of genes for understanding blood immune response mechanisms against LPS stress in Amphioctopus fangsiao

Gram-negative bacteria are significant pathogens in the ocean, posing serious threats to marine organisms. Lipopolysaccharide (LPS) is a characteristic chemical constituent in Gram-negative bacteria that can be recognized by the pattern recognition receptor (PRR) of immune cells. This system is often used to simulate the invasion of bacteria. Blood is a transport channel for immune cells, and its transcriptome information obtained from Amphioctopus fangsiao stimulated by LPS is essential for understanding the antibacterial biological mechanisms of this species. In this study, we analyzed the gene expression profiles of A. fangsiao blood within 24h under LPS stress and found 778 and 561 differentially expressed genes (DEGs) at 6 and 24h, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses were performed to search for immune-related DEGs. The relationships among immune genes were examined by constructing a protein-protein interaction (PPI) network. Finally, 16 hub genes were identified based on the PPI network and KEGG enrichment analysis. The expression profiles of these genes were verified using quantitative RT-PCR (qRT-PCR). This research provides valuable resources for the healthy culture of A. fangsiao and helps us understand the molecular mechanisms of innate immunity.

The first transcriptome sequencing and data analysis of the Javan mahseer (Tor tambra)

The Javan mahseer (Tor tambra) is one of the most valuable freshwater fish found in Tor species. To date, other than mitogenomic data (BioProject: PRJNA422829), genomic and transcriptomic resources for this species are still lacking which is crucial to understand the molecular mechanisms associated with important traits such as growth, immune response, reproduction and sex determination. For the first time, we sequenced the transcriptome from a whole juvenile fish using Illumina NovaSEQ6000 generating raw paired-end reads. De novo transcriptome assembly generated a draft transcriptome (BUSCO5 completeness of 91.2% [Actinopterygii_odb10 database]) consisting of 259,403 putative transcripts with a total and N50 length of 333,881,215 bp and 2283 bp, respectively. A total count of 77,503 non-redundant protein coding sequences were predicted from the transcripts and used for functional annotation. We mapped the predicted proteins to 304 known KEGG pathways with signal transduction cluster having the highest representation followed by immune system and endocrine system. In addition, transcripts exhibiting significant similarity to previously published growth-and immune-related genes were identified which will facilitate future molecular breeding of Tor tambra.

De Novo Transcriptomic Characterization Enables Novel Microsatellite Identification and Marker Development in Betta splendens

The wild populations of the commercially valuable ornamental fish species, Betta splendens, and its germplasm resources have long been threatened by habitat degradation and contamination with artificially bred fish. Because of the lack of effective marker resources, population genetics research projects are severely hampered. To generate genetic data for developing polymorphic simple sequence repeat (SSR) markers and identifying functional genes, transcriptomic analysis was performed. Illumina paired-end sequencing yielded 105,505,486 clean reads, which were then de novo assembled into 69,836 unigenes. Of these, 35,751 were annotated in the non-redundant, EuKaryotic Orthologous Group, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes and Gene Ontology databases. A total of 12,751 SSR loci were identified from the transcripts and 7970 primer pairs were designed. One hundred primer pairs were randomly selected for PCR validation and 53 successfully generated target amplification products. Further validation demonstrated that 36% (n = 19) of the 53 amplified loci were polymorphic. These data could not only enrich the genetic information for the identification of functional genes but also effectively facilitate the development of SSR markers. Such knowledge would accelerate further studies on the genetic variation and evolution, comparative genomics, linkage mapping and molecular breeding in B. splendens.

Comparison of Gonadal Transcriptomes Uncovers Reproduction-Related Genes with Sexually Dimorphic Expression Patterns in Diodon hystrix

Diodon hystrix is a new and emerging aquaculture species in south China. However, due to the lack of understanding of reproductive regulation, the management of breeding and reproduction under captivity remains a barrier for the commercial aquaculture of D. hystrix. More genetic information is needed to identify genes critical for gonadal development. Here, the first gonadal transcriptomes of D. hystrix were analyzed and 151.89 million clean reads were generated. All reads were assembled into 57,077 unigenes, and 24,574 could be annotated. By comparing the gonad transcriptomes, 11,487 differentially expressed genes were obtained, of which 4599 were upregulated and 6888 were downregulated in the ovaries. Using enrichment analyses, many functional pathways were found to be associated with reproduction regulation. A set of sex-biased genes putatively involved in gonad development and gametogenesis were identified and their sexually dimorphic expression patterns were characterized. The detailed transcriptomic data provide a useful resource for further research on D. hystrix reproductive manipulation.

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.

Interleukin-34, a Novel Paracrine/Autocrine Factor in Mouse Testis, and Its Possible Role in the Development of Spermatogonial Cells In Vitro

Spermatogenesis is the process of spermatogonial stem cell (SSC) proliferation and differentiation to generate sperm. This process is regulated by cell–cell interactions between Sertoli cells and developing SSCs by autocrine/paracrine and endocrine factors. It is also affected by cells in the interstitial compartment, such as Leydig cells and peritubular cells. Here, we demonstrate, for the first time, the presence of interleukin-34 (IL-34) in Leydig, Sertoli, and peritubular cells and in the premeiotic, meiotic, and postmeiotic cells. Its receptor, colony-stimulating factor-1 (CSF-1), has already been demonstrated in Leydig, Sertoli, premeiotic, and meiotic cells. IL-34 was detected in testicular homogenates and Sertoli cell-conditioned media, and was affected by mouse age. We showed that the addition of IL-34 in vitro to isolated cells from the seminiferous tubules of 7-day-old mice, using the methylcellulose culture system (MCS), increased the percentages and expression of the premeiotic cells (VASA), the meiotic cells (BOULE), and the meiotic/postmeiotic cells (ACROSIN) after four weeks of culture, when examined by immunofluorescence staining (IF) and qPCR analysis. It is possible to suggest that IL-34 is a novel paracrine/autocrine factor involved in the development of spermatogenesis. This factor may be used in future therapeutic strategies for the treatment of male infertility.

Transcriptional profiling of skeletal muscle reveals starvation response and compensatory growth in Spinibarbus hollandi

BACKGROUND

Spinibarbus hollandi is an economically important fish species in southern China. This fish is known to have nutritional and medicinal properties; however, its farming is limited by its slow growth rate. In the present study, we observed that a compensatory growth phenomenon could be induced by adequate refeeding following 7 days of fasting in S. hollandi. To understand the starvation response and compensatory growth mechanisms in this fish, the muscle transcriptomes of S. hollandi under control, fasting, and refeeding conditions were profiled using next-generation sequencing (NGS) techniques.

RESULTS

More than 4.45 × 108 quality-filtered 150-base-pair Illumina reads were obtained from all nine muscle samples. De novo assemblies yielded a total of 156,735 unigenes, among which 142,918 (91.18%) could be annotated in at least one available database. After 7 days of fasting, 2422 differentially expressed genes were detected, including 1510 up-regulated genes and 912 down-regulated genes. Genes involved in fat, protein, and carbohydrate metabolism were significantly up-regulated, and genes associated with the cell cycle, DNA replication, and immune and cellular structures were inhibited during fasting. After refeeding, 84 up-regulated genes and 16 down-regulated genes were identified. Many genes encoding the components of myofibers were significantly up-regulated. Histological analysis of muscle verified the important role of muscle hypertrophy in compensatory growth.

CONCLUSION

In the present work, we reported the transcriptome profiles of S. hollandi muscle under different conditions. During fasting, the genes involved in the mobilization of stored energy were up-regulated, while the genes associated with growth were down-regulated. After refeeding, muscle hypertrophy contributed to the recovery of growth. The results of this study may help to elucidate the mechanisms underlying the starvation response and compensatory growth.

Comparative transcriptomics identifies genes differentially expressed in the intestine of a new fast-growing strain of common carp with higher unsaturated fatty acid content in muscle

We have created a new, fast-growing strain of common carp with higher unsaturated fatty acid content in muscle. To better understand the impacts of gene regulation in intestinal tissue on growth and unsaturated fatty acid content, we conducted a comparative RNA-Seq transcriptome analysis between intestine samples of Selected and Control groups (and corroborated selected results by PCR). After eight weeks of cage culture, weight gain of the Selected group was 20.84% higher. In muscles of the control group, monounsaturated fatty acids (FAs) were more abundant, whereas polyunsaturated FAs were more abundant in muscles of the Selected group. In total, we found 106 differentially expressed genes (DEGs) between the two groups. Only the endocytosis pathway was significantly enriched in DEGs, with two upregulated genes: il2rb and ehd1. The latter is involved in the growth hormone/insulin-like growth factor (Gh/Igf) axis, which plays a key role in the regulation of growth in animals. tll2, which is known to be associated with intestinal regeneration, was extremely highly upregulated in both transcriptomic (infinite) and qPCR (610.70) analyses. Two of the upregulated genes are associated with the fatty acid metabolism, several genes are likely to be indicators of heightened transcription levels, several are associated with metabolic and developmental roles, several with neuronal functions (including two with vision), several with the immune system, and two downregulated genes with the development of vasculature. The higher growth rate of the Selected group is likely to be at least partially attributed to increased endocytosis efficiency and genetically-driven behavioural differences (higher aggression levels). There are some indications that this new strain might have slightly impaired immune responses, and a higher propensity for inherited diseases leading to sight impairment, as well for neurodegenerative diseases in general, but these indications still need to be confirmed.

Transcriptome Dynamics During Turbot Spermatogenesis Predicting the Potential Key Genes Regulating Male Germ Cell Proliferation and Maturation

Spermatogenesis is a dynamic developmental process in which spermatogonial stem cells proliferate, differentiate and mature into functional spermatozoa. These processes require an accurate gene regulation network. Here, we investigated the dynamic changes that occur during spermatogenesis through a combination of histological and transcriptome analyses of different developmental stages of the testis. We constructed 18 testis transcriptome libraries, and the average length, N50, and GC content of the unigenes were 1,795 bp; 3,240 bp and 49.25%, respectively. Differentially expressed genes (DEGs) that were related to germ cell proliferation and maturation, such as NANOS3, RARs, KIFs, steroid hormone synthesis-related genes and receptor genes, were identified between pairs of testis at different developmental stages. Gene ontology annotation and pathway analyses were conducted on DEGs with specific expression patterns involved in the regulation of spermatogenesis. Nine important pathways such as steroid hormone biosynthesis related to spermatogenesis were identified. A total of 21 modules that ranged from 49 to 7,448 genes were designed by a weighted gene co-expression network analysis. Furthermore, a total of 83 candidate miRNA were identified by computational methods. Our study provides the first transcriptomic evidence for differences in gene expression between different developmental stages of spermatogenesis in turbot (Scophthalmus maximus).

Transcriptome analysis for identification of candidate genes related to sex determination and growth in Charybdis japonica

Charybdis japonica is an important cultured crab in China and it exhibits sex differences in their growth. Growth is an important economic trait that is controlled by many genes. In order to discover the growth-related regulatory mechanisms, whole-body transcriptomic sequencing and comparative analyses in different genders of C. japonica were conducted based on Illumina RNA-seq technology. In the present study, we obtained 62,879,204 and 60,226,334 clean reads in female and male libraries, respectively. 25,000,000 clean reads of every library were randomly selected and compared with Nt database to examine the possible contamination. Results showed that all clean reads were distributed among C. japonica or other species that were closely relative to this species, indicating no-pollution. De novo assembly was performed and a total of 32,543 and 44,174 unigenes were produced in female and male of C. japonica, respectively. Among all the unigenes, 12,591 and 14,455 unigenes of female and male crabs were annotated based on protein databases. Moreover, a total of 33,926 unigenes were found to contain ORFs and 52,839 SSRs were detected. The contrast between male and female C. japonica identifying 1939 unigenes were significantly differentially expressed. In addition, we specifically discussed some gene functions and pathways potentially associated with sex determination and growth. This is the first systematic report of whole transcriptome in C. japonica. The transcriptome information provides a basic resource for further studies on understanding the molecular basis of biological processes in C. japonica and other crustaceans.

The identification of growth, immune related genes and marker discovery through transcriptome in the yellow drum (Nibea albiflora)

Yellow drum (Nibea albiflora) is a commercially important marine fish, which is widely distributed in the coastal waters of China, Japan and Korea. Wild yellow drum resources have dramatically declined due to overfishing and ocean pollution. Genetic data can contribute to biodiversity conservation and protection. And molecular markers can play important roles in genetic breeding and aid in germplasm preservation in fish. In this study, 11 tissues (brain, heart, liver, kidney, muscle, head kidney, skin, fin, spleen, gonad and air bladder) were collected for pooled RNA sequencing. The unigenes were assembled using Trinity and EvidentialGene, and were then aligned to nr, nt, Swiss-Prot GO, KEGG, and KOG for annotation. Molecular markers (e.g. simple sequence repeat, SSR and single nucleotide polymorphism, SNP) were detected using MIcroSAtellite identification tool (MISA) and Genome Analysis Tool Kit (GATK). All clean reads were assembled into 109,209 transcripts, and 31,183 unigenes were generated after pruning and classifying, ranging from 201 to 19,857 bp in length (1230 bp in average), and 26,728 (85.7%) assembled unigenes had significant hits in public databases. Total of 27 and 103 unigenes were respectively identified as involved in growth- and immune-related pathways in the N. albiflora transcriptome. In addition, we identified a considerable quantity of molecular markers, including 11,484 SSRs and 56,186 SNPs. The growth- and immune-relevant genes and the molecular markers identified here provided a meaningful reference gene set and laid a foundation for future genetic selection and breeding for this species.

De novoassembly ofSchizothorax waltonitranscriptome to identify immune-related genes and microsatellite markers

Schizothorax waltoni (S. waltoni) is one kind of the subfamily Schizothoracinae and an indigenous economic tetraploid fish to Tibet in China. It is rated as a vulnerable species in the Red List of China's Vertebrates, owing to overexploitation and biological invasion. S. waltoni plays an important role in ecology and local fishery economy, but little information is known about genetic diversity, local adaptation, immune system and so on. Functional gene identification and molecular marker development are the first and essential step for the following biological function and genetics studies. For this purpose, the transcriptome from pooled tissues of three adult S. waltoni was sequenced and analyzed. Using paired-end reads from the Illumina Hiseq4000 platform, 83 103 transcripts with an N50 length of 2337 bp were assembled, which could be further clustered into 66 975 unigenes with an N50 length of 2087 bp. The majority of the unigenes (58 934, 87.99%) were successfully annotated by 7 public databases, and 15 KEGG pathways of immune-related genes were identified for the following functional research. Furthermore, 19 497 putative simple sequence repeats (SSRs) of 1–6 bp unit length were detected from 14 690 unigenes (21.93%) with an average distribution density of 1 : 3.28 kb. We identified 3590 unigenes (5.36%) containing more than one SSR, providing abundant potential polymorphic markers in functional genes. This is the first reported high-throughput transcriptome analysis of S. waltoni, and it would provide valuable genetic resources for the functional genes involved in multiple biological processes, including the immune system, genetic conservation, and molecular marker-assisted breeding of S. waltoni.

De novo assembly of Schizothorax waltoni transcriptome.

Transcriptome profiling and expression analysis of immune responsive genes in the liver of Golden mahseer (Tor putitora) challenged with Aeromonas hydrophila

Transcriptome profiling has been used to decipher the novel mechanisms behind immune responses of the fishes. However, the molecular mechanism underlining immune response in mahseer is not studied so far. Fishes are greatly affected by bacterial pathogens such as Aeromonas hydrophila. In this study, transcriptome response of golden mahseer (Tor putitora) infected with A. hydrophila was examined using paired end Illumina sequencing of liver tissue to understand the immune response of the fish. The de novo assembly generated 61,042 unigenes ranging from 200 to 9322 bp in length and an average length of 463 bp. The gene ontology annotations resulted a total of 131,826 term assignments to the annotated transcriptome including 60,846 (46.16%) allocations from the biological process; 21,603 (16.39%) from molecular function and 49,377 (37.46%) from cellular components. Differential gene expression analysis of the transcriptome data from challenged and control group revealed 1104 upregulated and 1304 down-regulated unigenes. The differentially expressed genes were mainly involved in the pathways including cell surface receptor signaling, TH1 and TH2 cell differentiation, pathogen recognition, and immune system process/defense response especially complement cascade. Twelve unigenes including ankyrin, serum amyloid, hsp4b, STAT3, complement factor c3 and c7 were validated using qPCR and found differentially expressed in accordance with in silico expression analysis. The results obtained in this study will provide the first and crucial information on the molecular mechanism of mahseer fishes against bacterial infection.

Turbot (Scophthalmus maximus) vs. VHSV (Viral Hemorrhagic Septicemia Virus): A Review

Turbot (Scophthalmus maximus) is a very valuable fish species both in Europe and China. The culture of this flatfish is well-established but several bacteria, viruses, and parasites can produce mortality or morbidity episodes in turbot farms. Viral Hemorrhagic Septicemia Virus (VHSV) is one of the most threatening pathogens affecting turbot, because neither vaccines nor treatments are commercially available. Although the mortality in the turbot farms is relatively low, when this virus is detected all the stock have to be destroyed. The main goals that need to be improved in order to reduce the incidence of this disease is to know what are the strategies or molecules the host use to fight the virus and, in consequence, try to potentiate this response using different ways. Certain molecules can be selected as potential antiviral treatments because of their high protective effect against VHSV. On the other hand, the use of resistance markers for selective breeding is one of the most attractive approaches. This review englobes all the investigation concerning the immune interaction between turbot and VHSV, which until the last years was very scarce, and the knowledge about VHSV-resistance markers in turbot. Nowadays, the availability of abundant transcriptomic information and the recent sequencing of the turbot genome open the door to a more exhaustive and profuse investigation in these areas.