Transcriptome analysis provides insights into differentially expressed genes and long noncoding RNAs involved in sex‐related differences in Amur sturgeon ( Acipenser schrenckii )

Transcriptome Analyses Reveal Expression Profiles of Morphologically Undifferentiated and Differentiated Gonads of Yangtze Sturgeon Acipenser dabryanus

Sturgeon is known as a primitive fish with the ZZ/ZW sex determination system and is highly prized for its valuable caviar. Exploring the molecular mechanisms underlying gonadal differentiation would contribute to broadening our knowledge on the genetic regulation of sex differentiation of fish, enabling improved artificial breeding and management of sturgeons. However, the mechanisms are still poorly understood in sturgeons. This study aimed to profile expression patterns between female and male gonads at morphologically undifferentiated and early differentiated stages and identify vital genes involved in gonadal sex differentiation of sturgeons. The sexes of Yangtze sturgeon (Acipenser dabryanus) juveniles were identified via the sex-specific DNA marker and histological observation. Transcriptome analyses were carried out on female and male gonads at 30, 80 and 180 days post-hatching. The results showed that there was a total of 17 overlapped DEGs in the comparison groups of between female and male gonads at the three developmental stages, in which there were three DEGs related to ovarian steroidogenesis, including hsd17b1, foxl2 and cyp19a1. The three DEGs were highly expressed in the female gonads, of which the expression levels were gradually increased with the number of days after hatching. No well-known testis-related genes were found in the overlapped DEGs. Additionally, the expression levels of hsd17b1 and cyp19a1 mRNA were decreased with the knockdown of foxl2 mRNA via siRNA. The results further suggested that foxl2 should play a crucial role in the ovarian differentiation of sturgeons. In conclusion, this study showed that more genes involved in ovarian development than testis development emerged with sexually dimorphic expression during early gonadal sex differentiation, and it provided a preliminary understanding of the molecular regulation on gonadal differentiation of sturgeons.

The potential role of eyestalk in the immunity of Litopenaeus vannamei to Vibrio parahaemolyticus infection II. From the perspective of long non-coding RNA

Long non-coding RNAs (lncRNAs) have been linked to immunological modulation. Unfortunately, little is known about the processes of immune control in shrimp. In crustaceans such as Litopenaeus vannamei, a prominent aquaculture species, the X-organ-sinus gland complex (XO-SG) in the eyestalk is an essential neuroendocrine regulatory organ. Eyestalk ablation is commonly employed in aquaculture to accelerate ovarian maturation in shrimp. It does, however, have a negative impact on the shrimps' immunocompetence and causes death. As a result, we used RNA-seq to profile the transcriptomes of L. vannamei hemocytes infected with Vibrio parahaemolyticus after the eyestalk ablation. Following strict transcript screening procedures, 2307 lncRNAs were identified from L. vannamei hemocytes in this study. Pearson correlation analysis was finally used to uncover 535 DElncRNAs and 1566 DEmRNA targets. According to the Venn diagram analysis, 326 non-eyestalk regulatory lncRNAs (NElncRNAs) with a target of 1014 non-eyestalk regulatory genes (NEmRNAs), 47 eyestalk negative regulatory lncRNAs (ENRlncRNAs) with a target of 95 eyestalk negative regulatory genes (ENRmRNAs), and 162 eyestalk positive regulatory lncRNAs (EPRlncRNAs) with a target of 457 eyestalk positive regulatory genes (EPRmRNAs) were screened. The bioinformatics analysis revealed that lncRNAs were associated with Axon regeneration, Rap1 signaling pathway, Thyroid hormone signaling pathway, TGF-beta signaling pathway, and PI3K-Akt signaling pathway, implying that lncRNAs may play a role in the regulation of the neuroendocrine-immune (NEI) system. Furthermore, several lncRNAs targeting HSP70, YWHAZ, FER2, HIF1α, and Notch were discovered and verified by qRT-PCR. These findings showed that regulation of lncRNAs in hemocytes which were controlled by the eyestalk might be one of the impact variables in controlling the differential expression of mRNAs associated with immune response in L. vannamei infected with V. parahaemolyticus.

Differentially expressed lncRNAs involved in immune responses of Haliotis diversicolor and H. discus hannai challenged with Vibrio parahaemolyticus

