Transcriptome Bioinformatical Analysis of Vertebrate Stages of Schistosoma japonicum Reveals Alternative Splicing Events

Identification and expression profiles of olfactory-related genes in the antennal transcriptome of Graphosoma rubrolineatum (Hemiptera: Pentatomidae)

Graphosoma rubrolineatum (Hemiptera: Pentatomidae) is an important pest of vegetables and herbs (e.g., Umbelliferae and Cruciferae) in China, Siberia, Korea, and Japan. Insects are highly dependent on their olfactory system to detect odorants. However, no molecular-mediated olfactory genes in G. rubrolineatum have yet been identified. In this study, we first established the antennal transcriptome of G. rubrolineatum and identified 189 candidate olfactory genes, including 31 odorant-binding proteins (OBPs), 15 chemosensory proteins (CSPs), four sensory neuron membrane proteins (SNMPs),94 odorant receptors (ORs), 23 ionotropic receptors (IRs), and 22 gustatory receptors (GRs). Additionally, phylogenetic trees were constructed for olfactory genes between G. rubrolineatum and other hemipteran insects. We also detected the expression profiles of ten OBPs, five CSPs, two SNMPs, five ORs, four IRs, and four GRs by real-time quantitative PCR. The results revealed that most genes (GrubOBP1/11/31, GrubCSP3/8, GrubSNMP1a/1b, GrubOrco/OR9/11/13, GrubGR1/4/22, GrubIR25/75h/76b/GluR1) were highly expressed in the antennae, GrubOBP13/31 and GrubCSP4/11/12 were highly expressed in the legs, while GrubOBP20 and GrubGR19 were highly expressed in the wings. Our results will enrich the gene inventory of G. rubrolineatum and provide further insight into the molecular chemosensory mechanisms of G. rubrolineatum.

Genetic difference between two Schistosoma japonicum isolates with contrasting cercarial shedding patterns revealed by whole genome sequencing

Schistosoma japonicum is one of the major infectious agents of human schistosomiasis, mainly endemic in China and the Philippines. We have previously reported the finding of two schistosome isolates, each with a different cercarial emergence pattern adapted to their different hosts. However, there are currently no whole-genome sequencing studies to investigate the underlining genetics of the adaptive traits. We sampled schistosomes in 2013 and 2020 from a hilly area Shitai (ST) and a marshland area Hexian (HX) of Anhui, China. Ten to 15 male or female adult worms from each site/year were sent for whole genome sequencing. Genetics were analyzed, and selection signals along genomes were detected. Gene enrichment analysis was performed for the genome regions under selection. The results revealed considerable genetic differentiation between the two isolates. The genome "windows" affected by natural selection were fewer in ST (64 windows containing 78 genes) than in HX (318 windows containing 276 genes). Twelve significantly enriched genes were identified in ST, but none in HX. These genes were mainly related to specific DNA binding and intercellular signaling transduction. Some functional region changes identified along the genome of the hilly schistosome may be related to its unique late afternoon cercarial emergence.

Comparative transcriptome profiles of Schistosoma japonicum larval stages: Implications for parasite biology and host invasion

Schistosoma japonicum is prevalent in Asia with a wide mammalian host range, which leads to highly harmful zoonotic parasitic diseases. Most previous transcriptomic studies have been performed on this parasite, but mainly focus on stages inside the mammalian host. Moreover, few larval transcriptomic data are available in public databases. Here we mapped the detailed transcriptome profiles of four S. japonicum larval stages including eggs, miracidia, sporocysts and cercariae, providing a comprehensive development picture outside of the mammalian host. By analyzing the stage-specific/enriched genes, we identified functional genes associated with the biological characteristic at each stage: e.g. we observed enrichment of genes necessary for DNA replication only in sporocysts, while those involved in proteolysis were upregulated in sporocysts and/or cercariae. This data indicated that miracidia might use leishmanolysin and neprilysin to penetrate the snail, while elastase (SjCE2b) and leishmanolysin might contribute to the cercariae invasion. The expression profile of stem cell markers revealed potential germinal cell conversion during larval development. Additionally, our analysis indicated that tandem duplications had driven the expansion of the papain family in S. japonicum. Notably, all the duplicated cathepsin B-like proteases were highly expressed in cercariae. Utilizing our 3rd version of S. japonicum genome, we further characterized the alternative splicing profiles throughout these four stages. Taken together, the present study provides compressive gene expression profiles of S. japonicum larval stages and identifies a set of genes that might be involved in intermediate and definitive host invasion.

