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

2020 ◽  
Author(s):  
Xiujuan Zhang ◽  
Jiabin Zhou ◽  
Linmiao Li ◽  
Wenzhong Huang ◽  
Hafiz Ishfaq Ahmad ◽  
...  

Abstract 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 2,782 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. 3,028), 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 percent (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 4,339 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.

2019 ◽  
Author(s):  
Xiujuan Zhang ◽  
Jiabin Zhou ◽  
Linmiao Li ◽  
Wenzhong Huang ◽  
Hafiz Ishfaq Ahmad ◽  
...  

Abstract 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 2,782 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. 3,028), 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 percent (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 4,339 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.


2019 ◽  
Author(s):  
Xiujuan Zhang ◽  
Jiabin Zhou ◽  
Linmiao Li ◽  
Wenzhong Huang ◽  
Hafiz Ishfaq Ahmad ◽  
...  

Abstract 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 2,782 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. 3,028), and functional assignment indicated that 6 of 14 annotated KEGG pathways were directly ovary-related and had abundant transcripts and differential expression genes. Importantly, 60 early gametogenesis-related genes (involving 755 unigenes) were successfully identified, and exactly 50 percent (30/60) of those showed differential expression in testes and ovaries. The Amh and Gsdf with testis-biased expression, and 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 of 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 4,339 predicted transcript factors (TFs) were also respectively identified, which 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. These provide candidate genes and theoretical basis for further the mechanisms of reproduction regulation of sturgeon.


PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e9133 ◽  
Author(s):  
Hongjun Yang ◽  
Danping Xu ◽  
Zhihang Zhuo ◽  
Jiameng Hu ◽  
Baoqian Lu

Background Red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is one of the most destructive insects for palm trees in the world. However, its genome resources are still in the blank stage, which limits the study of molecular and growth development analysis. Methods In this study, we used PacBio Iso-Seq and Illumina RNA-seq to first generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th larva, female and male) to increase our understanding of the life cycle and molecular characteristics of R. ferrugineus. Results A total of 63,801 nonredundant full-length transcripts were generated with an average length of 2,964 bp from three developmental stages, including the 7th instar larva, pupa, female adult and male adult. These transcripts showed a high annotation rate in seven public databases, with 54,999 (86.20%) successfully annotated. Meanwhile, 2,184 alternative splicing (AS) events, 2,084 transcription factors (TFs), 66,230 simple sequence repeats (SSR) and 9,618 Long noncoding RNAs (lncRNAs) were identified. In summary, our results provide a new source of full-length transcriptional data and information for the further study of gene expression and genetics in R. ferrugineus.


Genes ◽  
2021 ◽  
Vol 12 (2) ◽  
pp. 233
Author(s):  
Daling Wang ◽  
Ying Li ◽  
Reyilamu Aierken ◽  
Qi Kang ◽  
Xianyan Wang ◽  
...  

Cetaceans are a group of secondary aquatic mammals whose ancestors returned to the ocean from land, and during evolution, their immune systems adapted to the aquatic environment. Their skin, as the primary barrier to environmental pathogens, supposedly evolved to adapt to a new living environment. However, the immune system in the skin of cetaceans and the associated molecular mechanisms are still largely unknown. To better understand the immune system, we extracted RNA from the sperm whale’s (Physeter macrocephalus) skin and performed PacBio full-length sequencing and RNA-seq sequencing. We obtained a total of 96,350 full-length transcripts with an average length of 1705 bp and detected 5150 genes that were associated with 21 immune-related pathways by gene annotation enrichment analysis. Moreover, we found 89 encoding genes corresponding to 33 proteins were annotated in the NOD-like receptor (NLR)-signaling pathway, including NOD1, NOD2, RIP2, and NF-κB genes, which were discussed in detail and predicted to play essential roles in the immune system of the sperm whale. Furthermore, NOD1 was highly conservative during evolution by the sequence comparison and phylogenetic tree. These results provide new information about the immune system in the skin of cetaceans, as well as the evolution of immune-related genes.


2021 ◽  
Author(s):  
Yu Bai ◽  
Yonglu Meng ◽  
Jianlin Luo ◽  
Hui Wang ◽  
Guoyong Li ◽  
...  

The Chinese giant salamander, Andrias davidianus, is the largest amphibian species in the world; it is thus an economically and ecologically important species. The skin of A. davidianus exhibits complex adaptive structural and functional adaptations to facilitate survival in aquatic and terrestrial ecosystems. Here, we report the first full-length amphibian transcriptome from the dorsal skin of A. davidianus, which was assembled using hybrid sequencing and the PacBio and Illumina platforms. A total of 153,038 transcripts were hybrid assembled (mean length of 2,039 bp and N50 of 2,172 bp), and 133,794 were annotated in at least one database (nr, Swiss-Prot, KEGG, KOGs, GO, and nt). A total of 58,732, 68,742, and 115,876 transcripts were classified into 24 KOG categories, 1,903 GO term categories, and 46 KEGG pathways (level 2), respectively. A total of 207,627 protein-coding regions, 785 transcription factors, 27,237 potential long non-coding RNAs, and 8,299 simple sequence repeats were also identified. The hybrid-assembled transcriptome recovered more full-length transcripts, had a higher N50 contig length, and a higher annotation rate of unique genes compared with that assembled in previous studies using next-generation sequencing. The high-quality full-length reference gene set generated in this study will help elucidate the genetic characteristics of A. davidianus skin and aid the identification of functional skin proteins.


Genes ◽  
2019 ◽  
Vol 10 (7) ◽  
pp. 507 ◽  
Author(s):  
Zhang

