- Research article
- Open Access
RNA-sequencing of the sturgeon Acipenser baeri provides insights into expression dynamics of morphogenic differentiation and developmental regulatory genes in early versus late developmental stages
© The Author(s). 2016
- Received: 17 December 2015
- Accepted: 15 June 2016
- Published: 8 August 2016
Acipenser baeri, one of the critically endangered animals on the verge of extinction, is a key species for evolutionary, developmental, physiology and conservation studies and a standout amongst the most important food products worldwide. Though the transcriptome of the early development of A. baeri has been published recently, the transcriptome changes occurring in the transition from embryonic to late stages are still unknown. The aim of this work was to analyze the transcriptomes of embryonic and post-embryonic stages of A. baeri and identify differentially expressed genes (DEGs) and their expression patterns using mRNA collected from specimens at big yolk plug, wide neural plate and 64 day old sturgeon developmental stages for RNA-Seq.
The paired-end sequencing of the transcriptome of samples of A. baeri collected at two early (big yolk plug (T1, 32 h after fertilization) and wide neural plate formation (T2, 45 h after fertilization)) and one late (T22, 64 day old sturgeon) developmental stages using Illumina Hiseq2000 platform generated 64039846, 64635214 and 75293762 clean paired-end reads for T1, T2 and T22, respectively. After quality control, the sequencing reads were de novo assembled to generate a set of 149,265 unigenes with N50 value of 1277 bp. Functional annotation indicated that a substantial number of these unigenes had significant similarity with proteins in public databases. Differential expression profiling allowed the identification of 2789, 12,819 and 10,824 DEGs from the respective T1 vs. T2, T1 vs. T22 and T2 vs. T22 comparisons. High correlation of DEGs’ features was recorded among early stages while significant divergences were observed when comparing the late stage with early stages. GO and KEGG enrichment analyses revealed the biological processes, cellular component, molecular functions and metabolic pathways associated with identified DEGs. The qRT-PCR performed for candidate genes in specimens confirmed the validity of the RNA-seq data.
This study presents, for the first time, an extensive overview of RNA-Seq based characterization of the early and post-embryonic developmental transcriptomes of A. baeri and provided 149,265 gene sequences that will be potentially valuable for future molecular and genetic studies in A. baeri.
- Acipenser baeri
- Transcriptome sequencing
- Differentially expressed genes
- Embryonic stages
- Late stages
During metazoan embryonic development, a totipotent zygote divides, grows and experiences intense post-embryonic anatomical and physiological changes resulting in an adult living organism constituted of different specific tissues . This process is naturally correlated with spatial or temporal changes in gene expression. Therefore, complete assessment of gene expression levels during ontogeny of an organism would be fundamental for giving a comprehensive dataset for genetic regulation of developmental processes. Up to date, a large body of literature has reported the use of in situ hybridization, microarrays, and more latterly transcriptome sequencing (RNA-seq) technologies for profiling gene expression during development of various model and non-model organisms [1–6]. Sturgeons, one of the earliest origins of vertebrate groups, constitute an important archetypal material for studying the origin of species and evolution [7–9]. Moreover, sturgeons are listed in the appendix of the endangered species by the Convention on International Trade in Endangered Species of wild fauna and flora (CITES) . However, scientific studies on molecular mechanisms controlling the development of sturgeons are scarce and generally focused on a narrow range of single genes or gene sets. Transcriptomics resources for sturgeons have emerged just recently and, until now, data of sturgeon transcriptome include those made available owing to recent studies on reproductive tissues [11, 12] obtained by next-generation pyrosequencing of gonad transcriptomes of Acipenser fulvescens, de novo assembly of the gonadal transcriptome of Acipenser sinensis and microRNA transcriptome and expression assay in Acipenser schrenckii . In addition, we have recently made available the largest sturgeons’ transcriptomics data using RNA-sequencing (RNA-seq) to generate the transcriptome for the early development of A. baeri . Nevertheless, little is known about late developmental stages of A. baeri and in regards to the molecular background concerning the transition from pre-larval to juvenile stages, even less information has been made accessible, thus hindering aquaculture practices for this species. Studying this undiscovered molecular areas of A. baeri’s developmental biology, especially gene regulation underlying the transformation of embryos into adult fish, would be vital for feeding, reproductive and fish health management purposes, and would give insights into the biology of sturgeons and other related fish species.
The aim of the present study was to assess the transcriptome and the gene expression dynamics of three developmental stages of A. baeri ranging from the embryonic up to the 64 days old sturgeon stages.
Sequencing data quality assessment and de novo assembly
Statistical results of raw and preprocessed sequences
Functional annotation of unigenes
To determine the function of de novo assembled transcripts, the whole set of sequences were aligned against the NCBI Uniprot protein databases using BLASTX with an E-value cutoff of 1E-3. The significant alignment results are reported in Additional file 3. The result showed that 57,346 unigenes (38.42 %) had noteworthy hits to Uniprot databases equivalent to 45,837 single known proteins and 11,509 homologous orthology clusters in Uniprot protein databases whereas the remaining 61.58 % unigenes could stand for UTRs, non-protein coding genes or A. baeri-specific genes which were too different to be annotated by homology search with the adopted E-value cutoff. We performed analysis of BLASTX results to determine best species hits (Additional file 4). The results showed the hits of A. baeri’s transcripts with 1312 distinct species including fish species such as Latimeria chalumnae (5940 transcripts), Danio rerio (4943 transcripts), Helobdella robusta (4591 transcripts) and Oreochromis niloticus (2325 transcripts) and in minor extent to other vertebrate species.
The GO annotation was performed by the mapping of unigene sequences with Uniprot database and the recovery of GO terms linked to protein sequences in Uniprot. The results (Additional file 5) showed that a total of 43,062 unigenes (28.85 %) were ascribed at least one GO term in the GO classes of “cellular component”, “biological process” and “molecular function”.
In the KEGG database (Additional file 6), 29,526 unigenes (19.78 %) were annotated into 329 pathways. “Metabolic pathways” (2010 transcripts) was the major pathway and was accompanied by “biosynthesis of secondary metabolites” (642 transcripts) and “microbial metabolism in diverse environments” (417 transcripts). In the pathway class of signal transduction, “PI3K-Akt signaling pathway” (347 transcripts) was the most exemplified. Other pathways related to translation and diseases were similarly present. Overall, 27,773 unigenes were annotated in both GO and KEGG databases.
Interproscan annotation (Additional file 7) was exploited to detect conserved domains associated with the protein sequences. A total of 16,735 unigenes (11.21 %) were annotated in Interproscan and produced 6751 domains. The statistics (Additional file 8) revealed “zinc finger, C2H2” as the main conserved domain (count = 3392) followed by “WD40 repeat” (count = 2312) and “zinc finger, C2H2-like” (count = 1799).
Unigene differential expression profiling
GO and KEGG pathway enrichment analysis of DEGs within contiguous development stages
In the KEGG pathway enrichment analysis (Additional file 14, Fig. 5b), the pathway class of “cell growth and death” with 12 downregulated unigenes participating in cell cycle was the most represented. Asthma (2 downregulated unigenes), leishmaniasis (1 upregulated and 4 downregulated unigenes) and ovarian steroidogenesis (1 upregulated and 3 downregulated unigenes) were the most significantly (p-value < 0.05) enriched pathways following cell cycle.
Study case 1: differences among morphogenic differentiation regulatory genes
In order to identify transcripts implicated in morphogenesis, we screened the GO terms by searching for terms in relation with morphogenesis using the GO enrichment file. On this basis, we found 260 DEGs associated with morphogenesis in the category of biological process. The GO terms associated with these DEGs in the three pairwise comparisons and their expression profiles are summarized in Additional file 19. Morphogenic differentiation DEGs screened between T1 and T2 were associated with the GO terms of anatomical structure morphogenesis, anatomical structure formation involved in morphogenesis, tissue morphogenesis, embryonic morphogenesis and morphogenesis of follicular epithelium. We observed that 38 of these DEGs including ACTA1, twn-A, AANF, buc, dlx5a, NKD1, LOC770168, LOC101155136, APLNR, APOD, tbx16, LOC101164950, sox2 and cbs were identified as up-regulated unigenes while 35 unigenes including bbs9, DUSP4, CABZ01041002.1, CDC42SE2, Ndp, robo4, robo4, ERRFI, LOC100722730, epha7, D623_10029598, epb41l5, LAMA1 and CTNNB1 were downregulated. In the transition from T1 to T22, ACTA1, twn-A, AANF, buc, dlx5a, NKD1, LOC770168, LOC101155136, APLNR, APOD, tbx16, LOC101164950 and sox2 were found among the upregulated DEGs while the 166 downregulated DEGs included robo4, ERRFI, LOC100722730, epha7, D623_10029598, epb41l5, LAMA1 and CTNNB1. These DEGs were associated with anatomical structure morphogenesis, embryonic morphogenesis, anatomical structure formation involved in morphogenesis, tissue morphogenesis, morphogenesis of an epithelium and blood vessel morphogenesis. The transition from T2 to T22 was characterized by the upregulation of 31 and the downregulation of 131 DEGs encoding for biological processes such as anatomical structure morphogenesis, blood vessel morphogenesis, anatomical structure formation involved in morphogenesis, cellular component morphogenesis, cell morphogenesis, embryonic morphogenesis, regulation of cell morphogenesis and regulation of anatomical structure morphogenesis. The highest number of morphogenic differentiation DEGs (213) was found in T1 vs. T22 comparison while the lowest one (73) was found between T1 and T2. The T2 vs. T22 comparison allowed the identification of 162 DEGs.
Study case 2: changes among developmental regulatory genes
The definitive objective of this study was to identify an assortment of potential genes involved in developmental processes. Through analysis of GO enrichment output file, we found that 517 DEGs were directly implicated in developmental processes (Additional file 20). In T1 to T2 comparison, 151 DEGs were enriched in 22 developmental process including embryo development, organ development, system development and tissue development with nos3, Schip1, fbn2b, LOC101073546 and LOC100546827 as the most upregulated DEGs while the most downregulated DEGs included birc5b, buc, LOC101155136, nanog, LSM14B, lft1 and LOC101064479. In the T1 vs. T22 analysis, 35 enriched GO terms (biological processes) were orchestrated by the expression of 416 DEGs. Genes POSTN, f3a, LOC101073546, SCEL, COL1A3 and gfra3 were found as the most upregulated whereas fgfrl1b, twn-A, lft1, nanog, buc, LOC101064479 and tbx16 were significantly downregulated. In the comparison of T2 and T22, 80 upregulated DEGs including COL1A3, HG2A, LOC100534402, MBP, NEBL and alas2, and 237 down-regulated DEGs containing nanog, FOXP1B, LOC100709614, alcam,fgfrl1b, twn-A, RIPP1, buc and tbx16 were in charge of 28 biological processes. Similarities between T1 vs. T22 and T2 vs. T22 were observed, which further confirmed that T2 presented intermediary transcriptional changes.
Candidate gene qRT-PCR validation
To validate the sequencing data, we randomly chose five unigenes, namely six3a, sox17, HOXD10, wnt11b and eya3 genes, to perform the qRT-PCR experiment using the same pooled RNA employed for generating our RNA-seq data. The expression levels of each unigene obtained by RNA-seq or RT-PCR are presented in Additional file 21. In the sequencing result, we found that six3a was upregulated in T2 vs. T1 and T22 vs. T2 comparisons. HOXD10 as well as sox17 were upregulated in T2 compared to other samples. Initially highly expressed in T1, wnt11b and eya3 were progressively downregulated until T22. The trends of expression of all five candidate genes, measured by qRT-PCR, were correlated with those of RNA-seq method. Altogether, qRT-PCR results largely supported RNA-seq results.
With the advent of next-generation sequencing and the development of bioinformatics systems , RNA sequencing, which is preferred compared with the conventional Sanger sequencing due to its low cost and high-throughput generation of quality transcriptome data, has turned into an essential instrument used in research and relatively short reads can be successfully assembled for non-model organisms [17–20]. In this study, we have produced the first broad map of A. baeri transcriptome using an Illumina paired end RNA-seq platform for inspecting gene expression dynamics in early and late developmental stages. Around 64039846, 64635214 and 75293762 high quality reads were generated respectively for T1, T2 and T22 from HiSeq 2000 and de novo assembled into 149,265 unigenes. The number of unigenes was higher than that reported in the testicular and ovarian transcriptomes of A. schrenckii containing respectively 122,381 and 114,527 unigenes  and that of the gonads transcriptome of A. sinensis constituted of 86,027 unigenes  or the 55,000 high quality ESTs of A. naccarii organized into a freely available AnaccariiBase [11, 23]. The higher number of unigenes generated here could be due to the fact that we performed sequencing using whole bodies of A. baeri samples, a deep sequencing coverage or differences in assembly softwares used. In addition, we remarked the increase of unigene amount in the transition from early developmental stages to subsequent late stages probably because of the activation of additional unigenes necessary for achieving new biological processes that arise harmoniously with the developmental course. This was in accordance with the fact that initial stages of vertebrate development depend dominatingly on maternal factors deposited in the egg with negligible zygotic translation until the complete activation of the embryonic genetic material [24, 25]. The validity of de novo assembly was assessed by examination of assembled unigenes using publicly available protein databases, functional annotation and validation of randomly selected unigenes by the qRT-PCR. Annotation of unigenes showed that 57,346 unigenes (38.42 %) had noteworthy homologs in the Uniprot databases while 28.85 % of unigenes were ascribed at least one GO term in the GO classes of “cellular component”, “biological process” and “molecular function”. In the KEGG database, 19.78 % of unigenes were annotated into 329 pathways. The Interproscan annotation produced 6751 domains corresponding to 16,735 unigenes (11.21 %). Though various technologically advanced de novo assemblers (e.g. Trans-AbySS, Oasis, SOAP2denovo and Trinity) have been established, scientific studies have proven the efficacy of Trinity since it is advantageous for annotation of transcriptomes for diverse vertebrate species [19, 26]. Here, a substantial amount of the Trinity de novo assembled unigenes was annotated in known protein databases. The unigenes without hits presumably fit in with untranslated regions, non-coding RNA, novel genes, short sequences without protein domain or assembly errors, and suggests that the A. baeri genome is still incomplete. Our study will help in the future completion and full annotation of the A. baeri’s genome.
We equally examined the differential expression of unigenes and identified 2789, 12,819 and 10,824 DEGs from T1 vs. T2, T1 vs. T22 and T2 vs. T22 comparisons according to our RNA-seq data and high correlation was found between T1 and T2 while no significant correlation was recorded for T1 vs. T22 and T2 vs. T22. This result showed that the DEGs match expectations in that the fewest number of differences were observed between the two early developmental stages (T1 vs. T2), and the greatest number was observed between T1 and T22, while T2 vs. T22 differences were intermediate. Based on GO enrichment analysis of DEGs, we identified 260 DEGs involved in morphogenic differentiation and 517 DEGs encoding for diverse developmental processes. Besides, the GO and KEGG pathway functional annotation of the overall set of DEGs allowed identification of multiple functions and pathways in which DEGs between samples were involved. These data give insight in the gene expression alterations occurring in the transition from the early stage to the late stage of A.baeri development.
Although the RNA-Seq experiment was performed using pooled RNA extracted from 3 specimens per development stage, the dataset presents some limitations. There were no biological or technical replicates for the pairwise comparison between the 3 samples. The dataset is therefore un-replicated and ‘sample’, but not really ‘development-stage’ specific. Consequently, although the statistical analysis of many transcripts yields statistically significant changes between these 3 pooled samples, the biological replicate is still N = 1, so there is no variation in gene expression addressed. In addition, the present RNASeq data is a transcriptome screen, not necessarily a “functional assay” and the fact that we did not perform some translational experiments could mislead our interpretation because the gene function is not necessarily determined using RNA-Seq given that, while it is regarded and experimentally proven that instances of up- or down-levels of expression generally translate (no pun intended) into higher or lower protein levels, respectively, this cannot be automatically assumed without further biochemical analyses. Post-translational modifications do occur that could alter anticipated final protein levels and hence influence pathway interactions.
Our research established, for the first time, an extensive overview of RNA-Seq based characterization of the early and post-feeding developmental transcriptome of A. baeri as well as significant data on differential expression among both late and early developmental stages. The RNA-seq provided considerable gene sequences that will be valuable for future molecular and genetic studies in A. baeri and other related sturgeon and fish species.
Fish specimens, mRNA extraction and Illumina sequencing
Main characteristics of studied specimens defining each developmental stage
Main morphological characteristics
Big yolk plug
The formation of large yolk plug:
About 32 h after fertilization, the embryo develops into yolk plug period: blastopore endodermal cells are still externally exposed at the vegetal pole, formation of a large yolk plug inserted in the germ ring, at this time the animal is very bright yellow and the plant polar pigmentation is very deep.
Wide neural plate
Obvious wide neural plate divided into both internal and external parts: about 45 h after fertilization, there is apparition of an obvious wide neural plate; around the neural plate in the head region, there is a clear shaped horseshoe formation; the neural plate grows longitudinally and is divided into inner and outer parts; in the center of the neural plate grows a longitudinal neural groove.
64-day-old fish with electrical receptors ganglion
Fully developed bone plate, 14–16 dorsal spines, 45–47 lateral bone plates, 7–9 abdominal bone plates. Presence of ampullary organs on the ventral face
Sequencing data quality assessment and de novo assembly
Prior to assembly and mapping, we applied filters for quality control of sequenced reads. Trim Galore software was adopted for quality trimming of raw reads and dynamic removal of adapters and low-quality fragments.
High-quality sequence reads stemming from quality control analysis were De novo assembled into transcripts using the Trinity platform (http://trinityrnaseq.sf.net) including Inchworm, Chrysalis and Butterfly as independent modules following the protocol described elsewhere . The evaluation of the efficacy of the assembly was performed taking into account the total number of transcripts, transcripts length distribution and N50 value. We calculated the N50 size according to the threshold of lengths of different transcripts and counted transcripts with lengths greater than or equal to the minimum threshold.
A unigene was defined as the longest transcript among the multitude of assembly transcripts isomers. The optimization of assembly led to different transcript isomers (isoforms) or paralogs. Each unigene (expressed in prefix comp + digital ID) corresponded to one or numerous transcripts isomers (comp*_c*_Seq*).
To discriminate between valid transcript sequences and incorrectly assembled sequences, we used TransDecoder program integrated in Trinity software (log likelihood ratio based on the ratio of coding to noncoding sequences) to extract open reading frames (ORFs) and predict potential protein coding domain sequences (CDS) according to Markov model principle. CDS were translated into amino acid sequences according to the standard codon table in order to obtain potential protein sequences coded by the transcripts.
Functional annotation of unigenes
De novo functional annotation of A.baeri transcriptome was obtained by similarity search against UniProt protein databases (Swissprot or Tremble) using BLASTX. Alignments with an E-value cut off of 1E-3 were considered significant and gene annotation information was assigned to transcripts based on the highest BLAST hit. The Blast2GO suite  was used for the Gene Ontology (GO) annotation of unigenes. Annotation via Blast2GO was done by first searching for matches to the Uniprot databases, then mapping the BLAST results to the GO database and finally retrieving GO annotation information corresponding to transcripts with BLAST hits. The WEGO software was applied for classifying and counting GO classes. Additionally, the KEGG annotation of unigenes was achieved using the online KEGG database (http://www.genome.jp/kegg/). Using the HMM algorithm, Interproscan (http://www.ebi.ac.uk/InterProScan/), including PRINTS, SMART, Pfam, Coils, SUPERFAMILY, Gene3D, ProSiteProfiles, Hamap, ProSitePatterns, TIGRFAM and PIRSF databases, was performed for searching for protein domains.
Unigenes abundance estimation
To compute abundance estimates of transcripts, the original reads were firstly aligned to the Trinity transcripts. Subsequently, the RNA-Seq by Expectation-Maximization (RSEM) software  (default parameter Settings) was implemented to determine the expression levels of transcripts or corresponding unigenes in FPKM (expected number of fragments per kilobase of transcript sequence per millions base pairs sequenced) units.
Identification and functional enrichment analysis of differentially expressed genes
The Bioconductor tool edgeR  (p-value < 0.001, FDR <0.001) was exerted for screening differentially expressed genes (DEGs) by pairwise comparison between samples.
Note: N stands for the number of genes with GO annotation among all genes, n is the number of DEGs in N, M is the number of all genes that are annotated to a certain GO term; and m is the number of DEGs in M.
The method used for KEGG enrichment analysis was similar to that used in the GO enrichment analysis. Here, N is the number of genes with a KEGG annotation, n is the number of DEGs in N, M is the number of genes assigned to a specific pathway, and m is the number of DEGs in M. Pathways with a P-value < 0.05 were defined as significantly enriched pathways.
Quantitative real-time PCR analysis
List of primers designed for qRT-PCR validation of RNA-seq data
Region Y-Box 17 (sox17)
SIX homeobox 3a (six3a)
homeobox D10 (HOXD10)
wingless-type MMTV integration site family, member 11B (wnt11b)
EYA transcriptional coactivator and phosphatase 3 (eya3)
actin, beta (Actb)
We show appreciation to Shanghai YingBiotech Company for the dynamic help provided during the analysis of the sequencing data.
This work was supported by the National Natural Science Foundation of China (No. 31302161) and the National Infrastructure of Fishery Germplasm Resources.
Availability of data and materials
The sequencing data were deposited in the NCBI Short Read Archive (SRA) database (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/sra/) under the accession number SRP053165. All data of performed analyses are included as Additional files
WS, LM, KJ conceived and designed the experiments, WS, KJ, FZ, YL performed the experiments, WS, LM, KJ, FZ, YL analyzed the data, WS, LM, KJ, FZ, YL contributed in reagents/materials/analysis tools and WS, KJ, LM wrote the manuscript. All authors have read and approved the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was approved by the Review Committee for the Use of Animal Subjects of Shanghai Ocean University without requirement for particular permits.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Schep AN, Adryan B. A comparative analysis of transcription factor expression during metazoan embryonic development. PLoS One. 2013;8(6):e66826.View ArticlePubMedPubMed CentralGoogle Scholar
- Meyer E, Aglyamova GV, Wang S, Buchanan-Carter J, Abrego D, Colbourne JK, Willis BL, Matz MV. Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Genomics. 2009;10(1):1.View ArticleGoogle Scholar
- Zeng V, Villanueva KE, Ewen-Campen BS, Alwes F, Browne WE, Extavour CG. De novo assembly and characterization of a maternal and developmental transcriptome for the emerging model crustacean Parhyale hawaiensis. BMC Genomics. 2011;12(1):1.View ArticleGoogle Scholar
- Grün D, Kirchner M, Thierfelder N, Stoeckius M, Selbach M, Rajewsky N. Conservation of mRNA and protein expression during development of C. elegans. Cell Rep. 2014;6(3):565–77.View ArticlePubMedGoogle Scholar
- Brekhman V, Malik A, Haas B, Sher N, Lotan T. Transcriptome profiling of the dynamic life cycle of the scypohozoan jellyfish Aurelia aurita. BMC Genomics. 2015;16(1):74.View ArticlePubMedPubMed CentralGoogle Scholar
- Ventura T, Manor R, Aflalo ED, Chalifa-Caspi V, Weil S, Sharabi O, Sagi A. Post-embryonic transcriptomes of the prawn Macrobrachium rosenbergii: multigenic succession through metamorphosis. PLoS One. 2013;8(1):e55322.View ArticlePubMedPubMed CentralGoogle Scholar
- Saito T, Pšenička M, Goto R, Adachi S, Inoue K, Arai K, Yamaha E. The origin and migration of primordial germ cells in sturgeons. PLoS One. 2014;9(2):e86861.View ArticlePubMedPubMed CentralGoogle Scholar
- Duong TY, Scribner KT, Forsythe PS, Crossman JA, Baker EA. Interannual variation in effective number of breeders and estimation of effective population size in long-lived iteroparous lake sturgeon (Acipenser fulvescens). Mol Ecol. 2013;22(5):1282–94.View ArticlePubMedGoogle Scholar
- Narum SR, Buerkle CA, Davey JW, Miller MR, Hohenlohe PA. Genotyping-by-sequencing in ecological and conservation genomics. Mol Ecol. 2013;22(11):2841–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Raymakers C. CITES, the Convention on International Trade in Endangered Species of Wild Fauna and Flora: its role in the conservation of Acipenseriformes. J Appl Ichthyol. 2006;22(s1):53–65.View ArticleGoogle Scholar
- Vidotto M, Grapputo A, Boscari E, Barbisan F, Coppe A, Grandi G, Kumar A, Congiu L. Transcriptome sequencing and de novo annotation of the critically endangered Adriatic sturgeon. BMC Genomics. 2013;14(1):407.View ArticlePubMedPubMed CentralGoogle Scholar
- Hagihara S, Yamashita R, Yamamoto S, Ishihara M, Abe T, Ijiri S, Adachi S. Identification of genes involved in gonadal sex differentiation and the dimorphic expression pattern in undifferentiated gonads of Russian sturgeon Acipenser gueldenstaedtii Brandt & Ratzeburg, 1833. J Appl Ichthyol. 2014;30(6):1557–64.View ArticleGoogle Scholar
- Yuan L, Zhang X, Li L, Jiang H, Chen J. High-throughput sequencing of MicroRNA transcriptome and expression assay in the sturgeon, Acipenser schrenckii. PLoS One. 2014;9(12):e115251.View ArticlePubMedPubMed CentralGoogle Scholar
- Song W, Jiang K, Zhang F, Lin Y, Ma L. Transcriptome sequencing, De Novo assembly and differential gene expression analysis of the early development of Acipenser baeri. PLoS One. 2015;10(9):e0137450.View ArticlePubMedPubMed CentralGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10.Google Scholar
- Ekblom R, Galindo J. Applications of next generation sequencing in molecular ecology of non-model organisms. Heredity. 2011;107(1):1–15.View ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512.View ArticlePubMedGoogle Scholar
- Strickler SR, Bombarely A, Mueller LA. Designing a transcriptome next-generation sequencing project for a nonmodel plant species1. Am J Bot. 2012;99(2):257–66.View ArticlePubMedGoogle Scholar
- Jin SB, Zhang Y, Dong XL, Xi QK, Song D, Fu HT, Sun DJ. Comparative transcriptome analysis of testes and ovaries for the discovery of novel genes from Amur sturgeon (Acipenser schrenckii). Genet Mol Res. 2015;14(4):18913–27.View ArticlePubMedGoogle Scholar
- Yue H, Li C, Du H, Zhang S, Wei Q. Sequencing and De Novo assembly of the gonadal transcriptome of the endangered Chinese sturgeon (Acipenser sinensis). PLoS One. 2015;10(6):e0127332.View ArticlePubMedPubMed CentralGoogle Scholar
- Vidotto M. Transcriptomic analysis of the polyploid Adriatic sturgeon (Acipenser naccarii). 2013.Google Scholar
- Grimanelli D, Perotti E, Ramirez J, Leblanc O. Timing of the maternal-to-zygotic transition during early seed development in maize. Plant Cell. 2005;17(4):1061–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Pelegri F. The role of maternal factors in early zebrafish development. Fish Development and Genetics, Singapore: World Scientific. 2004; 1-38.Google Scholar
- Moreton J, Dunham SP, Emes RD. A consensus approach to vertebrate de novo transcriptome assembly from RNA-seq data: assembly of the duck (Anas platyrhynchos) transcriptome. Front Genet. 2014;5.Google Scholar
- Langer E-M. Methods of experimental yolk removal from Brachydanio rerio eggs. Experientia. 1979;35(9):1269–70.View ArticleGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.View ArticlePubMedGoogle Scholar
- Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.View ArticlePubMedPubMed CentralGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3(7):1–12.View ArticleGoogle Scholar
- Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64(15):5245–50.View ArticlePubMedGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26(6):509–15.View ArticlePubMedGoogle Scholar