Genome-wide transcriptome analysis of the early developmental stages of Echinococcus granulosus protoscoleces reveals extensive alternative splicing events in the spliceosome pathway

BACKGROUND

The complex life cycle of Echinococcus granulosus involves numerous environmental conditions within different intermediate and definitive hosts. This requires adaptation at different levels of transcript regulation. Alternative splicing (AS) and the related cellular functions as one of the major fields of post-genomics has been poorly studied in tapeworms. In the present study, we investigated AS events and their potential biological effects in E. granulosus.

METHODS

Whole transcriptome sequencing data of four groups of protoscoleces were prepared for RNA-seq library construction. Fresh protoscoleces were either used as non-induced controls (NT group) or incubated for 15 min with pepsin (PEP group) and cultivated in a biphasic medium for 12 and 24 h (12 and 24 h groups). The frequency and different types of AS events were identified using rMATS software. Functional annotations and gene ontology of differential AS (DAS) genes were performed using Blast2GO software. AS events were experimentally validated by PCR on the protoscolex cDNAs using specific primers for each gene.

RESULTS

At least one AS event was found in 38.1% of the genes (3904 out of 10,245) in the protoscoleces during early strobilar development. The genes were associated primarily with cellular and metabolic processes and binding and catalytic activity. KEGG pathway analysis of DAS events revealed a number of genes belonging to different components of the spliceosome complex. These genes tended to belong to common SR proteins, U1-related factors, U2-related factors, complex A-specific factors and other splicing-related proteins.

CONCLUSIONS

The high number of AS events in the transcriptome regulatory mechanisms indicates the essential rapid molecular processes required by the parasite for adaptation in different environments.

Comparative analysis of transcriptional profiles of Schistosoma japonicum adult worms derived from primary-infected and re-infected water buffaloes

BACKGROUND

Schistosoma japonicum (S. japonicum) is an important zoonotic parasite that is prevalent in China and parts of Southeast Asia. Water buffaloes are an important reservoir and the main transmission sources of S. japonicum. However, self-curing and resistance to re-infection have been observed in water buffaloes.

RESULTS

In this study, we compared the morphometry and differences in transcriptional expression of adult S. japonicum worms recovered from primary-infected and re-infected water buffaloes using Illumina RNA-sequencing (RNA-Seq) technology. Results of morphometry analysis revealed that adult S. japonicum worms recovered from re-infected water buffaloes were runtish with smaller organs. The ventral length of male worms was shorter in re-infected buffaloes (328 ± 13 vs 273 ± 8 µm, P < 0.05), and in female worms the oral sucker length (44 ± 3 vs 33 ± 5 µm, P < 0.05), ovary length (578 ± 23 vs 297 ± 27 µm, P < 0.05) and width (150 ± 8 vs 104 ± 9 µm, P < 0.05) were shorter, with fewer eggs in the uteri (41 ± 2 vs 12 ± 1, P < 0.05). Of 13,605 identified genes, 112 were differentially expressed, including 51 upregulated and 61 downregulated genes, in worms from re-infected compared with primary-infected water buffaloes. Gene ontology (GO) enrichment analysis revealed that GO terms such as "oxidation-reduction process", "calcium-dependent phospholipid binding", "lipid binding" and "calcium ion binding" were significantly enriched in downregulated genes, whereas GO terms related to metabolism and biosynthesis were significantly enriched in upregulated genes. The results revealed that the downregulation of some important genes might contribute to a reduction in worm numbers and maldevelopment of surviving worms in re-infected water buffaloes. Furthermore, upregulation of genes related to metabolic processes and biosynthesis might be a compensatory mechanism of worms in disadvantageous environments.

CONCLUSIONS

To our knowledge, our results present the first large-scale transcriptional expression study identifying the differences between adult S. japonicum worms from primary-infected and re-infected water buffaloes, and particularly emphasize differential expression that may affect the survival and growth of worms in re-infected water buffalo. This will provide new insight into screening for anti-schistosome targets and vaccine candidates.