Although many studies have shown that lncRNA, a non-coding RNA with a length of more than 200 bases, is involved in various biological functions, including the immune process, stress process, and cell development process. However, the function of lncRNA in abalone, especially in immunity, has been rarely studied. H. discus hannai and H. diversicolor are two main aquaculture abalone, and their growth is easily affected by the main pathogen Vibrio parahaemolyticus. Through rigorous screening procedures for transcripts in this study, we found that lncRNAs were 34,240, 23,022 in Haliotis diversicolor and H. discus hannai injected with V. parahaemolyticus, respectively. We also identified the unique and common lncRNAs and mRNAs of two abalone species for the first time; the shared lncRNAs and mRNAs in Haliotis diversicolor and H. discus hannai were 2352 and 13,165, respectively. Then gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the differentially expressed target genes of common and unique lncRNAs has shown that common lncRNAs could be widely involved in the biological processes of stress and cell development in both abalone species. In contrast, unique lncRNAs are linked to the Toll-like receptor, NF-kappaB signaling pathway of H. diversicolor, and pattern recognition receptors and lectins immune-related pathways of H. discus hannai. The co-expression network shows that some immune-related genes, such as INFK1, INFK2, CASP2, CASP8, IRAK1, lectin C, were closely related to lncRNAs. Further, we identified the targeted relationship between some immune-related genes and lncRNAs by qRT-PCR, through which we showed that the expression trend between targeted genes, such as INFK1 and Lnc7057, lectin C and Lnc6943, Lnc5637, and PLCG1 and Lnc1692, were consistent. In general, our results showed that lncRNA expression was induced in the two species of abalone after being infected with V. parahaemolyticus, and lncRNA was involved in the immune response of abalone by targeting coding genes.

Siberian sturgeon multi-tissue reference transcriptome database

Motivation: Siberian sturgeon is a long lived and late maturing fish farmed for caviar production in 50 countries. Functional genomics enable to find genes of interest for fish farming. In the absence of a reference genome, a reference transcriptome is very useful for sequencing based functional studies.

Results: We present here a high-quality transcriptome assembly database built using RNA-seq reads coming from brain, pituitary, gonadal, liver, stomach, kidney, anterior kidney, heart, embryonic and pre-larval tissues. It will facilitate crucial research on topics such as puberty, reproduction, growth, food intake and immunology. This database represents a major contribution to the publicly available sturgeon transcriptome reference datasets.

Availability: The database is publicly available at http://siberiansturgeontissuedb.sigenae.org

Supplementary information:  Supplementary data are available at Database online.

Full-length transcriptome sequencing and comparative transcriptomic analysis to uncover genes involved in early gametogenesis in the gonads of Amur sturgeon (Acipenser schrenckii)

BACKGROUND

Sturgeons (Acipenseriformes) are polyploid chondrostean fish that constitute an important model species for studying development and evolution in vertebrates. To better understand the mechanisms of reproduction regulation in sturgeon, this study combined PacBio isoform sequencing (Iso-Seq) with Illumina short-read RNA-seq methods to discover full-length genes involved in early gametogenesis of the Amur sturgeon, Acipenser schrenckii.

RESULTS

A total of 50.04 G subread bases were generated from two SMRT cells, and herein 164,618 nonredundant full-length transcripts (unigenes) were produced with an average length of 2782 bp from gonad tissues (three testes and four ovaries) from seven 3-year-old A. schrenckii individuals. The number of ovary-specific expressed unigenes was greater than those of testis (19,716 vs. 3028), and completely different KEGG pathways were significantly enriched between the ovary-biased and testis-biased DEUs. Importantly, 60 early gametogenesis-related genes (involving 755 unigenes) were successfully identified, and exactly 50% (30/60) genes of those showed significantly differential expression in testes and ovaries. Among these, the Amh and Gsdf with testis-biased expression, and the Foxl2 and Cyp19a with ovary-biased expression strongly suggested the important regulatory roles in spermatogenesis and oogenesis of A. schrenckii, respectively. We also found the four novel Sox9 transcript variants, which increase the numbers of regulatory genes and imply function complexity in early gametogenesis. Finally, a total of 236,672 AS events (involving 36,522 unigenes) were detected, and 10,556 putative long noncoding RNAs (lncRNAs) and 4339 predicted transcript factors (TFs) were also respectively identified, which were all significantly associated with the early gametogenesis of A. schrenckii.

CONCLUSIONS

Overall, our results provide new genetic resources of full-length transcription data and information as a genomic-level reference for sturgeon. Crucially, we explored the comprehensive genetic characteristics that differ between the testes and ovaries of A. schrenckii in the early gametogenesis stage, which could provide candidate genes and theoretical basis for further the mechanisms of reproduction regulation of sturgeon.