https://orcid.org/0000-0002-2057-4023
Figures
Abstract
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.
Author summary
Schistosomes are parasitic flatworms that require a snail host and a mammalian host to complete their life cycle. Due to the difficulties in obtaining materials, little is known about the molecular aspects of this fluke’s larval stages. Based on RNA-Seq, we provide the first high-resolution, transcriptomic analysis of four larval stages of Schistosoma japonicum. The data showed the biological and physical features of each stage, also highlighted that miracidia and cercariae might use a different group of proteases for host invasion. Additionally, it indicated that different populations of germinal cells may exist in the larval stages. The high expression of tandem duplicated cathepsin B-like proteases at the cercariae stage may contribute to the wide definitive host range of S. japonicum. Additionally, we observed that alternative splicing plays a vital role in regulating gene expression in S. japonicum, among which skip exon was the most predominant. Our data provide valuable information on the expression and function of S. japonicum genes across their larval stages and will support basic and applied research for the community.
Citation: Cheng S, Zhu B, Luo F, Lin X, Sun C, You Y, et al. (2022) Comparative transcriptome profiles of Schistosoma japonicum larval stages: Implications for parasite biology and host invasion. PLoS Negl Trop Dis 16(1): e0009889. https://doi.org/10.1371/journal.pntd.0009889
Editor: James Cotton, Wellcome Trust Sanger Institute, UNITED KINGDOM
Received: April 7, 2021; Accepted: October 8, 2021; Published: January 13, 2022
Copyright: © 2022 Cheng et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All raw sequencing data are available via NCBI under SRA accessions PRJNA719283. Nucleotide sequences have been deposited in the Sequence Read Archive (SRA) of NCBI under accession codes SRR14133806-SRR14133817. The genome assembly and gene set of SjV3 are available through Zenodo (https://doi.org/10.5281/zenodo.5795038).
Funding: The project was foundeded by the National Key research and Development Project of China (Award Number: 2018YFA0507300) and the National Natural Science Foundation of China (Award Number: 31572513) to Wei Hu. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Schistosomiasis is a neglected tropical disease caused by Schistosoma spp., which is prevalent in approximately 78 countries and affects more than 240 million people worldwide [1]. Similar to other species, Schistosoma japonicum has a complex life cycle occurring in two hosts: snail and mammal. The life cycle outside the mammalian hosts involves eggs excreted with feces, asexually reproducing larvae in the snail host and two intermediate free-swimming stages. The mature egg is excreted with host feces and releases a miracidium in fresh water. Then the ciliated and free-swimming larva seeks the snail-Oncomelania genus guided by chemical attraction and further penetrates the snail host assisted by secretions, which are probably proteolytic enzymes [2]. Within the snail, the miracidium loses its ciliated plates and undergoes a dramatic developmental conversion into a mother sporocyst that contains a population of totipotent stem cells, called germinal cells. The germinal cell will proliferate and form the germinal balls (primordial daughter sporocyst embryos) to further produce daughter sporocysts. Upon maturation, the daughter sporocysts emerge from the mother sporocyst and migrate to the snail’s hepatopancreas where they undergo development to produce the free-swimming cercariae by asexual reproduction [3]. Under suitable light and temperature conditions, the fork-tailed cercariae will shed from the snail. There are five pairs of acetabular glands inside the cercarial head which contain many kinds of proteases [4]. Once contact with the skin of human and mammals, the cercariae complete the penetration process by the mechanical activity of the tail and the hydrolytic activities of the proteolytic enzymes.
Transcriptomic studies across the life cycle of schistosomes have been extensively conducted, but mainly by ESTs, SAGE or microarray approaches [5–10]. These techniques lack the accuracy and sensitivity of a more contemporary RNA-seq approach, which is a powerful tool for delivering genome-wide transcription profiles unconstrained by genomic annotation. Previous studies applying RNA-seq mainly focused on adult stages, gonads [11–14] or parasites derived from different hosts [15,16]. Dynamic transcriptome profiles of S. japonicum and S. mansoni from juvenile schistosomulae to sex mature adult worms offered many insights into the reproductive development of the parasites during intra-mammalian development [17,18]. Recently, the intramolluscan transcriptomes of S. mansoni were investigated [19], but the samples used for RNA-Seq were infected whole snails. Thus, the landscape of gene expression in ex vivo larval stages, especially the molecular basis behind the stage transitions, is still lacking.
To obtain the gene expression information in the larval stages, we performed RNA-seq analysis on S. japonicum eggs, miracidia, sporocysts and cercariae. We used the most recent 3rd version of S. japonicum genome for mapping and annotation. Gene expression information correlated well with the biology of each life stage. Miracidia and cercariae showed high motor and proteolysis activity, ready for the host invasion. DNA replication and cell division only occurred in the sporocysts. We identified genes that are stage specifically expressed or with enriched expression that could thus be vital for the dominant functions of the parasite in those life stages. We found that each larval stage has germinal cells and there may be germinal cell conversion during the larval development. Interestingly, we discovered the tandem duplication events drove the expansion of the papain gene family in S. japonicum. Furthermore, we identified a large number of alternative splicing (AS) events in each stage, indicating that AS is a widespread process for generating protein isoform diversity in S. japonicum. This study provides rich and valuable resources for the community to understand the larval biology and will assist in the exploration of novel anti-schistosome targets and vaccine candidates.
Materials and methods
Ethics statement
All experiments involving animals were carried out in accordance with the guidelines for the Care and Use of Laboratory Animals of the Ministry of Science and Technology of the People’s Republic of China (2006398) and approved by the Ethics and Animal Welfare Committee of the National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention, Shanghai, China (IPD2008-4).
Parasite material
All parasite material was from an Anhui isolate of S. japonicum maintained in the National Institute of Parasitic Diseases, Chinese Center for Diseases Control and Prevention, Shanghai.
Isolation of eggs
Eggs were isolated by an improved enzymatic method. Three New Zealand white rabbits were infected with 800–1000 cercariae each. Six weeks later, the liver tissues were chopped with a scalpel blade and homogenized in 500 mL saline solution. The suspension was successively passed through 80 and 180 mesh metal sieves. After repeated centrifugation and removing the tissue debris, the pellet was resuspended in 10 mL saline solution containing 100 μg collagenase IV (Solarbio Life Sciences, Beijing), then incubated at 37°C for 30 min with gentle shaking. The sample was then centrifuged at 2,000 rpm at 4°C for 8 min, and the residues after digestion were removed. The egg pellet was then washed twice with saline solution. Finally, the eggs were washed by pipetting on a 300-mesh nylon screen, then collected and stored in 1.2% NaCl solution at 4°C under dark. Eggs isolated from one rabbit were used as a biological replicate.
Isolation of the miracidia
Purified eggs were transferred into a 200 mL hatching measuring cylinder wrapped completely in light-blocking black tape with the exclusion of the top 4 cm from the lip, thereby producing a light gradient. The hatching cylinder was topped with artificial pond water (0.46 μM FeCl3 ·6 H2O, 220 μM CaCl2 ·2 H2O, 100 μM MgSO4 ·7 H2O, phosphate buffer [313 μM KH2PO4, 14 μM (NH4)2SO4] pH 7.2) [4] until above the tape-covered area ~1.5cm and exposed to bright light at 28°C. Eggs were incubated for 2 h post-hatch, and the top 10 mL of miracidia-containing water (MCW) was collected for miracidia isolation. Hatched miracidia were collected by centrifugation at 8000 × g for 1 min at 4°C, and were then washed twice with water. Miracidia hatched from the eggs purified from the liver of one rabbit were used as one replicate.
Isolation of the daughter sporocysts
Daughter sporocysts were separated from the hepatopancreas of Oncomelania hupensis snails after ~5–6 weeks post-infection with dissecting needles under an optical microscope and purified after washing three times with sterilized PBS (pH 7.4). Daughter sporocysts collected from 6~10 infected snails were used as one replicate.
Collection of cercariae
To obtain S. japonicum cercariae, O. hupensis snails ~7–8 weeks post-infection were exposed to light in artificial pond water at 26°C for 2 h. The emerging cercariae were gravity-concentrated by cooling on ice for two hours, which prevented swimming, then concentrated by centrifugation. Cercariae collected from 50~80 infected snails were used as one replicate.
RNA isolation, library preparation and sequencing
Each larval stage has three biological replicates. For each replicate, we used 100~150 mg eggs, 15, 000~20, 000 miracidia, 1~2 mg daughter sporocysts and 15, 000~20, 000 cercariae. Parasites were homogenized in a 1 mL sterilized glass tissue grinder (Solarbio Life Sciences, Beijing) and total RNA was isolated using Qiagen RNeasy Micro Kit (Valencia, CA). RNA quality was assessed by 1% agarose gel electrophoresis and a NanoPhotometer spectrophotometer (Implen, Westlake Village, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA). RNA-Seq libraries were generated with the NEBNext Ultra Directional RNA Library Prep Kits (NEB, USA) according to the manufacturer’s protocol. After clusters generation on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia), the libraries were sequenced on an Illumina Novaseq platform (Novogene, Tianjin, China) with paired-end 150 bp.
Read mapping and data processing
Quality control (QC) of the raw sequencing data was performed using the FASTQC program. Low-quality reads and adapter sequences were trimmed using fastp tool v.0.20.1. (parameters: -q 15 -u 40 -n 5 -l 15) [20]. The clean reads were mapped to the chromosome-level S. japonicum reference genome (SjV3) using STAR 2.4.2a in twopassMode (parameters:—outFilterMultimapScoreRange 1—outFilterMultimapNmax 10—outFilterMismatchNmax 10—alignIntronMax 500000—sjdbScore 2—alignSJDBoverhangMin 3) [21] and further used to estimate the transcript abundance in TPM (Transcripts Per Kilobase million) using RSEM v.1.3.1 [22] with default parameters. These transcript abundances were imported into R and summarized with tximport v1.18.0 [23]. Principal Component Analysis (PCA) was performed using the prcomp function in the stats (v3.6.0) R package. Hierarchical clustering analysis (HCA) was performed with pheatmap (https://cran.r-project.org/web/packages/pheatmap/index.html). The R package DEseq2 v1.26.0 [24] was used to perform differential expression analysis. Gene Ontology (GO) enrichment analysis was performed with the R package clusterProfiler [25]. The P-values were corrected for multiple hypothesis testing with the Benjamini–Hochberg false-discovery rate procedure (adjusted P-value).
Gene clustering
To compare time-series gene expression data, the TPM values were clustered using R package Mfuzz v 2.46.0. Soft clustering was run with fuzzifier parameter set to m = 2.51, which was estimated using the mestimate() function. The cluster number was manually set at c = 8. Cluster members at a filter of 0.5 were used in subsequent enrichment analyses.
Identification and phylogenetic analysis of papain gene family
The protein and genome sequences of S. japonicum (SjV3) were downloaded from the Zenodo website (https://doi.org/10.5281/zenodo.5795038). The protein and genome sequences of S. mansoni were downloaded from the WormBase ParaSite (https://parasite.wormbase.org/index.html). The peptidase C1 (papain) family is part of clan CA of cysteine peptidases containing catalytic Cys25 (hereinafter papain numbering) and His159 residues in the active site [26]. To identify the papain genes in S. japonicum and S. mansoni, the HMMER profile of peptidase C1 domain (PF00112) was from the Pfam database and searched against the protein sequences of S. japonicum and S. mansoni with an E value cutoff of 1 × 10−40. To study the phylogenetic relationships of the papain genes between S. japonicum and S. mansoni, multiple sequence alignments of amino acid sequences were performed using the ClustalW program with default parameters. An unrooted neighbor-joining tree was constructed with 500 bootstrap replications using MEGA 7.0 software based on the full-length protein sequence alignment. All identified papain genes were classified into different groups based on gene annotation and the alignment results. Consensus Newick format trees were compiled with MEGA 7.0 software and edited with Adobe Illustrator.
Characterization of the gene structures, conserved motifs and chromosomal distributions of the papain genes
Gene annotations of the identified papain genes were extracted from the genome reference GFF files. The conserved domains were analyzed using MEME 5.1.1 software and the maximum number of motifs was set to 10. The gene structures and motif patterns were drawn using TBtools software [27]. Collinearity analysis of papain protein sequences from S. japonicum and S. mansoni were performed with BLASTP and MCScanX software [28]. Duplication of papain genes was analyzed using MCScanX with E-value < 1 × 10−5. The R package "pheatmap" was used to draw the heatmap of papain genes based on their expression levels across the four S. japonicum life stages.
Alternative splicing analysis
Alternative splicing (AS) events in these four S. japonicum life stages were detected and quantified using rMATs-turbo [29] with the reference transcript annotation. Percent spliced in (PSI) values were calculated for five classes of alternative splicing events, including skip exons (SE), retained introns (RI), mutually exclusive exons (MXE), alternative 5’ splice sites (A5SS), and alternative 3’ splice sites (A3SS). In total, AS events with a false discovery rate (FDR) < 0.05 between two adjacent stages were identified as differential AS events. Sashimi plots used to show splicing events were generated by rmats2sashimiplot (https://github.com/Xinglab/rmats2sashimiplot).
Results and discussion
An overview of the transcriptomes
Twelve digital expression (DGE) libraries were constructed and sequenced using the total RNA isolated from four S. japonicum stages (eggs, miracidia, sporocysts and cercariae) with three biological replicates (Fig 1A). Using RNA-Seq, we obtained 40,499,690 to 57,349,268 raw reads and 40,383,584 to 57,153,380 clean reads for each library. 85.84% to 92.31% of the clean reads were mapped to the 3rd version of the S. japonicum reference genome (S1 Table). To study the relationship between samples and evaluate the reproducibility of the biological replicates, we performed Principal component analysis (PCA), Pearson’s rank correlation analysis, and hierarchical clustering analysis (HCA). The PCA plot shows the top two principal components that explain most of the variance between samples in the data set, 70% and 17% for PC1 and PC2. Replicates from eggs and miracidia were clustered closely, and the sporocysts and cercariae were clustered away from other stages (Fig 1B). Pearson’s rank correlation analysis confirmed high reproducibility and consistent quality among the biological replicates (S1 Fig). Based on HCA, the expression profiles of the four life stages are clearly separated (S2 Fig). Additionally, these analyses revealed a high similarity between eggs and miracidia, which is reasonable because some eggs contain well-developed miracidia. To obtain the co-expressed and differentially expressed genes between the four stages, an inter-sample Venn diagram analysis was performed. The mRNA with an expression value equal to or greater than 1 TPM in at least two of the three biological replicates was considered expressed. A total of 8,732 genes were identified, and 7,925, 7,681, 8,223, and 7,801 genes were found in eggs, miracidia, sporocysts and cercariae, respectively. The number of genes identified in each larval stage was much higher than that reported by Gobert et al. [10] and Cai et al. [9]. Most of the genes (7,072) were identified in all four life stages, but 94, 21, 251, and 141 genes were exclusively detected in eggs, miracidia, sporocysts and cercariae (Fig 1C and S2 Table).
RNA-Seq was performed on replicate samples of eggs (Egg), miracidia (Mir), sporocysts (Spo) and cercariae (Cer). (A) Four S. japonicum larval stages used for RNA-Seq. (B) PCA results. Each symbol indicates an individual sample. (C) Venn diagram showing differentially transcribed genes among the four life stages. Circles of different colors represent a set of genes transcribed in one stage, where the values represent the number of uniquely transcribed genes in one stage, or common transcribed genes between two, three or four life stages.
Stage-specific genes (SSG) and stage-enriched genes (SEG)
Since the parasite at each stage shows distinct biological characteristics, we attempted to identify genes expressed specifically at each stage, or significantly more highly expressed at each stage. We thus defined SSG (genes expressed at only one stage) and SEG (genes have significantly higher expression at one stage compared to the other three with a FC > 5 and FDR < 0.05) to further analyze the genes related to their stage-related features. We obtained 94 SSG and 245 SEG at the egg stage, 21 SSG and 2 SEG in miracidia, 251 SSG and 255 SEG in sporocyst as well as 141 SSG and 209 SEG in cercaria (Fig 2A and S3 Table). By performing GO analysis, we explored the enriched functions of SSG and SEG at each life stage (Fig 2B and S4 Table).
(A) The number of SSG and SEG in each stage. (B) GO enrichment analysis based on the SSG and SEG in each stage. Gene ratio is the percentage of total SEG or SSG in the given GO term.
SSG and SEG in eggs
SEA (soluble egg antigen) secreted by egg induces granulomas in the mammalian host. At this stage, four members of T2 ribonucleases (Sjc_0002258, Sjc_0002263, Sjc_0002281, and Sjc_0002296) associated with ribonucleases T2 activity (GO: 0033897) showed specific expression. T2 ribonuclease was identified as one of the top 25 most highly enriched genes in eggs in a previous study [9]. T2 ribonucleases are found in the genome of protozoans, plants, animals and viruses. A broad range of biological roles for these ribonucleases have been revealed, including scavenging of nucleic acids, degradation of self-RNA, serving as extra- or intracellular cytotoxins, and modulating host immune responses [30]. Previous study reported that S. mansoni omega-1, one major component of the SEA that was specifically expressed in eggs [31], was able to drive Th2 cell responses [32]. Omega-1 is a member of T2 ribonuclease family, both its RNase activity and glycosylation are essential for Th2 skewing [31]. Omega-1 can be taken up by dendritic cells and alter their cytoskeletal structure and function [33,34]. The protein structure analysis revealed that S. mansoni Omega-1 and the four S. japonicum T2 ribonucleases all contained ribonuclease_T2 domain (S3 Fig). Therefore, we propose that these T2 ribonucleases specifically expressed in eggs may participate in the regulation of host immune response.
SSG and SEG in miracidia
Miracidia are free-swimming and penetrate snail host. We noticed that ten members of cilia- and flagella- associated protein (CFAP) (Sjc_0004811, Sjc_0003323, Sjc_0007156, Sjc_0000523, Sjc_0002745, Sjc_0000347, Sjc_0007012, Sjc_0000734, Sjc_0006213, and Sjc_0000078) related to motile cilium (GO: 0031514) were enriched at this stage. CFAP plays a role in the reproduction of the size and morphology of cilia [35]. Therefore, it may be involved in the assembly of the cilium during the development of miracidia. Alternatively, CFAP was reported to regulate olfactory transduction in mice [36]. Since miracidia are attracted by the chemical substances secreted by the snail [37] to locate the host, CFAP may participate in this process. In addition, a 5-hydroxytryptamine (5-HT, serotonin) receptor (Sjc_0004493) was only detected at this stage. The neurotransmitter molecule 5-HT regulates diverse physiological processes in both invertebrates and vertebrates [38]. In schistosomes, 5-HT treatment significantly stimulated motility of the in vitro cultured sporocysts or adult worms [39,40]. And its receptor Sm5HTR was required for the proper control of motility in S. mansoni [41]. Thus, we reasoned that this 5-HT receptor may be responsible for the control of movement in miracidia.
SSG and SEG in sporocysts
Sporocysts are resident in snail hosts and produce cercariae. One cercarial elastase (SjCE2b, Sjc_0008947) and nine leishmanolysins (Sjc_0006363, Sjc_0006650, Sjc_0006649, Sjc_0006646, Sjc_0006648, Sjc_0006647, Sjc_0006120, Sjc_0006219, and Sjc_0006399) associated with proteolysis (GO: 0006508) were specifically expressed at this stage. There was an expanded family of elastases in S. mansoni [42], but only one member in S. japonicum [43]. Cercarial elastases are the major invasive proteases in S. mansoni [44], and are considered to be involved in S. japonicum cercariae invasion as well [45]. Leishmanolysin (also called GP63) is a critical virulence factor in various Leishmania species. This important metalloprotease manipulates the host immune system by allowing the parasite to establish, survive and propagate within mammalian macrophages [46]. Recently, it was reported that a leishmanolysin derived from S. mansoni excretory/secretory products could interfere snail haemocyte morphology and migration [47]. It’s worth noting that the cercaria head is mostly transcriptionally and translationally quiescent [48]. The elastase and leishmanolysins transcribed in the daughter sporocysts may be prepared for the cercariae penetration. Additionally, genes necessary for DNA replication (GO: 0006160) and DNA replication initiation (GO: 0006270) were enriched in daughter sporocysts, such as DNA replication licensing factor MCM2 (MCM2, Sjc_0000705), G2/mitotic-specific cyclin-B2 (ccnb2, Sjc_0005264) and mitotic spindle assembly checkpoint protein MAD2A (MAD2L1, Sjc_0002379). These genes play pivotal roles in the mitosis process, which is consistent with active cell division to produce cercariae at this stage. We also noted that more hypothetical proteins were enriched in daughter sporocyst than those from other stages, which may contribute to the biological functions of this poorly studied stage.
SSG and SEG in cercariae
Cercariae are free-swimming and invade the mammalian host. One rhodopsin-type GPCR (Sjc_0002522) and two octopamine receptors (Sjc_0007275, Sjc_0000926) were specifically transcribed at this stage. Rhodopsin-type GPCRs were known to be involved in photoreception, typically thought of as light sensors in animals [49]. Given that the release of cercariae from snails is triggered by the sunlight, we speculated that it may participate in cercariae photokinesis. Octopamine (OA) is one of the invertebrate-specific biogenic amines. In locusts and mollusks, OA is involved in motor control [50]. In S. mansoni adult worms, OA is widely distributed in neurons of the peripheral nerve net that innervate muscle [51]. The octopamine receptors that specifically expressed in cercariae may play a role in movement control. Besides, calcium binding protein (CaBP, Sjc_0006317) related to calcium ion binding (GO: 0005509) exhibited the most enriched expression in the cercariae stage. The divalent cation calcium is used as a cellular signal or ionic cofactor involved in diverse metabolic processes, including secretion, metabolism, muscle movement and neuronal function [52]. In Paragonimus ohirai, the excystment of metacercaria is a calcium-dependent process [53]. In Trichobilharzia regenti, CaBP represented the fifth most differentially transcribed gene between cercariae and schistosomulae [54]. In S. mansoni, the preacetabular glands of cercariae contain a high concentration of calcium [55], and there is a downregulation of CaBP in cercariae following epidermal penetration [56]. These findings suggest that the function of calcium may be conserved among trematode species by regulating the larval physiology, while CaBP may be important for the infectious cercariae.
Differentially expressed genes (DEG) between adjacent life stages
During development, S. japonicum undergoes dramatic morphological changes as well as physiological changes. To view the transcriptional changes during the transition of these four stages, pairwise differential gene expression analysis was performed for adjacent life stages. A list of differentially expressed genes (DEG) with a P-value < 0.05 was generated for each of the different pairwise comparisons of the life stages. A total of 1,184 DEGs were associated with the development from egg to miracidium, of which 190 and 994 were up-and downregulated in miracidium, respectively. 3,530 DEG between miracidium and sporocyst while 3,066 DEG between sporocyst and cercaria were detected (Fig 3 and S5 Table). Besides, Gene Ontology (GO) enrichment analysis was performed based on the DEG to analyze the enriched gene functions, the fold change cutoff value was set at > 5 to obtain more representative GO terms (S4 Fig and S6 Table).
Volcano plots showing differentially expressed genes (DEGs) in miracidium compared to egg (A), in sporocyst compared to miracidium (B), and in cercaria compared to sporocyst (C). Genes with p-value near 0 were adjusted to the p-value = 10−500. The dot line on the y axis refers to p-value = 0.05. The dotted lines on the x axis refer to fold change = -2 or 2.
During the development from egg to miracidium, two GO categories were upregulated, including magnesium ion binding (GO: 000287). The associated gene was a 7-methylguanosine phosphate-specific 5’-nucleotidase (Sjc_0000850). 5′-nucleotidases are enzymes catalyzing the hydrolytic dephosphorylation of nucleoside monophosphates to nucleosides and orthophosphate [57]. As catabolic enzymes, they play important roles in the regulation of nucleotide levels in living cells [58]. Compared to miracidia, fifteen GO categories were upregulated in eggs, including those related to lipid metabolic process (GO: 0006629), pyridoxal phosphate binding (GO: 0030170) and heme binding (GO: 0020037). Schistosomes are reported to uptake the basic sterols and fatty acids from host [59]. The pyridoxal phosphate (the active form of vitamin B6) [60] as well as heme, the important cofactor for oxygen transfer [61] were also derived from the host. Since the egg is produced and developed in the host, the high lipid metabolism activity and binding activity of vitamin B6 and heme may contribute to nutrients and oxygen that are required for the embryonic development in the egg.
During the development from miracidium to sporocyst, ten GO categories were enriched from the upregulated genes in miracidia. The enrichment of protein glycosylation (GO: 0006486) and fructosyltransferase activity (GO: 0008417) may be associated with the enhanced formation of glycocalyx that covers the surface of the miracidium [62]. Besides, a group of leishmanolysin and neprilysin linked to metalloendopeptidase activity (GO: 0004222) were upregulated. Notably, these leishmanolysins were different from those specifically expressed in sporocysts. Neprilysin (NEP) is a zinc-metalloenzyme belonging to the M13 family [63], which plays an important role in the interactions between host and parasite. Neplilysin was shown to involve in the production of immunoactive peptides, which could inactive the immunocytes from the snail host [64]. These data indicate that leishmanolysin and neprilysin may help the parasite escaping the immune attack of the Oncomelania hupensis hemocytes after miracidia penetration. On the opposite, ten GO categories were upregulated in sporocysts, including DNA replication (GO: 0006260), DNA replication initiation (GO: 0006270) and DNA binding (GO: 0003677), which are consistent with the asexual reproduction of sporocysts that undergoes active cell division [19].
Ten GO categories were upregulated during the development from sporocyst to cercaria. Go terms related to glycolytic process (GO: 0006096), transmembrane transport (GO: 0055085) and microtubule-based process (GO: 0007017) represent the high motility and metabolic activity of cercaria. Besides, proteolysis (GO: 0006508) related cysteine-type (GO: 0004197), aspartic-type (GO: 0004190) and serine-type endopeptidase activity (GO: 0004252) were also enriched. Many members of these important protease family have been studied, such as cathepsin B and cathepsin L from the cysteine protease family [65], cathepsin D from the aspartic protease family [66] and trypsin and elastase from the serine protease family [67]. These proteases participate in numerous biological processes, such as the parasite’s invasion, survival and longevity in the definitive host.
Gene clusters based on expression pattern across the four life stages
To further explore the transcriptomic changes across the four life stages, genes were clustered into 8 groups based on their expression profile (Fig 4). Each cluster showed a specific expression pattern. Genes from clusters 1, 2 and 5 showed highest expression in eggs, miracidia and cercariae respectively, while those from clusters 3 and 4 were highly expressed in sporocysts. GO categories enriched in these clusters were generally consistent with those based on SSG and SEG of each stage (S7 Table). We noticed that cluster 5 contained 17 genes that are linked to G protein-coupled receptor (GPCR) signaling pathway (GO: 0007186), including one metabotropic glutamate receptor (Sjc_0001309), one tachykinin-like peptides receptor (Sjc_0003380), one rhodopsin GPCR (Sjc_0003938), two 5-HT receptors (Sjc_0001525, Sjc_0007274) and three octopamine receptors (Sjc_0005446, Sjc_0007275, and Sjc_0006297). Glutamate is a neurotransmitter that is involved in controlling parasite locomotion [68]. Tachykinins are peptides that play a role in the processing of sensory information and control of motor activities [69]. These GPCRs may coordinate to regulate cercariae movement, light response, host sensory and adaption to a new osmotic pressure environment after host penetration. Genes in cluster 6, 7 and 8 showed high expression at two stages. Among them, cluster 6 showed increased expression in both the miracidia and cercariae. The enriched GO categories were mainly related to protein synthesis and transport, such as endoplasmic reticulum (GO: 0005783), intracellular protein transport (GO:0006886) and vesicular mediated transport (GO: 0016192). This may reflect that parasites at these two free-swimming stages are preparing proteins used for the host invasion.
Membership value was calculated using R package Mfuzz, which indicated the degree of membership of this gene for the corresponding cluster. The y-axes were scaled independently to emphasize the differences between clusters.
Genes associated with germinal cell and asexual reproduction in the larval stages
In planarians (close relatives of schistosomes), a population of pluripotent stem cells called neoblasts can regenerate injured tissues and replenish a whole animal from a single cell [70,71]. In recent years, it was shown that schistosomes, like the planarians, also contain stem cells. There are two major stem cells that play in the different life stages during intramolluscan and intramammalian development, including the germinal cells of sporocysts involved in asexual reproduction and the neoblast of adult worms involved in sexual reproduction [72]. Germinal cells have similar morphology with the neoblasts of planarians, they have a high nuclear to cytoplasmic ratio, an open chromatin structure and a large nucleolus [73]. These cells seem to proliferate indefinitely, evidenced by the serial transplantation of sporocysts into naïve molluscan hosts that led to continuous propagation of the parasites [74].
Recently, many methods were developed to study schistosome stem cells, including 5-ethynyl-2′-deoxyuridine (EdU) labeling, whole-mount in situ hybridization (WISH), and RNA interference (RNAi) [75]. Wang et al. compared gene expression profiles of miracidia and sporocysts 48 h post-transformation with the transcripts enriched in FACS-purified planarian neoblasts [76]. The same group also transcriptionally profiled stem cells from Schistosoma mansoni in vitro transformed mother sporocysts at single-cell levels. Three major stem cells classes were identified upon their respective markers: k-cells (that transcribe klf and nanos-2); ϕ-cells (that transcribe fgfrA and fgfrB); and δ-cells that produce transcripts of nanos-2 and fgfrA, B [77]. ScRNAseq analysis was also performed in the 2-day old schistosomulae to better understand the cell types and tissue differentiation [78]. These works provided many valuable markers of schistosome stem cells. We picked 18 markers and identified the homologs in S. japonicum using reciprocal BLAST comparisons (E-value < e−10) (S8 Table). Then we described the expression profiles of these marker genes in the four larval stages (Fig 5A). Nanos-2 expression was observed in the eggs and miracidia, but declined in the sporocysts. The RNA-binding protein nanos-2 has been showed to function in schistosomes as a conserved regulator of adult stem cells [79]. It has been proposed that there are two germinal cell subpopulations (nanos-2+ and nanos-2- cells) in the mother sporocysts, with the later subpopulation proliferating more rapidly [76]. The two fibroblast growth factor receptors (fgfr)—fgfrA and fgfrB, showed different expression patterns across these stages. FgfrA is essential for maintenance for adult stem cells [79], it’s expression in miracidia, sporocysts and cercariae was stable (Fig 5B). However, fgfrB showed relatively high expression in the sporocysts (Fig 5B), suggestive of its important role in germinal cells. A polo-like kinase (plk) gene highly expressed in sporocysts, polo-like kinases are important regulators of cell cycle progression and mitosis [80]. Cell cycle related transcripts (h2a, cyclinB, and PCNA) all showed relatively high expression in sporocysts (Fig 5B). Ago2-1, a stem cell marker and enriched in germinal cells [76], was also upregulated in sporocysts (Fig 5B). ZFP-1 is a zinc finger protein which is essential for the specification of new tegument cells [81]. This gene was highly expressed at the cercariae stage (Fig 5B). Soria et al. discovered a novel stem/germ cell marker-calmodulin (cam) in the schistosomulae [78]. Calmodulin is a Ca2+ transporter and is critical for miracidium-to-sporocyst transformation [82], for sporocyst growth and for egg hatching [83]. We found this gene was upregulated in the cercariae stage (Fig 5B). In conclusion, these stem cell marker genes displayed differential expression patterns across the four larval stages. This indicated that germinal cells existed in every larval stage and stem cell conversion may have occurred during the larval development.
(A) Heatmaps showing the relative gene expression of the eighteen stem cell markers in egg, miracidium, sporocyst and cercaria, as indicated. (B) Expression profiles of the nine stem cell markers in the four larval stages.
Highly expressed cathepsin-B like proteases in cercariae were produced by tandem duplications
Peptidase C1 genes are cysteine proteases of the papain superfamily. These enzymes play key roles in the pathogenesis of both protozoan and metazoan parasites, including skin or tissue penetration, hydrolysis of host or parasite proteins, and evasion or modulation of the host immune response [65,84]. The multiple functions of cysteine proteases make them attractive chemotherapeutic and vaccine targets [85].
Through a Hidden Markov Model (KMM) search, we identified 14 and 30 papain genes in S. japonicum and S. mansoni (S1 Dataset), respectively. To obtain the phylogenetic relationship of Schistosoma papain proteins, an unrooted tree was constructed. Papains were classified into six groups, including cathepsin B, cathepsin B-like, cathepsin C, cathepsin L, cathepsin L-like and cathepsin S (Fig 6A). A schematic representing the structure of all papain proteins was constructed from the MEME motif analysis results. As exhibited in Fig 6B, papain members within the same groups were usually found to share a similar motif composition. The exon-intron organizations of all the identified papain genes were examined to gain more insight into the evolution of the papain family. As showed in Fig 6B, all papain genes contain two to eleven exons (1with two exons, 5 with three exons, 27 with four exons, 1 with five exons, 5 with seven exons, 2 with eight exons, 2 with nine exons, and 1 with eleven exons). The papain genes in the same group shared similar gene structures.
(A) Phylogenetic relationship of papain. An un-rooted phylogenic tree was constructed in MEGA 7 using multiple alignments of S. japonicum and S. mansoni whole sequences data. (B) The motif composition and exon-intron structure of papain genes. For the motif pattern, the motifs, numbers 1–10, are displayed in different colored boxes. For the gene structure, green boxes indicate untranslated 5’- and 3’-regions; yellow boxes indicate exons; black lines indicate introns. (C) Synteny analysis of papain genes between S. japonicum and S. mansoni. Gray lines in the background show the collinear blocks within S. japonicum and S. mansoni genomes, while the red lines highlight the syntenic papain gene pairs. The tandem duplicated genes are marked in red. (D) Expression profiles of the papain genes at the four life stages.
To explore the relationship between S. japonicum and S. mansoni papain genes, we determined their chromosomal locations and whether they originated from gene duplication events. We identified papain genes on five chromes, most of them are located on chromosome 8 (for S. japonicum) or chromosome 1 (for S. mansoni). Tandem-duplicated genes are defined as two paralogous genes that are separated by fewer than 10 intervening genes [86], and it is one of the main sources of diversity for the evolution of gene families in eukaryotic organisms [87]. We identified 17 S. japonicum papain genes (10 pairs) that correspond to the tandem duplication events. In S. mansoni, we only detected 2 tandem-duplicated papain genes, comprising 1 gene pair (Fig 6C). We also added papain genes of Schistosoma haematobium and Schistosoma bovis into phylogenetic analysis. There are 14 papain genes in S. haematobium and 15 in S. bovis. The unrooted tree showed that a branch of S. japonicum cathepsin B-like proteases (14 members) separated from others, ten of them were produced by tandem duplication (S5 Fig). These results suggested that tandem duplication played important roles in the expansion of the papain gene family in S. japonicum. We then evaluated the papain gene expression profiles in different S. japonicum life stages via the RNA-seq data, and found that 87.5% of papain genes, including all the tandem-duplicated cathepsin B-like proteases, were highly expressed in the cercaria stage (Fig 6D).
Here, for the first time, we characterized the phylogenetic relationships, gene and protein structures and chromosome locations of the papain gene family in S. japonicum. We found that tandem duplication events drove the expansion of the papain gene family in S. japonicum. Besides, all the duplicated cathepsin B-like proteases showed the highest expression at the cercariae stage. S. mansoni and S. haematobium primarily infect humans, but S. japonicum is a zoonotic specie that infects more than 40 different mammalian species [88]. Thus, we speculate that these tandem duplicated and cercariae highly-expressed cathepsins may play important roles in assisting S. japonicum in establishing infections in broad definitive hosts.
Genome-wide identifications and comparative analysis of alternative splicing events within the four life stages
Alternative splicing (AS) is the process that enables one gene to encode two or more mature mRNAs through the differential utilization of splice sites [89]. AS greatly expands the transcriptome and proteome diversity; it is widespread in the genomes of humans and other species [90,91]. In parasites, AS may play fundamental roles in the host-parasite interactions by producing alternative isoforms with different functions [92]. The discovery of AS giving rise to different isoforms of antigenic proteins is indicative of immune evasion strategies by the parasites [93]. AS has been investigated in S. japonicum schistosomulae and adult worms [94,95], but they were based on the early RNA-Seq technologies, the sequencing depth and accuracy were far inferior to the one used in our study. Here, based on the most accurate SjV3 genome, we characterized the detailed AS in the four larval stages of S. japonicum.
In this study, five major AS types were considered as described before [96], including skip exons (SE), retained introns (RI), mutually exclusive exons (MXE), alternative 5’ splice sites (A5SS) and alternative 3’ splice sites (A3SS). 6,099, 6,101, 5,949 and 6,282 AS genes were determined in eggs, miracidia, sporocysts and cercariae (S2 Dataset), accounting for 62%, 62%, 61%, and 64% expressed genes, respectively (Fig 7A). The distributions of AS types in eggs, miracidia, sporocysts and cercariae are comparable; that is, SE was the most abundant AS events (44%-49%), followed by RI (21%-27%), A5SS (10%-12%), A3SS (10%-11%) and MXE (6%-9%) (Fig 7B).
(A) Number and percentage of AS events in S. japonicum eggs, miracidia, sporocysts and cercariae. (B) Proportions of different AS types detected in S. japonicum eggs, miracidia, sporocysts and cercariae. (C) The Venn diagram shows the overlap of AS genes in the four life stages. (D) Interactions between the five types of detected AS genes in the cercaria stage were visualized using an UpSet plot. (E) Venn diagram of the overlap of the DE and DAS genes between the eggs and miracidia stages. (F) An example of a gene (Sjc_0007531) displayed different SE at the eggs and miracidia stages.
The overall statistics of shared/unique AS genes in the four life stages are shown in Fig 7C. Interestingly, 5,374 genes were undergone alternative splicing in all four stages, and there are 143, 117, 124, 241 genes uniquely spliced at the stage of eggs, miracidia, sporocysts and cercariae, respectively. Notably, different AS patterns may occur for a single gene. Thus, UpSet plots were used to depict the intersections between AS types for each stage (Figs 7D and S6). For example, at the cercariae stage, a total of 2,240 genes only occurred in one type of AS event. The number of the SE is the most, with 1,562 genes, and the least is MXE, with only 13 genes. Furthermore, 519 genes contained up to five types of AS events (Fig 7D).
To explore whether gene expression and AS acted cooperatively or independently to regulate S. japonicum development, both gene expression levels and gene with AS in the three stage transitions were simultaneously analyzed (S2 Dataset). Among the genes identified to be differentially expressed, only a small number overlapped with the genes had differential AS in the three comparisons: 42 (8.6%) between eggs and miracidia, 142 (11.5%) between miracidia and sporocysts, 106 (9.4%) between sporocysts and cercariae (Figs 7E and S7). Sjc_0007531, a cell cycle control protein, showed different SE in the eggs and miracidia (Fig 7F).
In conclusion, more than 60% of the expressed genes showed alternative splicing in each stage, much higher than previously detected in schistosomulae, which was 42.14%. The most common type of AS events detected in the larval stages was the SE, while the MXE was the least. In the schistosomulae and adult worms, the most common were SE and A3SS, the least were RI and MXE, respectively [94,95]. Among all events, RI is the predominant mode of AS in plants [97], whereas ES is the major type in humans [98]. However, in the cestodes Echinococcus granulosus and Echinococcus multilocularis or the free-living flatworm Schmidtea mediterranea, the major type AS was RI [99,100]. Hence, our results suggest that the gene regulation pattern of S. japonicum is much closer to its human hosts. Besides, we found that, similar to Trypanosoma cruzi and Arabidopsis [101,102], the co-regulated genes account only for a relatively small portion of all DAS or DE genes, which indicated that AS and gene activation could be separately regulated.
Conclusions
This work presents the most thorough examination to date of the transcriptomes of S. japonicum larval stages. Evidence of DNA replication and cell division was only seen and confirmed in the sporocysts, while each stage upregulated different genes involved in development, morphogenesis, movement and host invasion. Our data indicated that neprilysins and leishmanolysins might play a role in the penetration of the snail by the miracidia. It’s known that S. japonicum cathepsin B2 (SjCB2) played fundamental roles in skin penetration [103]. Our analysis indicated that cercarial elastase (SjCE2b) and leishmanolysins might also be involved in the process of cercariae invasion and the tandem duplications of cathepsin B-like proteases probably contributed to the wide mammalian host range of the S. japonicum. These genes should be targeted in the future for hypothesis-driven functional studies. The expression profile of stem cell markers revealed that different populations of germinal cells exist in the larval stages. We also performed the most comprehensive AS analysis in S. japonicum. We found that the AS prevalence was 61–64% at the genome-wide level, and ES was the predominant AS type in the larval transcriptomes, which revealed an affinity with its mammalian hosts in gene regulation patterns. The transcriptome profiles of S. japonicum larval stages provide new insights into host invasion, and the landscape of AS will not only facilitate future investigations on transcriptome complexity and AS regulation during the life cycle of Schistosoma species but also offered an invaluable resource for future functional and evolutionary studies of AS in platyhelminth parasites.
Supporting information
S1 Fig. A heatmap showing Person’s correlation coefficient among different samples, with correlation levels indicated by colors.
The scores calculated by the R function () indicated the correlation levels between two samples.
https://doi.org/10.1371/journal.pntd.0009889.s001
(TIF)
S2 Fig. Hierarchical clustering analysis (HCA) of transcriptional profiles from 12 S. japonicum samples with 8,732 genes.
Egg, egg; Mir, miracidium; Spo, sporocyst; Cer, cercaria.
https://doi.org/10.1371/journal.pntd.0009889.s002
(TIF)
S3 Fig. Domain organizations of the S. mansoni Omega-1 and the four S. japonicum T2 ribonucleases.
The signal peptide and ribonuclease_T2 domain are depicted in red and blue, respectively.
https://doi.org/10.1371/journal.pntd.0009889.s003
(TIF)
S4 Fig.
GO enrichment for differentially expressed genes (DEGs) in miracidium compared to egg (A), in sporocyst compared to miracidium (B), and in cercaria compared to sporocyst (C).
https://doi.org/10.1371/journal.pntd.0009889.s004
(TIF)
S5 Fig. Phylogenetic relationship of papain in the four Schistosoma species.
The protein and genome sequences of Schistosoma haematobium SchHae_2.0 [104] and Schistosoma bovis ASM395894v1 [105] were downloaded from the WormBase ParaSite (https://parasite.wormbase.org/index.html). An un-rooted phylogenic tree was constructed in MEGA 7 on the basis of multiple alignment of full-sequences from S. japonicum, S. mansoni, S. haematobium and S. bovis. Tandem duplicated cathepsin B-like cysteine proteases of S. japonicum and S. mansoni were indicated by bold black lines.
https://doi.org/10.1371/journal.pntd.0009889.s005
(TIF)
S6 Fig.
Interactions between the five types of detected AS genes in the (A) egg stage, (B) miracidium stage, and (C) sporocyst stage were visualized using an UpSet plot.
https://doi.org/10.1371/journal.pntd.0009889.s006
(TIF)
S7 Fig.
Venn diagram of the overlap of the DE and DAS genes between the (A) miracidium and sporocyst stages, and (B) sporocyst and cercaria stages.
https://doi.org/10.1371/journal.pntd.0009889.s007
(TIF)
S1 Table. Summary of sequence statistics for S. japonicum RNA-seq data.
https://doi.org/10.1371/journal.pntd.0009889.s008
(XLSX)
S2 Table. Illumina RNA-Seq S. japonicum transcript TPM (Transcripts Per Kilobase million) values.
https://doi.org/10.1371/journal.pntd.0009889.s009
(CSV)
S3 Table. List of stage-specific genes (SSG) or stage-enriched genes (SEG) in the four S. japonicum larval stages.
https://doi.org/10.1371/journal.pntd.0009889.s010
(XLSX)
S4 Table. Enriched GO terms of SSG and SEG in the four S. japonicum larval stages.
Gene ratio is the percentage of total SEG or SSG in the given GO term.
https://doi.org/10.1371/journal.pntd.0009889.s011
(XLSX)
S5 Table. Differentially expressed genes (DEG) between the adjacent life stages.
https://doi.org/10.1371/journal.pntd.0009889.s012
(XLSX)
S6 Table. Enriched GO terms of genes differentially expressed between the adjacent life stages.
https://doi.org/10.1371/journal.pntd.0009889.s013
(XLSX)
S7 Table. Enriched GO terms of genes in the eight clusters.
https://doi.org/10.1371/journal.pntd.0009889.s014
(XLSX)
S8 Table. 18 homologs of stem cell markers between S. japonicum and S. mansoni.
https://doi.org/10.1371/journal.pntd.0009889.s015
(CSV)
S1 Dataset. The HMMER output file and multiple sequence alignment results of peptidase C1 (papain) in the four Schistosoma species, and the GFF file of peptidase C1 (papain) genes in the SjV3 genome.
https://doi.org/10.1371/journal.pntd.0009889.s016
(ZIP)
S2 Dataset. AS events in each life stage and the differential AS events between the adjacent life stages.
https://doi.org/10.1371/journal.pntd.0009889.s017
(ZIP)
References
- 1. McManus DP, Dunne DW, Sacko M, Utzinger J, Vennervald BJ, Zhou X-N. Schistosomiasis. Nat Rev Dis Primer. 2018;4: 13. pmid:30093684
- 2. Faust EC. The Reactions of the Miracidia of Schistosoma japonicum and S. haematobium in the Presence of Their Intermediate Hosts. J Parasitol. 1924;10: 199.
- 3. Mouahid G, Rognon A, de Carvalho Augusto R, Driguez P, Geyer K, Karinshak S, et al. Transplantation of schistosome sporocysts between host snails: A video guide. Wellcome Open Res. 2018;3. pmid:29487916
- 4. Iii JJC, King RS, Cogswell A, Williams DL, Newmark PA. An Atlas for Schistosoma mansoni Organs and Life-Cycle Stages Using Cell Type-Specific Markers and Confocal Microscopy. PLoS Negl Trop Dis. 2011;5: e1009. pmid:21408085
- 5. Fung M-C, Lau M-T, Chen X-G. Expressed sequence tag (EST) analysis of a Schistosoma japonicum cercariae cDNA library. Acta Trop. 2002;82: 215–224. pmid:12020895
- 6. Williams DL, Sayed AA, Bernier J, Birkeland SR, Cipriano MJ, Papa AR, et al. Profiling Schistosoma mansoni development using serial analysis of gene expression (SAGE). Exp Parasitol. 2007;117: 246–258. pmid:17577588
- 7. Taft AS, Vermeire JJ, Bernier J, Birkeland SR, Cipriano MJ, Papa AR, et al. Transcriptome analysis of Schistosoma mansoni larval development using serial analysis of gene expression (SAGE). Parasitology. 2009;136: 469–485. pmid:19265565
- 8. Leutner S, Oliveira KC, Rotter B, Beckmann S, Buro C, Hahnel S, et al. Combinatory Microarray and SuperSAGE Analyses Identify Pairing-Dependently Transcribed Genes in Schistosoma mansoni Males, Including Follistatin. PLoS Negl Trop Dis. 2013;7: e2532. pmid:24244773
- 9. Cai P, Liu S, Piao X, Hou N, You H, McManus DP, et al. A next-generation microarray further reveals stage-enriched gene expression pattern in the blood fluke Schistosoma japonicum. Parasit Vectors. 2017;10: 19. pmid:28069074
- 10. Gobert GN, Moertel L, Brindley PJ, McManus DP. Developmental gene expression profiles of the human pathogen Schistosoma japonicum. BMC Genomics. 2009;10: 128. pmid:19320991
- 11. Phuphisut O, Ajawatanawong P, Limpanont Y, Reamtong O, Nuamtanong S, Ampawong S, et al. Transcriptomic analysis of male and female Schistosoma mekongi adult worms. Parasit Vectors. 2018;11: 504. pmid:30201055
- 12. Anderson L, Amaral MS, Beckedorff F, Silva LF, Dazzani B, Oliveira KC, et al. Schistosoma mansoni Egg, Adult Male and Female Comparative Gene Expression Analysis and Identification of Novel Genes by RNA-Seq. PLoS Negl Trop Dis. 2015;9: e0004334. pmid:26719891
- 13. Almeida GT, Amaral MS, Beckedorff FCF, Kitajima JP, DeMarco R, Verjovski-Almeida S. Exploring the Schistosoma mansoni adult male transcriptome using RNA-seq. Exp Parasitol. 2012;132: 22–31. pmid:21745473
- 14. Lu Z, Sessler F, Holroyd N, Hahnel S, Quack T, Berriman M, et al. Schistosome sex matters: a deep view into gonad-specific and pairing-dependent transcriptomes reveals a complex gender interplay. Sci Rep. 2016;6: 31150. pmid:27499125
- 15. Liu S, Zhou X, Piao X, Wu C, Hou N, Chen Q. Comparative Analysis of Transcriptional Profiles of Adult Schistosoma japonicum from Different Laboratory Animals and the Natural Host, Water Buffalo. PLoS Negl Trop Dis. 2015;9: e0003993. pmid:26285138
- 16. Liu R, Cheng W-J, Ye F, Zhang Y-D, Zhong Q-P, Dong H-F, et al. Comparative Transcriptome Analyses of Schistosoma japonicum Derived From SCID Mice and BALB/c Mice: Clues to the Abnormality in Parasite Growth and Development. Front Microbiol. 2020;11: 274. pmid:32218772
- 17. Wang J, Yu Y, Shen H, Qing T, Zheng Y, Li Q, et al. Dynamic transcriptomes identify biogenic amines and insect-like hormonal regulation for mediating reproduction in Schistosoma japonicum. Nat Commun. 2017;8: 14693. pmid:28287085
- 18. Wangwiwatsin A, Protasio AV, Wilson S, Owusu C, Holroyd NE, Sanders MJ, et al. Transcriptome of the parasitic flatworm Schistosoma mansoni during intra-mammalian development. PLoS Negl Trop Dis. 2020;14: e0007743. pmid:32374726
- 19. Buddenborg SK, Kamel B, Hanelt B, Bu L, Zhang S-M, Mkoji GM, et al. The in vivo transcriptome of Schistosoma mansoni in the prominent vector species Biomphalaria pfeifferi with supporting observations from Biomphalaria glabrata. PLoS Negl Trop Dis. 2019;13: e0007013. pmid:31568484
- 20. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34: i884–i890. pmid:30423086
- 21. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29: 15–21. pmid:23104886
- 22. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12: 323. pmid:21816040
- 23. Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Research. 2016;4. pmid:26925227
- 24. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15: 550. pmid:25516281
- 25. Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics J Integr Biol. 2012;16: 284–287. pmid:22455463
- 26. Rawlings ND, Barrett AJ, Thomas PD, Huang X, Bateman A, Finn RD. The MEROPS database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database. Nucleic Acids Res. 2018;46: D624–D632. pmid:29145643
- 27. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, et al. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol Plant. 2020;13: 1194–1202. pmid:32585190
- 28. Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40: e49. pmid:22217600
- 29. Shen S, Park JW, Lu Z, Lin L, Henry MD, Wu YN, et al. rMATS: Robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proc Natl Acad Sci. 2014;111: E5593–E5601. pmid:25480548
- 30. Luhtala N, Parker R. T2 Family Ribonucleases: Ancient enzymes with diverse roles. Trends Biochem Sci. 2010;35: 253–259. pmid:20189811
- 31. Ittiprasert W, Mann VH, Karinshak SE, Coghlan A, Rinaldi G, Sankaranarayanan G, et al. Programmed genome editing of the omega-1 ribonuclease of the blood fluke, Schistosoma mansoni. eLife. 2019;8: e41337. pmid:30644357
- 32. Mouser EE, Pollakis G, Smits HH, Thomas J, Yazdanbakhsh M, Jong EC de, et al. Schistosoma mansoni soluble egg antigen (SEA) and recombinant Omega-1 modulate induced CD4+ T-lymphocyte responses and HIV-1 infection in vitro. PLOS Pathog. 2019;15: e1007924. pmid:31487324
- 33. Everts B, Perona-Wright G, Smits HH, Hokke CH, van der Ham AJ, Fitzsimmons CM, et al. Omega-1, a glycoprotein secreted by Schistosoma mansoni eggs, drives Th2 responses. J Exp Med. 2009;206: 1673–1680. pmid:19635864
- 34. Steinfelder S, Andersen JF, Cannons JL, Feng CG, Joshi M, Dwyer D, et al. The major component in schistosome eggs responsible for conditioning dendritic cells for Th2 polarization is a T2 ribonuclease (omega-1). J Exp Med. 2009;206: 1681–1690. pmid:19635859
- 35. Maia TM, Gogendeau D, Pennetier C, Janke C, Basto R. Bug22 influences cilium morphology and the post-translational modification of ciliary microtubules. Biol Open. 2014;3: 138–151. pmid:24414207
- 36. Talaga AK, Dong FN, Reisert J, Zhao H. Cilia- and Flagella-Associated Protein 69 Regulates Olfactory Transduction Kinetics in Mice. J Neurosci. 2017;37: 5699. pmid:28495971
- 37. Wright CA. Host-Location by Trematode Miracidia. Ann Trop Med Parasitol. 1959;53: 288–292. pmid:13846053
- 38. Weiger WA. Serotonergic modulation of behaviour: a phylogenetic overview. Biol Rev Camb Philos Soc. 1997;72: 61–95. pmid:9116165
- 39. Boyle JP, Zaide JV, Yoshino TP. Schistosoma mansoni: Effects of Serotonin and Serotonin Receptor Antagonists on Motility and Length of Primary Sporocysts in Vitro. Exp Parasitol. 2000;94: 217–226. pmid:10831389
- 40. Boyle JP, Yoshino TP. Serotonin-Induced Muscular Activity in Schistosoma mansoni Larval Stages: Importance of 5-HT Transport and Role in Daughter Sporocyst Production. J Parasitol. 2005;91: 542–550. pmid:16108544
- 41. Patocka N, Sharma N, Rashid M, Ribeiro P. Serotonin Signaling in Schistosoma mansoni: A Serotonin–Activated G Protein-Coupled Receptor Controls Parasite Movement. PLoS Pathog. 2014;10. pmid:24453972
- 42. Ingram JR, Rafi SB, Eroy-Reveles AA, Ray M, Lambeth L, Hsieh I, et al. Investigation of the Proteolytic Functions of an Expanded Cercarial Elastase Gene Family in Schistosoma mansoni. PLoS Negl Trop Dis. 2012;6: e1589. pmid:22509414
- 43. Zhou Y, Zheng H, Chen Y, Zhang L, Wang K, Guo J, et al. The Schistosoma japonicum genome reveals features of host–parasite interplay. Nature. 2009;460: 345–351. pmid:19606140
- 44. Ingram J, Knudsen G, Lim KC, Hansell E, Sakanari J, McKerrow J. Proteomic Analysis of Human Skin Treated with Larval Schistosome Peptidases Reveals Distinct Invasion Strategies among Species of Blood Flukes. PLoS Negl Trop Dis. 2011;5: e1337. pmid:21980548
- 45. Zhang T, Mo X-J, Xu B, Yang Z, Gobert GN, Yan S, et al. Enzyme activity of Schistosoma japonicum cercarial elastase SjCE-2b ascertained by in vitro refolded recombinant protein. Acta Trop. 2018;187: 15–22. pmid:30017495
- 46. Neto ALC, Brito ANALM, Rezende AM, Magalhães FB, Neto OP de M. In silico characterization of multiple genes encoding the GP63 virulence protein from Leishmania braziliensis: identification of sources of variation and putative roles in immune evasion. BMC Genomics. 2019;20. pmid:30732584
- 47. Hambrook JR, Kaboré AL, Pila EA, Hanington PC. A metalloprotease produced by larval Schistosoma mansoni facilitates infection establishment and maintenance in the snail host by interfering with immune cell function. PLoS Pathog. 2018;14: e1007393. pmid:30372490
- 48. Hagerty JR, Jolly ER. Heads or tails? Differential translational regulation in cercarial heads and tails of schistosome worms. PLOS ONE. 2019;14: e0224358. pmid:31658287
- 49. Y S H I. Visual pigment: G-protein-coupled receptor for light signals. Cell Mol Life Sci CMLS. 1998;54: 1299–1315. pmid:9893707
- 50. Burrows M, Pfluger HJ. Action of locust neuromodulatory neurons is coupled to specific motor patterns. J Neurophysiol. 1995;74: 347–357. pmid:7472336
- 51. El-Sakkary N, Chen S, Arkin MR, Caffrey CR, Ribeiro P. Octopamine signaling in the metazoan pathogen Schistosoma mansoni: localization, small-molecule screening and opportunities for drug development. Dis Model Mech. 2018;11. pmid:29925529
- 52. Bhattacharya A, Padhan N, Jain R, Bhattacharya S. Calcium-binding proteins of Entamoeba histolytica. Arch Med Res. 2006;37: 221–225. pmid:16380322
- 53. Ikeda T. Effects of blockers of Ca2+ channels and other ion channels on in vitro excystment of Paragonimus ohirai metacercariae induced by sodium cholate. Parasitol Res. 2004;94: 329–331. pmid:15449179
- 54. Leontovyč R, Young ND, Korhonen PK, Hall RS, Tan P, Mikeš L, et al. Comparative Transcriptomic Exploration Reveals Unique Molecular Adaptations of Neuropathogenic Trichobilharzia to Invade and Parasitize Its Avian Definitive Host. PLoS Negl Trop Dis. 2016;10: e0004406. pmid:26863542
- 55. Dorsey CH, Stirewalt MA. Schistosoma mansoni: localization of calcium-detecting reagents in electron-lucent areas of specific preacetabular gland granules. Z Parasitenkd Berl Ger. 1977;54: 165–173. pmid:605648
- 56. Ram D, Grossman Z, Markovics A, Avivi A, Ziv E, Lantner F, et al. Rapid changes in the expression of a gene encoding a calcium-binding protein in Schistosoma mansoni. Mol Biochem Parasitol. 1989;34: 167–175. pmid:2710168
- 57. Bianchi V, Spychala J. Mammalian 5’-nucleotidases. J Biol Chem. 2003;278: 46195–46198. pmid:12947102
- 58. Rampazzo C, Miazzi C, Franzolin E, Pontarin G, Ferraro P, Frangini M, et al. Regulation by degradation, a cellular defense against deoxyribonucleotide pool imbalances. Mutat Res. 2010;703: 2–10. pmid:20561600
- 59. Brito CFA, Oliveira GC, Oliveira SC, Street M, Riengrojpitak S, Wilson RA, et al. Sm14 gene expression in different stages of the Schistosoma mansoni life cycle and immunolocalization of the Sm14 protein within the adult worm. Braz J Med Biol Res. 2002;35: 377–381. pmid:11887217
- 60. Liang J, Han Q, Tan Y, Ding H, Li J. Current Advances on Structure-Function Relationships of Pyridoxal 5′-Phosphate-Dependent Enzymes. Front Mol Biosci. 2019;6. pmid:30891451
- 61. Shimizu T. Heme: emergent roles of heme in signal transduction, functional regulation and as catalytic centres. Chem Soc Rev. 2019; 35. pmid:31748766
- 62. Chiang C-P, Caulfield JP. Schistosoma mansoni: Ultrastructural demonstration of a miracidial glycocalyx that cross-reacts with antibodies raised against the cercarial glycocalyx. Exp Parasitol. 1988;67: 63–72. pmid:2458959
- 63. Bland ND, Pinney JW, Thomas JE, Turner AJ, Isaac RE. Bioinformatic analysis of the neprilysin (M13) family of peptidases reveals complex evolutionary and functional relationships. BMC Evol Biol. 2008;8: 16. pmid:18215274
- 64. Duvaux-Miret O, Stefano GB, Smith EM, Dissous C, Capron A. Immunosuppression in the definitive and intermediate hosts of the human parasite Schistosoma mansoni by release of immunoactive neuropeptides. Proc Natl Acad Sci U S A. 1992;89: 778–781. pmid:1309957
- 65. Grote A, Caffrey CR, Rebello KM, Smith D, Dalton JP, Lustigman S. Cysteine proteases during larval migration and development of helminths in their final host. PLoS Negl Trop Dis. 2018;12: e0005919. pmid:30138448
- 66. Sojka D, Hartmann D, Bartošová-Sojková P, Dvořák J. Parasite Cathepsin D-Like Peptidases and Their Relevance as Therapeutic Targets. Trends Parasitol. 2016;32: 708–723. pmid:27344362
- 67. Yang Y, Wen Y jun, Cai YN, Vallée I, Boireau P, Liu MY, et al. Serine Proteases of Parasitic Helminths. Korean J Parasitol. 2015;53: 1–11. pmid:25748703
- 68. de Bono M, Maricq AV. Neuronal substrates of complex behaviors in C. elegans. Annu Rev Neurosci. 2005;28: 451–501. pmid:16022603
- 69. Vanden Broeck J, Torfs H, Poels J, Van Poyer W, Swinnen E, Ferket K, et al. Tachykinin-like peptides and their receptors. A review. Ann N Y Acad Sci. 1999;897: 374–387. pmid:10676464
- 70. Newmark PA, Alvarado AS. Not your father’s planarian: a classic model enters the era of functional genomics. Nat Rev Genet. 2002;3: 210–219. pmid:11972158
- 71. Wagner DE, Wang IE, Reddien PW. Clonogenic neoblasts are pluripotent adult stem cells that underlie planarian regeneration. Science. 2011;332: 811–816. pmid:21566185
- 72. You H, Jones MK, Whitworth DJ, McManus DP. Innovations and Advances in Schistosome Stem Cell Research. Front Immunol. 2021;12. pmid:33746946
- 73. Wendt GR, Collins JJ. Schistosomiasis as a disease of stem cells. Curr Opin Genet Dev. 2016;40: 95–102. pmid:27392295
- 74. Jourdane J, Theron A. Schistosoma mansoni: cloning by microsurgical transplantation of sporocysts. Exp Parasitol. 1980;50: 349–357. pmid:7428911
- 75. Collins JNR, Collins JJ. Methods for Studying the Germline of the Human Parasite Schistosoma mansoni. Methods Mol Biol Clifton NJ. 2017;1463: 35–47. pmid:27734345
- 76. Wang B, Collins JJ III, Newmark PA. Functional genomic characterization of neoblast-like stem cells in larval Schistosoma mansoni. eLife. 2013;2: e00768. pmid:23908765
- 77. Wang B, Lee J, Li P, Saberi A, Yang H, Liu C, et al. Stem cell heterogeneity drives the parasitic life cycle of Schistosoma mansoni. eLife. 2018;7: e35449. pmid:29988015
- 78. Diaz Soria CL, Lee J, Chong T, Coghlan A, Tracey A, Young MD, et al. Single-cell atlas of the first intra-mammalian developmental stage of the human parasite Schistosoma mansoni. Nat Commun. 2020;11: 6411. pmid:33339816
- 79. Collins Iii JJ, Wang B, Lambrus BG, Tharp ME, Iyer H, Newmark PA. Adult somatic stem cells in the human parasite Schistosoma mansoni. Nature. 2013;494: 476–479. pmid:23426263
- 80. Dissous C, Grevelding CG, Long T. Schistosoma mansoni Polo-like kinases and their function in control of mitosis and parasite reproduction. An Acad Bras Cienc. 2011;83: 627–635. pmid:21670883
- 81. Wendt GR, Collins JN, Pei J, Pearson MS, Bennett HM, Loukas A, et al. Flatworm-specific transcriptional regulators promote the specification of tegumental progenitors in Schistosoma mansoni. eLife. 2018;7: e33221. pmid:29557781
- 82. Taft AS, Yoshino TP. CLONING AND FUNCTIONAL CHARACTERIZATION OF TWO CALMODULIN GENES DURING LARVAL DEVELOPMENT IN THE PARASITIC FLATWORM SCHISTOSOMA MANSONI. J Parasitol. 2011;97: 72–81. pmid:21348610
- 83. Katsumata T, Kohno S, Yamaguchi K, Hara K, Aoki Y. Hatching of Schistosoma mansoni eggs is a Ca2+/calmodulin-dependent process. Parasitol Res. 1989;76: 90–91. pmid:2622901
- 84. Siqueira-Neto JL, Debnath A, McCall L-I, Bernatchez JA, Ndao M, Reed SL, et al. Cysteine proteases in protozoan parasites. PLoS Negl Trop Dis. 2018;12: e0006512. pmid:30138453
- 85. El Ridi R, Tallima H, Selim S, Donnelly S, Cotton S, Gonzales Santana B, et al. Cysteine Peptidases as Schistosomiasis Vaccines with Inbuilt Adjuvanticity. PLoS ONE. 2014;9: e85401. pmid:24465551
- 86. Wang L, Bei X, Gao J, Li Y, Yan Y, Hu Y. The similar and different evolutionary trends of MATE family occurred between rice and Arabidopsis thaliana. BMC Plant Biol. 2016;16: 207. pmid:27669820
- 87. Zhang J. Evolution by gene duplication: an update. Trends Ecol Evol. 2003;18: 292–298.
- 88. Zhou X-N, Bergquist R, Leonardo L, Yang G-J, Yang K, Sudomo M, et al. Schistosomiasis japonica control and research needs. Adv Parasitol. 2010;72: 145–178. pmid:20624531
- 89. Ramani AK, Calarco JA, Pan Q, Mavandadi S, Wang Y, Nelson AC, et al. Genome-wide analysis of alternative splicing in Caenorhabditis elegans. Genome Res. 2011;21: 342. pmid:21177968
- 90. Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. Alternative isoform regulation in human tissue transcriptomes. 2008;456: 7. pmid:18978772
- 91. Fagnani M, Barash Y, Ip JY, Misquitta C, Pan Q, Saltzman AL, et al. Functional coordination of alternative splicing in the mammalian central nervous system. Genome Biol. 2007;8: R108. pmid:17565696
- 92. Hull R, Dlamini Z. The role played by alternative splicing in antigenic variability in human endo-parasites. Parasit Vectors. 2014;7: 53. pmid:24472559
- 93. DeMarco R, Mathieson W, Manuel SJ, Dillon GP, Curwen RS, Ashton PD, et al. Protein variation in blood-dwelling schistosome worms generated by differential splicing of micro-exon gene transcripts. Genome Res. 2010;20: 1112–1121. pmid:20606017
- 94. Piao X, Hou N, Cai P, Liu S, Wu C, Chen Q. Genome-wide transcriptome analysis shows extensive alternative RNA splicing in the zoonotic parasite Schistosoma japonicum. BMC Genomics. 2014;15: 715. pmid:25156522
- 95. Wang X, Xu X, Lu X, Zhang Y, Pan W. Transcriptome Bioinformatical Analysis of Vertebrate Stages of Schistosoma japonicum Reveals Alternative Splicing Events. PloS One. 2015;10: e0138470. pmid:26407301
- 96. Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008;40: 1413–1415. pmid:18978789
- 97. Chaudhary S, Khokhar W, Jabre I, Reddy ASN, Byrne LJ, Wilson CM, et al. Alternative Splicing and Protein Diversity: Plants Versus Animals. Front Plant Sci. 2019;10. pmid:31244866
- 98. Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456: 470–476. pmid:18978772
- 99. Resch AM, Palakodeti D, Lu Y-C, Horowitz M, Graveley BR. Transcriptome Analysis Reveals Strain-Specific and Conserved Stemness Genes in Schmidtea mediterranea. PLoS ONE. 2012;7: e34447. pmid:22496805
- 100. Liu S, Zhou X, Hao L, Piao X, Hou N, Chen Q. Genome-Wide Transcriptome Analysis Reveals Extensive Alternative Splicing Events in the Protoscoleces of Echinococcus granulosus and Echinococcus multilocularis. Front Microbiol. 2017;8: 929. pmid:28588571
- 101. Jung H, Han S, Lee Y. Transcriptome analysis of alternative splicing in the pathogen life cycle in human foreskin fibroblasts infected with Trypanosoma cruzi. Sci Rep. 2020;10. pmid:33060827
- 102. Ding F, Cui P, Wang Z, Zhang S, Ali S, Xiong L. Genome-wide analysis of alternative splicing of pre-mRNA under salt stress in Arabidopsis. BMC Genomics. 2014;15: 431. pmid:24897929
- 103. Zhu B, Luo F, Shen Y, Yang W, Sun C, Wang J, et al. Schistosoma japonicum cathepsin B2 (SjCB2) facilitates parasite invasion through the skin. PLoS Negl Trop Dis. 2020;14: e0008810. pmid:33104723
- 104. Stroehlein AJ, Korhonen PK, Chong TM, Lim YL, Chan KG, Webster B, et al. High-quality Schistosoma haematobium genome achieved by single-molecule and long-range sequencing. GigaScience. 2019;8. pmid:31494670
- 105. Oey H, Zakrzewski M, Gravermann K, Young ND, Korhonen PK, Gobert GN, et al. Whole-genome sequence of the bovine blood fluke Schistosoma bovis supports interspecific hybridization with S. haematobium. PLOS Pathog. 2019;15: e1007513. pmid:30673782