Transcriptomic analysis of male and female Schistosoma mekongi adult worms

Background

Schistosoma mekongi is one of five major causative agents of human schistosomiasis and is endemic to communities along the Mekong River in southern Lao People’s Democratic Republic (Laos) and northern Cambodia. Sporadic cases of schistosomiasis have been reported in travelers and immigrants who have visited endemic areas. Schistosoma mekongi biology and molecular biology is poorly understood, and few S. mekongi gene and transcript sequences are available in public databases.

Results

Transcriptome sequencing (RNA-Seq) of male and female S. mekongi adult worms (a total of three biological replicates for each sex) were analyzed and the results demonstrated that approximately 304.9 and 363.3 million high-quality clean reads with quality Q30 (> 90%) were obtained from male and female adult worms, respectively. A total of 119,604 contigs were assembled with an average length of 1273 nt and an N50 of 2017 nt. From the contigs, 20,798 annotated protein sequences and 48,256 annotated transcript sequences were obtained using BLASTP and BLASTX searches against the UniProt Trematoda database. A total of 4658 and 3509 transcripts were predominantly expressed in male and female worms, respectively. Male-biased transcripts were mostly involved in structural organization while female-biased transcripts were typically involved in cell differentiation and egg production. Interestingly, pathway enrichment analysis suggested that genes involved in the phosphatidylinositol signaling pathway may play important roles in the cellular processes and reproductive systems of S. mekongi worms.

Conclusions

We present comparative transcriptomic analyses of male and female S. mekongi adult worms, which provide a global view of the S. mekongi transcriptome as well as insights into differentially-expressed genes associated with each sex. This work provides valuable information and sequence resources for future studies of gene function and for ongoing whole genome sequencing efforts in S. mekongi.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3086-z) contains supplementary material, which is available to authorized users.

Comparative Study of Transcriptome Profiles of Mouse Livers and Skins Infected by Fork-Tailed or Non-Fork-Tailed Schistosoma japonicum

Schistosoma japonicum (S. japonicum) is a worldwide spread pathogen which penetrates host skin and then induces several diseases in infected host, such as fibrosis, formation of granulomas, hepatocirrhosis, and hepatomegaly. In present study, for the first time, transcriptomic profiles of mouse livers and skins infected by fork-tailed S. japonicum cercaria or non-fork-tailed S. japonicum cercaria were analyzed by using RNA-seq. The present findings demonstrated that transcriptomic landscapes of livers and skins infected by fork-tailed S. japonicum cercaria or non-fork-tailed S. japonicum cercaria were different. S. japonicum has great influence on hepatic metabolic processes. Fork-tailed S. japonicum cercaria upregulated hepatic metabolic processes, while non-fork-tailed S. japonicum cercaria downregulated hepatic metabolic processes. For the metabolism process or the metabolism enzyme expressional change, the pharmacokinetics of host could be changed during S. japonicum infection, regardless the biotypes of S. japonicum cercariae. The changes of infected skins focused on upregulation of immune response. During the S. japonicum skin infection period, fork-tailed S. japonicum cercaria infection induced stronger immune response comparing with that immune response triggered by non-fork-tailed S. japonicum cercaria. The transcription factor enrichment analysis showed that Irf7, Stat1 and Stat2 could play important roles in gene regulation during fork-tailed S. japonicum cercaria infection.

Integrated analysis of differential expression and alternative splicing of non-small cell lung cancer based on RNA sequencing