Taenia pisiformis is a tapeworm causing economic losses in the rabbit breeding industry worldwide. Due to the absence of genomic data, our knowledge on the developmental process of T. pisiformis is still inadequate. In this study, to better characterize differential and specific genes and pathways associated with the parasite developments, a comparative transcriptomic analysis of the larval stage (TpM) and the adult stage (TpA) of T. pisiformis was performed by Illumina RNA sequencing (RNA-seq) technology and de novo analysis. In total, 68,588 unigenes were assembled with an average length of 789 nucleotides (nt) and N50 of 1485 nt. Further, we identified 4093 differentially expressed genes (DEGs) in TpA versus TpM, of which 3186 DEGs were upregulated and 907 were downregulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes (KEGG) analyses revealed that most DEGs involved in metabolic processes and Wnt signaling pathway were much more active in the TpA stage. Quantitative real-time PCR (qPCR) validated that the expression levels of the selected 10 DEGs were consistent with those in RNA-seq, indicating that the transcriptomic data are reliable. The present study provides comparative transcriptomic data concerning two developmental stages of T. pisiformis, which will be of great value for future functional studies on the regulatory mechanisms behind adult worm pathogenesis and for developing drugs and vaccines against this important parasite.


Author(s):  
Jie Yang ◽  
Chi Zhang ◽  
Wei-Hong Li ◽  
Tian-Er Zhang ◽  
Guang-Zhong Fan ◽  
...  

Background:: In Traditional Chinese Medicine (TCM), the heads and tails of Angelica sinensis (Oliv.) Diels (AS) is used in treating different diseases due to their different pharmaceutical efficacies. The underline mechanisms, however, have not been fully explored. Objective:: Novel mechanisms responsible for the discrepant activities between AS heads and tails were explored by a combined strategy of transcriptomes and metabolomics. Method:: Six pairs of the heads and tails of AS roots were collected in Min County, China. Total RNA and metabolites, which were used for RNA-seq and untargeted metabolomics analysis, were respectively isolated from each AS sample (0.1 g) by Trizol and methanol reagent. Subsequently, differentially expressed genes (DEGs) and discrepant pharmaceutical metabolites were identified for comparing AS heads and tails. Key DEGs and metabolites were quantified by qRT-PCR and targeted metabolomics experiment. Results:: Comprehensive analysis of transcriptomes and metabolomics results suggested that five KEGG pathways with significant differences included 57 DEGs. Especially, fourteen DEGs and six key metabolites were relation to the metabolic regulation of Phenylpropanoid biosynthesis (PB) pathway. Results of qRT-PCR and targeted metabolomics indicated that higher levels of expression of crucial genes in PB pathway, such as PAL, CAD, COMT and peroxidase in the tail of AS were positively correlated with levels of ferulic acid-related metabolites. The average content of ferulic acid in tails (569.58162.39 nmol/g) was higher than those in the heads (168.73  67.30 nmol/g) (P˂0.01); Caffeic acid in tails (3.82  0.88 nmol/g) vs heads (1.37  0.41 nmol/g) (P˂0.01), and Cinnamic acid in tails (0.24  0.09 nmol/g) vs heads (0.14  0.02 nmol/g) (P˂0.05). Conclusion:: Our work demonstrated that overexpressed genes and accumulated metabolites derived from PB pathway might be responsible for the discrepant pharmaceutical efficacies between AS heads and tails.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Haifeng Yan ◽  
Huiwen Zhou ◽  
Hanmin Luo ◽  
Yegeng Fan ◽  
Zhongfeng Zhou ◽  
...  

Abstract Background Although extensive breeding efforts are ongoing in sugarcane (Saccharum officinarum L.), the average yield is far below the theoretical potential. Tillering is an important component of sugarcane yield, however, the molecular mechanism underlying tiller development is still elusive. The limited genomic data in sugarcane, particularly due to its complex and large genome, has hindered in-depth molecular studies. Results Herein, we generated full-length (FL) transcriptome from developing leaf and tiller bud samples based on PacBio Iso-Seq. In addition, we performed RNA-seq from tiller bud samples at three developmental stages (T0, T1 and T2) to uncover key genes and biological pathways involved in sugarcane tiller development. In total, 30,360 and 20,088 high-quality non-redundant isoforms were identified in leaf and tiller bud samples, respectively, representing 41,109 unique isoforms in sugarcane. Likewise, we identified 1063 and 1037 alternative splicing events identified in leaf and tiller bud samples, respectively. We predicted the presence of coding sequence for 40,343 isoforms, 98% of which was successfully annotated. Comparison with previous FL transcriptomes in sugarcane revealed 2963 unreported isoforms. In addition, we characterized 14,946 SSRs from 11,700 transcripts and 310 lncRNAs. By integrating RNA-seq with the FL transcriptome, 468 and 57 differentially expressed genes (DEG) were identified in T1vsT0 and T2vsT0, respectively. Strong up-regulation of several pyruvate phosphate dikinase and phosphoenolpyruvate carboxylase genes suggests enhanced carbon fixation and protein synthesis to facilitate tiller growth. Similarly, up-regulation of linoleate 9S-lipoxygenase and lipoxygenase genes in the linoleic acid metabolism pathway suggests high synthesis of key oxylipins involved in tiller growth and development. Conclusions Collectively, we have enriched the genomic data available in sugarcane and provided candidate genes for manipulating tiller formation and development, towards productivity enhancement in sugarcane.


Sign in / Sign up

Export Citation Format

Share Document