Non-small cell lung cancer (NSCLC) is the most common type of lung cancer, with high morbidity and mortality rates. Numerous diagnosis and treatment methods have been proposed, and the prognosis of NSCLC has improved to a certain extent. However, the mechanisms of NSCLC remain largely unknown, and additional studies are required. In the present study, the RNA sequencing dataset of NSCLC was downloaded from the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/). The clean reads obtained from the raw data were mapped to the University of California Santa Cruz human genome (hg19), based on TopHat, and were assembled into transcripts via Cufflink. The differential expression (DE) and differential alternative splicing (DAS) genes were screened out through Cuffdiff and rMATS, respectively. The significantly enriched gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes pathways were obtained through the Database of Annotation, Visualization and Integrated Discovery (DAVID). Different numbers of DE and DAS genes were identified in different types of NSCLC samples, but a number of common functions and pathways were obtained, including biological processes associated with abnormal immune and cell activity. GO terms and pathways associated with substance metabolism, including the insulin signaling pathway and oxidative phosphorylation, were enriched in DAS genes rather than DE genes. Integrated analysis of differential expression and alternative splicing may be helpful in understanding the mechanisms of NSCLC, in addition to its early diagnosis and treatment.

Genome-Wide Transcriptome Analysis Reveals Extensive Alternative Splicing Events in the Protoscoleces of Echinococcus granulosus and Echinococcus multilocularis

Alternative splicing (AS), as one of the most important topics in the post-genomic era, has been extensively studied in numerous organisms. However, little is known about the prevalence and characteristics of AS in Echinococcus species, which can cause significant health problems to humans and domestic animals. Based on high-throughput RNA-sequencing data, we performed a genome-wide survey of AS in two major pathogens of echinococcosis-Echinococcus granulosus and Echinococcus multilocularis. Our study revealed that the prevalence and characteristics of AS in protoscoleces of the two parasites were generally consistent with each other. A total of 6,826 AS events from 3,774 E. granulosus genes and 6,644 AS events from 3,611 E. multilocularis genes were identified in protoscolex transcriptomes, indicating that 33–36% of genes were subject to AS in the two parasites. Strikingly, intron retention instead of exon skipping was the predominant type of AS in Echinococcus species. Moreover, analysis of the Kyoto Encyclopedia of Genes and Genomes pathway indicated that genes that underwent AS events were significantly enriched in multiple pathways mainly related to metabolism (e.g., purine, fatty acid, galactose, and glycerolipid metabolism), signal transduction (e.g., Jak-STAT, VEGF, Notch, and GnRH signaling pathways), and genetic information processing (e.g., RNA transport and mRNA surveillance pathways). The landscape of AS obtained in this study will not only facilitate future investigations on transcriptome complexity and AS regulation during the life cycle of Echinococcus species, but also provide an invaluable resource for future functional and evolutionary studies of AS in platyhelminth parasites.

De novo transcriptomic analysis of the female and male adults of the blood fluke Schistosoma turkestanicum

Background

Schistosoma turkestanicum is a parasite of considerable veterinary importance as an agent of animal schistosomiasis in many countries, including China. The S. turkestanicum cercariae can also infect humans, causing cercarial dermatitis in many countries and regions of the world. In spite of its significance as a pathogen of animals and humans, there is little transcriptomic and genomic data in the public databases.

Methods

Herein, we performed the transcriptome Illumina RNA sequencing (RNA-seq) of adult males and females of S. turkestanicum and de novo transcriptome assembly.

Results

Approximately 81.1 (female) and 80.5 (male) million high-quality clean reads were obtained and then 29,526 (female) and 41,346 (male) unigenes were assembled. A total of 34,624 unigenes were produced from S. turkestanicum females and males, with an average length of 878 nucleotides (nt) and N50 of 1480 nt. Of these unigenes, 25,158 (72.7 %) were annotated by blast searches against the NCBI non-redundant protein database. Among these, 21,995 (63.5 %), 22,189 (64.1 %) and 13,754 (39.7 %) of the unigenes had significant similarity in the NCBI non-redundant protein (NR), non-redundant nucleotide (NT) and Swiss-Prot databases, respectively. In addition, 3150 unigenes were identified to be expressed specifically in females and 1014 unigenes were identified to be expressed specifically in males. Interestingly, several pathways associated with gonadal development and sex maintenance were found, including the Wnt signaling pathway (103; 2 %) and progesterone-mediated oocyte maturation (77; 1.5 %).

Conclusions

The present study characterized and compared the transcriptomes of adult female and male blood fluke, S. turkestanicum. These results will not only serve as valuable resources for future functional genomics studies to understand the molecular aspects of S. turkestanicum, but also will provide essential information for ongoing whole genome sequencing efforts on this pathogenic blood fluke.