Skip to content


  • Research article
  • Open Access

Comprehensive analysis of miRNA profiles reveals the role of Schistosoma japonicum miRNAs at different developmental stages

Veterinary Research201950:23

  • Received: 1 November 2018
  • Accepted: 19 February 2019
  • Published:


Schistosomiasis is an important zoonotic disease affecting up to 40 kinds of animals and 250 million people. It has been reported that the miRNAs play a role in the metabolism, differentiation, development and reproduction in many organisms. However, the roles of miRNAs regulating the development, maturation and production in schistosome in both females and males remains unclear. Here we present the dynamic transcriptome analysis of all 79 known Schistosoma japonicum miRNAs from pairing to production, including 14 days post-infection (dpi), 16, 18, 20, 22, 24, 26, 28 dpi female and male, by small RNA sequencing. The miRNA expression profiles showed time-related characteristics in male and female from paring to production, which could be clustered into three patterns, characterized by pairing stage highly expressed (cluster 1), maturating stage highly expressed (cluster 2), and egg producing stage highly expressed (cluster 3). The enrichment of miRNA cluster targeted genes in female and male were distinctly different. Network analysis of miRNAs and their target regulation showed that cluster 1 had 15 miRNAs involved in the regulation of interaction, communication, immune response in female–male and parasite–host. The other 11 miRNAs were involved in gender differentiation and the meiotic cell cycle process. In cluster 2, 11 miRNAs were involved in development and sexual maturation. In cluster 3, 45 miRNAs possibly regulate metabolism and synthesis of the substance for egg production. Analysis of the miRNA regulation network would contribute to understanding the molecular mechanism in S. japonicum development and egg production.


Schistosomiasis, a zoonotic parasitic disease, is caused by blood flukes of Schistosoma spp. Schistosomiasis has a wide range of definitive hosts, not only more than 40 kinds of animals including cattle, goats, sheep, camels, canine, swine, equine, equus ferus × asinus, macaca, felidae and rodents and avian [14], but also affecting 258 million people in 76 countries in tropical and subtropical regions. Schistosoma japonicum (S. japonicum) is prevalent in Asia, primarily China and the Philippines [5, 6]. After infection, schistosomula migrate into the host’s bloodstream, where they develop into sexual maturation and lay eggs. Dioecy is a unique characteristic of schistosomes, in which the adult male and female live in a copula, while the female does not mature without a male [79]. Each paired female S. japonicum can deposit 1500–3000 eggs per day, which is about 5 to 10 times more than that of Schistosoma mansoni (S. mansoni). Thus, S. japonicum causes more severe morbidity than other schistosome species [9, 10]. Schistosome sex differentiation, maturation and egg production are the key events contributing to consequent morbidity and transmission. Understanding the molecular mechanisms of these processes is extremely urgent in current investigations on trematodes. MicroRNAs (miRNAs) are a class of small non-coding RNAs (~22 nucleotides in length), which play a crucial role in regulation of gene expression by binding to target messenger RNA (mRNA) and triggering translation repression or mRNA degradation [11, 12].

MicroRNAs have been reported to be involved in post-transcriptional regulation of gene expression during development, differentiation, proliferation, death and metabolism in a variety of organisms [13]. To date an increasing number of schistosome miRNAs from various development stages have been identified with transcriptome analysis and deep-sequencing technique in S. japonicum [1418] and S. mansoni [1921]. miRNA and endogenous siRNA expression in S. japonicum were profiled at the stages of cercariae, lung-stage schistosomula, eggs sequestered in tissues, adults of both sexes by deep sequencing and qRT-PCR methods, which provided a broader view of small RNAs of the parasite [22, 23]. A study on sex differential expression of miRNAs in S. mansoni found that 10 of 13 microRNAs were more abundant in females than in males, exhibiting sex-biased expression patterns, and thus pointing towards a plausible involvement of miRNAs in the sex differentiation and maintenance [21]. Sun et al. investigated the differential expression profiles of S. japonicum miRNAs in mature and immature worms, and they also analyzed miRNAs’ target genes. The results showed that more genes and metabolic pathways were regulated by miRNAs in paired mature females than in unpaired immature females [17]. Sequencing of small RNAs from 16, 22 and 28 dpi S. japonicum showed that 18 miRNAs are expressed differently between male and female, and that by knocking down the expression of female enriched miRNAs sja-bantam and sja-miR-31, morphological alterations were observed in female ovaries [18]. Recently, Protasio et al. reported that a novel S. mansoni miRNA family (sma-miR-277/4989) might play a dominant role in post-translational regulation in the sexual development of female schistosomes [24]. Wang’s report [25] showed that S. japonicum pairing began at as early as 16 dpi. At 20 dpi, oogonium and spermatogonium were initially differentiated and mature vitelline cells firstly appeared in females. At 24 dpi, mature germ cells and vitelline glands had appeared in most of females. Meanwhile, at 20–22 dpi and 24–26 dpi, thousands of transcripts were down regulated in female worms. Therefore, 14–16 dpi, 20–24 dpi, 24–26 dpi were crucial stages for S. japonicum pairing, development and maturation [25]. However, all known 79 S. japonicum miRNAs expression profiles and the role of miRNAs during those developmental stages (14–16 dpi, 20–24 dpi, 24–26 dpi) are currently unknown.

To gain insight into the role of S. japonicum miRNAs in the developmental process, particularly in gender differentiation and maturation, we utilized high-throughput small RNA sequencing to quantify the dynamic expression profile of miRNAs at 14, 16, 18, 20, 22, 24, 26, 28 dpi S. japonicum worms. Our results demonstrated that 79 known miRNAs had time-related characteristics in male and female. Enrichments of each miRNA cluster targets in female and male were distinctly different. Network analysis showed different miRNAs might play critical roles in the pairing, developmental and egg production stages. Our results also showed that different transferase activities were enriched in F-cluster 3, M-cluster 1 and M-cluster 3, which suggested that miRNAs might participate in the post-translational modification during different developmental stages. Global view and specific understanding of miRNA roles would contribute to study the molecular mechanism of S. japonicum developmental and maturation events.

Materials and methods

Animals, parasite and worm collection

Schistosoma japonicum was harvested by hepatic-portal perfusion from mice (C57BL/6, female, 5–6 weeks) infected with cercariae released from Oncomelania hupensis snails (Anhui isolate). The infected snails were provided by the National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention, Shanghai. Worms were collected at 8 time points on every other day from 14 dpi through to 28 dpi. The sampled worms were washed with sterile DMEM (Dulbecco’s Modified Eagle’s medium) three times. Male and female worms were separated under light microscopy and stored in RNAlater stabilization solution (Invitrogen, USA) at −80 °C before use. Assays were processed with worms from one mouse representing one biological replicate, all assays in our study had three biological replicates.

RNA extraction, small RNA library construction and sequencing

Schistosoma japonicum total RNA was isolated using Qiagen miRNeasy Mini Kit (Valencia, CA, USA) according to the manufacturer’s instruction. The concentration and purity were evaluated spectrophotometrically using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and a NanoDrop ND2000 spectrophotometer at 260 nm and 280 nm. The small RNA library was constructed as previously described [2628]. Briefly, 10 µg total RNA was separated by electrophoresis on a Novex 15% TBE-Urea gel (Invitrogen Co. Ltd) and RNA in length of 18–30 nt was purified and ligated to 3′ and 5′ adapters. Then for reverse transcription, we used a RT-PCR kit (Invitrogen Co. Ltd) to produce complementary DNA (cDNA). Finally, the cDNA was amplified using adapter primers, the products were purified using a 6% TBE PAGE gel (Invitrogen Co. Ltd) and sequenced with Illumina Genome Analyzer (Illumina HiSeq 2000) at the BGI (Beijing Genomics Institute, Shenzhen, China).

Analysis of sequencing data

The quality of the RNA-Seq data was first examined using the package FastQC [29]. After a trimming process using Trimmomatic [30] adapters, PCR sequences and bases with low quality in the reads were clipped. A sequence read was discarded if it was below 17 bases in length. Using miRDeep v.2.0 [31], reads were aligned to the reference S. japonicum genome [32] from wormbase and known miRNAs database (downloaded from miRBase v20 [33]). Candidate miRNAs were identified with the miRDeep2 module using the core algorithm in RNAfold tool [31] to predict RNA secondary structures and evaluate the structure and signature of each putative miRNA precursor. Predicted miRNAs with significant Rand fold α level (P < 0.05) and miRDeep2 score (≥ 5) were considered as candidate miRNAs. Normalized expression profile (log scale) was generated by the quantifier module in miRDeep2 v.2.0.

Validation of miRNAs expression by qRT-PCR

RNA extraction, and reverse transcription

Total RNA of S. japonicum male or female worms were isolated using Trizol reagent (Life Technologies, USA). Genome DNA was removed from total RNA with the DNase I (Takara, Japan). For polyadenylation and reverse transcription, poly-A tails were added to total RNA (1 μg) uncontaining genome DNA with ATP by Escherichia coli (E. coli) poly (A) polymerase using miScript PCR Starter Kit (Qiagen, Germany). Controls without reverse transcriptase were included for each sample.

Real-time PCR amplification of miRNAs

We applied poly-A quantitative reverse transcription PCR (qRT-PCR) method [34] to amplify miRNAs. For each PCR amplification, approximately 50 ng cDNA was mixed with 10 μM forward primer and 10 μM reverse primer (universal primer) (Qiagen, Germany), 10 μL 2× QuantiTect SYBR Green qPCR master mix (Qiagen, Germany) and sterile water in a final volume of 20 μL. Cycling conditions were: 95 °C for 5 min, followed by 45 cycles of 95 °C for 15 s, 55 ± 2 °C for 30 s and 72 °C for 30 s. Dissociation curve analysis was carried out at the end of each PCR run to verify amplification specificity for each miRNA. All reactions were run in triplicate, each sample has three technical parallel replicates and each group has three biological replicates. Sj-U6 snRNA was used in the qPCR reaction as an endogenous reference with forward primer 5′-GCAAGGATGACACGCAAAT-3′ and reverse primer 5′-ATGGAACGCTTCACGAAT-3′ [35]. Data analysis for Real-time PCR amplification efficiency and for quantifying relative miRNA expression were calculated using the 2−ΔΔCT method [36].

MiRNA target prediction

We used two software packages, RNAhybrid [37] and miRanda [38], to predict miRNA targets. In this study, we predicted the potential miRNA target sites not only with 3′ untranslated region (3′ UTR), but also with coding sequence (CDS) and 5′ UTR as the Ref. [11, 39, 40]. The mRNAs were based off the latest dataset of S. japonicum whole transcriptome (GenBank accession no. PRJNA343582). The cutoff was set: RNAhybrid with e-value ≤ 1e−20, minimum free energy (MFE) ≤ −20 kcal/mol. miRanda was used with a cutoff score of ≥ 140, energy cutoff ≤  −20 kcal/mol, gap opening = − 9.0 and gap extension = −4.0. Target score (i.e. the overall score of a miRNA-3′ UTRs pair in case that multiple binding sites occur) is the sum of scores for all binding sites. The miRNA potential targeted genes predicted by both tools were selected for the following network analysis.

Luciferase system reporter assay and cell culture

MiRNA mimics and NC mimics (2′-O-methyl oligonucleotides) were synthesized chemically by Genepharma (Shanghai, China), the miRNA target binding site fragments were inserted into the multiple cloning sites of pmirGLO vector dual-luciferase report system (Promega, USA). HEK293T cells were cultured in 10% fetal bovine serum in DMEM (Gibco, USA), 800 nanograms of recombinant plasmid and 20 pmol miRNA mimics (GenePharma, China) were co-transfected into HEK293T cells by Lipofectamine 2000 (Invitrogen, USA) in 24-well plates. Firefly and renilla luciferase activity were measured after 24 h transfection by the dual luciferase assay system (Promega, USA).

MiRNA cluster and GO analysis

The miRNA clusters were processed by MeV 4.9 software based on their normalized expression level. MiRNAs expression profile in female were patterned into three clusters: F-cluster 1, F-cluster 2 and F-cluster 3. MiRNAs expression profile in male also were patterned into three clusters: M-cluster 1, M-cluster 2 and M-cluster 3.

GO IDs of target genes were obtained with the BioProject no. PRJNA343582. The Gene Ontology (GO) enrichment analysis of the 6 miRNA clusters were processed by Blast2GO [41] (v4.2). All potential targets from each miRNA cluster in female and male were pooled together and enrichment performed for male and female separately.

MiRNAs network analysis

Firstly, we sorted out GO terms enriched in female and male clusters, respectively. Secondly, we listed the targets enriched in each GO enrichment. Thirdly, we sorted out the miRNAs regulating those targets (from RNAhybrid and miRanda prediction) in each GO enrichment in female and male cluster. Lastly, Network was used to show the relationship between enrichment and miRNAs in cytoscape3.6.0 [42].


Characteristics of S. japonicum small RNAs sequencing

To identify and profile S. japonicum miRNAs expression from pairing, maturation to egg production, we deep sequenced small RNAs in 48 S. japonicum samples (two genders, eight developmental stages, in triplicate). Flow charts of sequencing data analysis are shown in Additional file 1. Results (Figure 1, Additional file 2) showed the highest proportion of total small RNAs were miRNAs (~70%) at 20–26 dpi. The miRNA proportion was around 30–50% at 14–18 dpi and 28 dpi. Other small RNAs were composed by repeat sequences (3.87–23.39%), unannotated sequences (11.68–41.63%), intron (sense) (1.23–4.74%), intron (antisense) (0.44–1.90%), exon (sense) (0.37–1.88%), exon (antisense) (0.28–1.29%), rRNA (0.12–4.43%), tRNA (0.22–3.83%), snRNA (< 0.01%), snoRNA (< 0.01%) (Figure 1). The unannotated sequences might be endogenous siRNAs (small interfering RNAs) derived from the natural antisense transcripts (NAT) [22]. The ratios of small RNAs mapping to the genome were around 75%, ratios of mapping to miRNAs database were around 60% (Additional file 3A). Length of small RNAs were mainly between 20 and 23 nt (Additional file 3B). Our study detected all the 79 known S. japonicum miRNAs (miRBase 20.0 dataset) expression profiles from 14 to 28 dpi in different genders, as well as predicted 108 candidate miRNAs (Additional file 4) from unannotated sequences, in which two candidate miRNAs had the homologues with S. mansoni miRNAs (Additional file 4). The 108 candidate miRNAs still needed to be validated with experiments, e.g.: northern blot, stem-loop qRT-PCR [23, 35, 43].
Figure 1
Figure 1

Characterization of S. japonicum small RNA sequencing. The percentages of miRNAs, rRNA, tRNA, snRNA, snoRNA, repeat, intron (sense), intron (antisense), exon (sense), exon (antisense), unannotated sequences in small RNA sequencing at eight developmental stages (14 dpi, 16 dpi, 18 dpi, 20 dpi, 22 dpi, 24 dpi, 26 dpi, 28 dpi). A Female, B male. rRNA: ribosomal RNA; tRNA: transfer RNA; snRNA: small nuclear RNA; snoRNA: small nucleolar RNAs.

Dynamic profiles of S. japonicum miRNAs in female and male

The expression profiles of 79 known miRNAs (Figure 2, Additional file 5) showed gender and stage differentiation. We validated 12 randomly selected miRNAs for expression by qRT-PCR, which showed miRNA relative expression by qRT-PCR was consistent with the profiles showed by miRNA sequencing (Additional files 6 and 7). The expression profile of all 79 known miRNAs (Figure 2) were grouped into three clusters in each gender (Table 1). In the female cluster 1 (F-cluster 1) (Figure 2B, C), the expression of 20 miRNAs elevated continuously from 14 to 18 dpi, then declined from 20 to 28 dpi. In F-cluster 2 (Figure 2B and C), 25 miRNAs maintained high expression from 18 to 22 dpi. In F-cluster 3 (Figure 2B and C), most of 34 miRNAs presented an abrupt high expression profile at 26 dpi and 28 dpi. In the male expression profile (Figure 2D) the miRNA expression was also grouped into three similar clusters. In M-cluster 1 (Figure 2E and F) 37 miRNAs were expressed at high levels from 14 to 20 dpi. In M-cluster 2 (Figure 2E and F), 22 miRNAs were expressed at high levels from 20 to 26 dpi. For M-cluster 3 (Figure 2E and F), most of those 20 miRNAs were expressed at low levels from 14 to 26 dpi while at 28 dpi their expression increased abruptly. These results provided a straightforward overview of similar miRNAs expression patterns of cluster 1, 2 and 3 in male and female. Cluster 1 was corresponding to the pairing stage (14–20 dpi), cluster 2 to the germ cell development and sexual maturation stage (20–26 dpi) and cluster 3 to egg production stage (24–28 dpi).
Figure 2
Figure 2

Heatmap of 79 S. japonicum known miRNAs expression profile. MiRNAs were significantly high expression (red) or low expression (green). A Heatmap of 79 S. japonicum known miRNAs expression profile in female; B Three clusters of 79 S. japonicum known miRNAs expression profile in female; C Fitting chart of three clusters miRNAs expression profile in female; D Heatmap of 79 S. japonicum known miRNAs expression profile in male; E Three clusters of 79 S. japonicum known miRNAs expression profile in male; F Fitting chart of three clusters miRNAs expression profile in male.Cluster 1: miRNAs had high expression from 14 to 18 dpi. Cluster 2: miRNAs had high expression from 18 to 22 dpi. Cluster 3: miRNAs had high expression from 26 to 28 dpi. MiRNAs expression profile (log scale) data was normalized reads, which were generated by the quantifier module in miRDeep2 from sequencing reads and formed by the fold change formula.

Table 1

The 79 known S. japonicum miRNAs clusters in female and male

MiRNAs cluster



Cluster 1

Total: 20

sja-miR-10-5p, sja-miR-124-5p, sja-miR-2162-5p, sja-miR-2c-5p, sja-miR-2c-3p, sja-miR-2d-5p, sja-miR-3479-5p, sja-miR-3482-5p, sja-miR-3484-3p, sja-miR-3494, sja-miR-7-3p, sja-miR-190-3p, sja-miR-2a-3p, sja-miR-2e-5p, sja-miR-310, sja-miR-3483-3p, sja-miR-3498, sja-miR-3501, sja-miR-3503, sja-miR-71b-3p

Total: 37

sja-miR-2b-3p, sja-miR-2c-3p, sja-miR-3485-5p, sja-miR-3485-3p, sja-miR-3489, sja-miR-3491, sja-miR-3498, sja-miR-3500, sja-miR-71b-3p, sja-miR-2162-5p, sja-miR-2d-3p, sja-miR-2e-5p, sja-miR-3483-3p, sja-miR-3486-5p, sja-miR-277, sja-miR-3490, sja-miR-3492, sja-miR-3493, sja-miR-3494, sja-miR-3497, sja-miR-3499, sja-miR-3501, sja-miR-3502, sja-miR-3504, sja-miR-3506, sja-miR-36-5p, sja-miR-71b-5p, sja-miR-7-3p, sja-miR-124-5p, sja-miR-190-3p, sja-miR-2a-3p, sja-miR-2c-5p, sja-miR-2d-5p, sja-miR-310, sja-miR-3479-5p, sja-miR-3481-5p, sja-miR-3485

Cluster 2

Total: 25

sja-let-7, sja-miR-219-5p, sja-miR-2a-5p, sja-miR-2e-3p, sja-miR-31-5p, sja-miR-31-3p, sja-miR-3507, sja-miR-36-5p, sja-miR-10-3p, sja-miR-124-3p, sja-miR-133, sja-miR-1, sja-miR-219-3p, sja-miR-2b-3p, sja-miR-2d-3p, sja-miR-307, sja-miR-3480-5p, sja-miR-3482-3p, sja-miR-3484-5p, sja-miR-36-3p, sja-miR-61, sja-miR-71a, sja-miR-71b-5p, sja-miR-8-5p, sja-miR-8-3p

Total: 22

sja-miR-124-3p, sja-miR-125a, sja-miR-190-5p, sja-miR-1, sja-miR-3488, sja-miR-3496, sja-miR-3505, sja-miR-61, sja-miR-8-3p, sja-let-7, sja-miR-10-3p, sja-miR-133, sja-miR-2162-3p, sja-miR-219-5p, sja-miR-219-3p, sja-miR-2e-3p, sja-miR-307, sja-miR-3479-3p, sja-miR-3482-3p, sja-miR-3487, sja-miR-7-5p, sja-bantam

Cluster 3

Total: 34

sja-miR-125b, sja-miR-3481-5p, sja-miR-3485-3p, sja-miR-3486-5p, sja-miR-3486-3p, sja-miR-3487, sja-miR-3488, sja-miR-3489, sja-miR-3492, sja-miR-3495, sja-miR-3499, sja-miR-3506, sja-miR-7-5p, sja-bantam, sja-miR-125a, sja-miR-190-5p, sja-miR-2162-3p, sja-miR-2b-5p, sja-miR-3479-3p, sja-miR-3480-3p, sja-miR-3481-3p, sja-miR-3483-5p, sja-miR-3485-5p, sja-miR-3490, sja-miR-3491, sja-miR-3496, sja-miR-3497, sja-miR-3500, sja-miR-3502, sja-miR-3504, sja-miR-277, sja-miR-3505, sja-miR-8185, sja-miR-3493

Total: 20

sja-miR-10-5p, sja-miR-2a-5p, sja-miR-2b-5p, sja-miR-31-3p, sja-miR-3480-5p, sja-miR-3482-5p, sja-miR-3483-5p, sja-miR-3484-3p, sja-miR-71a, sja-miR-8185, sja-miR-125b, sja-miR-3481-3p, sja-miR-3503, sja-miR-36-3p, sja-miR-8-5p, sja-miR-3480-3p, sja-miR-3507, sja-miR-31-5p, sja-miR-3484-5p, sja-miR-3486-3p

Functional enrichment of predicted miRNA targets

Overlap targets were predicted in RNAhybrid and miRanda considered as the potential targets of miRNA (Additional file 8), we selected 3 miRNAs (kept high expression level during 18 to 22 dpi) which had the potential regulation ability to the reproductive development. For identifying target genes of 3 miRNAs, we selected the genes containing the combine sites with miRNA in the 3′UTR, it is common view of miRNA potential target combine sites. The validation were performed in dual luciferase reporter system (Additional file 9). Then, miRNA targets were used to process the functional enrichments. The major enrichments showed in Figure 3 (female) and Figure 4 (male). In female, as to biological process assignment: F-cluster 1 was mainly enriched in regulation of biosynthetic processes, cellular macromolecule localization, and signal transduction by protein phosphorylation (Figure 3). F-cluster 2 was mainly enriched in metabolic processes. F-cluster 3 was mainly enriched in metabolic process and biosynthetic process (Figure 3). In males, as to the biological process assignment: In male, M-cluster 1 was mainly enriched in sex differentiation, meiotic nuclear division, regulation of viral process, regulation of symbiosis (encompassing mutualism through parasitism) (Figure 4). M-cluster 2 was mainly enriched in regulation of the biosynthetic, metabolic and catabolic processes. F-cluster 3 was mainly enriched in cellular macromolecule metabolic processes (Figure 4). As to the molecular function assignment, transferase activity was an arresting enrichment, which were enriched in F-cluster 3 (188 targets), M-cluster 1 (79 targets) and M-cluster 3 (49 targets).
Figure 3
Figure 3

GO enrichments of miRNA target genes of clusters in female. Each cluster GO terms were processed with Blast2GO (v4.2) GO terms were divided into three categories: biological process, molecular function, cellular component. The number of targets in each cluster category was shown in bar charts.

Figure 4
Figure 4

GO enrichments of miRNA target genes of clusters in male. Each cluster GO terms were processed with Blast2GO v4.2. GO terms were divided into three categories: biological process, molecular function, cellular component. The number targets in each cluster category were shown in bar charts.

Except for the major enrichment, there were also some specific enrichments in female and male, which were interesting and worthy of being concerned. Detailed information was showed in Table 2. In a summary, enrichments of each cluster in female and male were different. We discriminated and concluded the enrichments in female and male, respectively. Then we listed which miRNAs regulated those enrichments and showed their network in Figure 5 and Additional file 10.
Table 2

Difference of miRNAs cluster targets enrichment between female and male



Specific functional enrichments


Cluster 1 (pairing stage)

Regulation of chemotaxis

Leukocyte proliferation

Cell junction organization

Regulation of epithelial to mesenchymal transition

Cluster 2 (developmental and mature stage)

Response to BMP (bone morphogenetic protein)

Female pregnancy

Spermatid differentiation

Regulation of meiotic cell cycle

Cluster 3 (egg production stage)

Nerve development

Chemosensory behavior


Cluster 1 (pairing stage)

Sex differentiation

Meiotic nuclear division

Cluster 2 (developmental and mature stage)

Regulation of immune effector process

Regulation of cell growth

Regulation of response to biotic stimulus

Cluster 3 (egg production stage)

Cell communication by chemical coupling

Regulation of cholesterol homeostasis

Regulation of protein activation cascade

Body morphogenesis

Figure 5
Figure 5

Network of cluster miRNA-targeted genes GO terms in female and male. The relationship between miRNAs and their targets’ GO enrichment in each cluster were processed in cytoscape3.6.0. Purple, square represented enrichment in female. Blue, square represented enrichment in male. Turquoise, triangle represented miRNA. Cluster 1: pairing stage; cluster 2: developmental and sexual maturation stage; cluster 3: egg production stage.

Network analysis between miRNAs and their targets enrichment

In order to further understand the multiple regulation roles of miRNAs during different developmental stages, we analyzed the network between miRNAs and their target enrichments in cluster 1, cluster 2 and cluster 3 (Figure 5, Additional file 10), respectively. The network showed 31 miRNAs participated in the regulation in cluster 1 (pairing stage). Among them, 20 miRNAs were specifically enriched in females, 7 miRNAs were specifically enriched in males, and 4 miRNAs were enriched both in females and males. Notably in females, 5 miRNAs regulated the leukocyte proliferation, 7 miRNAs regulated the chemotaxis, and 9 miRNAs regulated the signal transduction by protein phosphorylation. These 15 miRNAs (Figure 5, Additional file 10) might play a key role in the stimulation or attraction of female-male pairing, communication and immune response in parasite–host during pairing stage. In males, 7 miRNAs regulated sex differentiation, and 7 miRNAs regulated meiotic nuclear division. It was suggested that miRNAs may be involved in male sexual development regulation at the early paring stage [44]. In cluster 2 (developmental and mature stage), 29 miRNAs participated the developmental and sexual maturation stage. 15 miRNAs were only enriched in females, 10 miRNAs were only enriched in males, and 4 miRNAs were enriched both in females and males. Notably, 6 miRNAs regulated meiotic cell cycle, pregnancy process in females. 5 miRNAs regulated the cell growth in males. In cluster 3 (egg production stage), 45 miRNAs participated in regulation during the egg production stage, the items of enrichments were less in female and male comparing with cluster 1 and 2. However, 34 miRNAs in female and 21 miRNAs in male were enriched in metabolic processes, which was the largest number of miRNAs observed in a process (Figure 5, Additional file 10). These miRNAs might regulate numerous material synthesis and metabolism for egg production. Notably, 2 miRNAs were enriched in chemosensory behavior, 3 miRNAs were enriched in nerve development, 3 miRNAs were enriched in body morphogenesis and 3 miRNAs were enriched in regulation of epithelial to mesenchymal transition. These 7 miRNAs in cluster 3 might play specific role during interplay in female-male and female specialized function of egg production (Figure 5, Additional file 10). Above information may provide many new clues to investigating the regulation mechanism of male and female sexual development, which deserves further attention and experimental validation.


This study focused on S. japonicum 79 known miRNAs profile (14 to 28 dpi), and the results suggested that S. japonicum miRNAs had three similar expression clusters in male and female from pairing, development, sexual maturation to egg production. The stage related miRNA clusters may regulate specific functions and play distinct roles in order to accomplish S. japonicum parasitism in host during pairing, sex-maturation and production. MiRNA is a potential factor in male–female and parasite–host interplay. Our study provided a comprehensive miRNA analysis with the goal of understanding the molecular role for miRNAs during developmental events of the schistosome life cycle.

We also looked at the target functional enrichments of “the spermatid differentiation” in females, which was interesting. It has previously been reported [44] that cells in the male seminal vesicle were identified as early stages of spermatogenesis, which were immature sperm, while the mature sperms were found in female seminal vesicle and oviduct. This phenomenon indicated that the spermatid maturation process might occur in female. Another possibility of this phenomenon was the sperm maturation occurred in male and the mature sperms were stored in the end of excretory bladder and then given to the female. We couldn’t observed the whole process and couldn’t improve this suppose.

Transferase activity are involved in hundreds of different biochemical pathways throughout the whole biological process, and are integral to some of life’s most important processes. Our results in transferase activity enrichments in F-cluster 3 (188, 39%), M-cluster 1 (72, 99%), M-cluster 3 (49, 100%). showed different transferase activity groups enriched in female and male. Results (Additional files 11, 12) showed 53.1% (26 of 49) of the targets were specific in M-cluster 3, 18.1% (34 of 188) targets were specific in F-cluster 3, 33.3% targets (24 of 72) were specific in M-cluster 1. Next, we also analyzed which miRNAs regulated those transferase group, results (Additional file 13) showed 10 miRNAs regulated transferase activity were specific in M-cluster 3, 7 miRNAs were specific in F-cluster 3, 12 miRNAs were specific in M-cluster 1. There were also different subcategories (Additional file 14) of transferase were enriched in F-cluster 3, M-cluster 1 and M-cluster 3. The rich transferase family member of the targets suggested that miRNA can regulate many important biological processes through post-translational modification [4548] which included of methyltransferase, acetyltransferase, glucosyltransferase, galactosyltransferase, phosphoribosyl transferase, aminotransferase, mannosyl transferase, ribosyl transferase, cytidylyl transferase, pyrophospho kinase, et al. Especially, methyltransferase and acetyltransferase were key regulatory factors in epigenetics via DNA methylation, histone methylation and histone acetylation. Okada et al. [49] reported that histone demethylase JHDM2A is critical for Tnp1 and Prm1 transcription and spermatogenesis, O’Carroll et al. [50] reported that histone H3 methyltransferase gene-Suv39h2 that displays testis-specific expression, Tachibana et al. [51] reported G9a histone methyltransferase plays a dominant role in euchromatic histone H3 lysine 9 methylation and is essential for early embryogenesis. Disorder in histone methylation or acetylation during spermatogenesis directly affected the establishment and maintenance of epigenetic modifications, leading to spermatogenic cell abnormalities and even infertility. Even though those reports were all carried out in mammalian systems, our results indicated methyltransferase activity regulated by miRNA in schistosome might play important developmental role via post-translational modification, e.g. methylation and acetylation.


Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

WH conceived and designed this study, JY performed the analysis of data, experimental assay and drafted the manuscript. YY, HS, RZ and QL participated in the data process, MC participated in the animal experiments and miRNAs sequencing. QL and WH modified the manuscript. WH and MS coordinated preparation of the manuscript. All authors read and approved the final manuscript.


This work was funded by the National Natural Science Foundation of China (Nos. 81271867 and 31725025), International Science and Technology Cooperation Program of China (No. 2014DFA31130).

Availability of data and materials

All raw data of 48 samples small RNA sequencing were submitted to GenBank with the project accession number of PRJNA512425. Nucleotide sequences have been deposited in the Sequence Read Archive (SRA) of NCBI. The dataset supporting the conclusions of this article is included within the article and in additional files.

Ethics approval and consent to participate

All laboratory animal tests were approved by the Ethics Committee of the National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention in Shanghai, China (Ref. No.: 20100525-1). The use of mice in these experiments was conducted in adherence to the guidelines for the Care and Use of Laboratory Animals of the Ministry of Science and Technology of People’s Republic of China ((2006)398).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

College of Veterinary Medicine, Northeast Agricultural University, Harbin, 150030, China
State Key Laboratory of Genetic Engineering, Ministry of Education Key Laboratory of Contemporary Anthropology, Collaborative Innovation Center for Genetics and Development, School of Life Sciences, Fudan University, Shanghai, 200438, China
Joint Research Laboratory of Genetics and Ecology on Parasite-host Interaction, Chinese Center for Disease Control and Prevention & Fudan University, WHO Collaborating Centre for Tropical Diseases, National Center for International Research on Tropical Diseases, Key Laboratory of Parasite and Vector Biology Ministry of Health, National Institute of Parasitic Diseases, Shanghai, 200025, China
Heilongjiang Key Laboratory for Zoonosis, Harbin, 150030, China


  1. You H, Cai P, Tebeje BM, Li Y, McManus DP (2018) Schistosome vaccines for domestic animals. Trop Med Infect Dis 3:E68View ArticleGoogle Scholar
  2. Devkota R, Brant SV, Thapa S, Loker ES (2014) Two avian schistosome cercariae from Nepal, including a Macrobilharzia-like species from Indoplanorbis exustus. Parasitol Int 63:374–380View ArticleGoogle Scholar
  3. Charlier J, van der Voort M, Kenyon F, Skuce P, Vercruysse J (2014) Chasing helminths and their economic impact on farmed ruminants. Trends Parasitol 30:361–367View ArticleGoogle Scholar
  4. Webster BL, Diaw OT, Seye MM, Webster JP, Rollinson D (2013) Introgressive hybridization of Schistosoma haematobium group species in Senegal: species barrier break down between ruminant and human schistosomes. PLoS Negl Trop Dis 7:e2110View ArticleGoogle Scholar
  5. Colley DG, Bustinduy AL, Secor WE, King CH (2014) Human schistosomiasis. Lancet 383:2253–2264View ArticleGoogle Scholar
  6. WHO (2016) Schistosomiasis: number of people treated worldwide in 2014. Wkly Epidemiol Rec 91:53–60Google Scholar
  7. Basch PF, Basch N (1984) Intergeneric reproductive stimulation and parthenogenesis in Schistosoma mansoni. Parasitology 89:369–376View ArticleGoogle Scholar
  8. Gupta BC, Basch PF (1987) The role of Schistosoma mansoni males in feeding and development of female worms. J Parasitol 73:481–486View ArticleGoogle Scholar
  9. Lu Z, Sessler F, Holroyd N, Hahnel S, Quack T, Berriman M, Grevelding CG (2016) Schistosome sex matters: a deep view into gonad-specific and pairing-dependent transcriptomes reveals a complex gender interplay. Sci Rep 6:31150View ArticleGoogle Scholar
  10. Davis A (2003) Schistosomiasis. In: Cook GC, Zumla AI (eds) Manson’s tropical diseases, 21st edn. WB Saunders, London, pp 1431–1469Google Scholar
  11. Lytle JR, Yario TA, Steitz JA (2007) Target mRNAs are repressed as efficiently by microRNA-binding sites in the 5′ UTR as in the 3′ UTR. Proc Natl Acad Sci USA 104:9667–9672View ArticleGoogle Scholar
  12. Bartel DP (2009) MicroRNAs: target recognition and regulatory functions. Cell 136:215–233View ArticleGoogle Scholar
  13. Bushati N, Cohen SM (2007) microRNA functions. Annu Rev Cell Dev Biol 23:175–205View ArticleGoogle Scholar
  14. Huang J, Hao P, Chen H, Hu W, Yan Q, Liu F, Han ZG (2009) Genome-wide identification of Schistosoma japonicum microRNAs using a deep-sequencing approach. PLoS One 4:e8206View ArticleGoogle Scholar
  15. Cai P, Piao X, Hao L, Liu S, Hou N, Wang H, Chen Q (2013) A deep analysis of the small non-coding RNA population in Schistosoma japonicum eggs. PLoS One 8:e64003View ArticleGoogle Scholar
  16. Cheng G, Luo R, Hu C, Cao J, Jin Y (2013) Deep sequencing-based identification of pathogen-specific microRNAs in the plasma of rabbits infected with Schistosoma japonicum. Parasitology 140:1751–1761View ArticleGoogle Scholar
  17. Sun J, Wang S, Li C, Ren Y, Wang J (2014) Novel expression profiles of microRNAs suggest that specific miRNAs regulate gene expression for the sexual maturation of female Schistosoma japonicum after pairing. Parasit Vectors 7:177View ArticleGoogle Scholar
  18. Zhu L, Zhao J, Wang J, Hu C, Peng J, Luo R, Zhou C, Liu J, Lin J, Jin Y, Davis RE, Cheng G (2016) MicroRNAs are involved in the regulation of ovary development in the pathogenic blood fluke Schistosoma japonicum. PLoS Pathog 12:e1005423View ArticleGoogle Scholar
  19. de Souza Gomes M, Muniyappa MK, Carvalho SG, Guerra-Sa R, Spillane C (2011) Genome-wide identification of novel microRNAs and their target genes in the human parasite Schistosoma mansoni. Genomics 98:96–111View ArticleGoogle Scholar
  20. Oliveira KC, Carvalho ML, Maracaja-Coutinho V, Kitajima JP, Verjovski-Almeida S (2011) Non-coding RNAs in schistosomes: an unexplored world. An Acad Bras Cienc 83:673–694View ArticleGoogle Scholar
  21. Marco A, Kozomara A, Hui JH, Emery AM, Rollinson D, Griffiths-Jones S, Ronshaugen M (2013) Sex-biased expression of microRNAs in Schistosoma mansoni. PLoS Negl Trop Dis 7:e2402View ArticleGoogle Scholar
  22. Hao L, Cai P, Jiang N, Wang H, Chen Q (2010) Identification and characterization of microRNAs and endogenous siRNAs in Schistosoma japonicum. BMC Genomics 11:55View ArticleGoogle Scholar
  23. Cai P, Hou N, Piao X, Liu S, Liu H, Yang F, Wang J, Jin Q, Wang H, Chen Q (2011) Profiles of small non-coding RNAs in Schistosoma japonicum during development. PLoS Negl Trop Dis 5:e1256View ArticleGoogle Scholar
  24. Protasio AV, van Dongen S, Collins J, Quintais L, Ribeiro DM, Sessler F, Hunt M, Rinaldi G, Collins JJ, Enright AJ, Berriman M (2017) MiR-277/4989 regulate transcriptional landscape during juvenile to adult transition in the parasitic helminth Schistosoma mansoni. PLoS Negl Trop Dis 11:e0005559View ArticleGoogle Scholar
  25. Wang J, Yu Y, Shen H, Qing T, Zheng Y, Li Q, Mo X, Wang S, Li N, Chai R, Xu B, Liu M, Brindley PJ, McManus DP, Feng Z, Shi L, Hu W (2017) Dynamic transcriptomes identify biogenic amines and insect-like hormonal regulation for mediating reproduction in Schistosoma japonicum. Nat Commun 8:14693View ArticleGoogle Scholar
  26. Chen X, Gao C, Li H, Huang L, Sun Q, Dong Y, Tian C, Gao S, Dong H, Guan D, Hu X, Zhao S, Li L, Zhu L, Yan Q, Zhang J, Zen K, Zhang CY (2010) Identification and characterization of microRNAs in raw milk during different periods of lactation, commercial fluid, and powdered milk products. Cell Res 20:1128–1137View ArticleGoogle Scholar
  27. McGinn J, Czech B (2014) Small RNA library construction for high-throughput sequencing. Methods Mol Biol 1093:195–208View ArticleGoogle Scholar
  28. Ow MC, Lau NC, Hall SE (2014) Small RNA library cloning procedure for deep sequencing of specific endogenous siRNA classes in Caenorhabditis elegans. Methods Mol Biol 1173:59–70View ArticleGoogle Scholar
  29. Brown J, Pirrung M, McCue LA (2017) FQC Dashboard: integrates FastQC results into a web-based, interactive, and extensible FASTQ quality control tool. Bioinformatics 33(9):3137–3139View ArticleGoogle Scholar
  30. Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120View ArticleGoogle Scholar
  31. Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N (2012) miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res 40:37–52View ArticleGoogle Scholar
  32. The Schistosoma japonicum Genome Sequencing, Functional Analysis Consortium (2009) The Schistosoma japonicum genome reveals features of host–parasite interplay. Nature 460:345–351View ArticleGoogle Scholar
  33. Kozomara A, Griffiths-Jones S (2014) miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res 42:D68–D73View ArticleGoogle Scholar
  34. da Silveira JC, Veeramachaneni DN, Winger QA, Carnevale EM, Bouma GJ (2012) Cell-secreted vesicles in equine ovarian follicular fluid contain miRNAs and proteins: a possible new form of cell communication within the ovarian follicle. Biol Reprod 86:71View ArticleGoogle Scholar
  35. Wang Z, Xue X, Sun J, Luo R, Xu X, Jiang Y, Zhang Q, Pan W (2010) An “in-depth” description of the small non-coding RNA population of Schistosoma japonicum schistosomulum. PLoS Negl Trop Dis 4:e596View ArticleGoogle Scholar
  36. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods 25:402–408View ArticleGoogle Scholar
  37. Kruger J, Rehmsmeier M (2006) RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res 34:W451–W454View ArticleGoogle Scholar
  38. Peterson SM, Thompson JA, Ufkin ML, Sathyanarayana P, Liaw L, Congdon CB (2014) Common features of microRNA target prediction tools. Front Genet 5:23View ArticleGoogle Scholar
  39. Hausser J, Syed AP, Bilen B, Zavolan M (2013) Analysis of CDS-located miRNA target sites suggests that they can effectively inhibit translation. Genome Res 23:604–615View ArticleGoogle Scholar
  40. Guo ZW, Xie C, Yang JR, Li JH, Yang JH, Zheng L (2015) MtiBase: a database for decoding microRNA target sites located within CDS and 5′UTR regions from CLIP-Seq and expression profile datasets. Database. View ArticlePubMedPubMed CentralGoogle Scholar
  41. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M (2005) Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21:3674–3676View ArticleGoogle Scholar
  42. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T (2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13:2498–2504View ArticleGoogle Scholar
  43. Xue X, Sun J, Zhang Q, Wang Z, Huang Y, Pan W (2008) Identification and characterization of novel microRNAs from Schistosoma japonicum. PLoS One 3:e4034View ArticleGoogle Scholar
  44. Nollen PM, Floyd RD, Kolzow RG, Deters DL (1976) The timing of reproductive cell development and movement in Schistosoma mansoni, S. japonicum, and S. haematobium, using techniques of autoradiography and transplantation. J Parasitol 62:227–231View ArticleGoogle Scholar
  45. Garcia-Estrada C, Ullan RV, Velasco-Conde T, Godio RP, Teijeira F, Vaca I, Feltrer R, Kosalkova K, Mauriz E, Martin JF (2008) Post-translational enzyme modification by the phosphopantetheinyl transferase is required for lysine and penicillin biosynthesis but not for roquefortine or fatty acid formation in Penicillium chrysogenum. Biochem J 415:317–324View ArticleGoogle Scholar
  46. Fung AW, Fahlman RP (2015) The molecular basis for the post-translational addition of amino acids by L/F transferase in the N-end rule pathway. Curr Protein Pept Sci 16:163–180View ArticleGoogle Scholar
  47. Fabrini R, Bocedi A, Camerini S, Fusetti M, Ottaviani F, Passali FM, Topazio D, Iavarone F, Francia I, Castagnola M, Ricci G (2014) Inactivation of human salivary glutathione transferase P1-1 by hypothiocyanite: a post-translational control system in search of a role. PLoS ONE 9:e112797View ArticleGoogle Scholar
  48. Puglisi I, Lo Cicero L, Lo Piero AR (2013) The glutathione S-transferase gene superfamily: an in silico approach to study the post translational regulation. Biodegradation 24:471–485View ArticleGoogle Scholar
  49. Okada Y, Scott G, Ray MK, Mishina Y, Zhang Y (2007) Histone demethylase JHDM2A is critical for Tnp1 and Prm1 transcription and spermatogenesis. Nature 450:119–123View ArticleGoogle Scholar
  50. O’Carroll D, Scherthan H, Peters AH, Opravil S, Haynes AR, Laible G, Rea S, Schmid M, Lebersorger A, Jerratsch M, Sattler L, Mattei MG, Denny P, Brown SD, Schweizer D, Jenuwein T (2000) Isolation and characterization of Suv39h2, a second histone H3 methyltransferase gene that displays testis-specific expression. Mol Cell Biol 20:9423–9433View ArticleGoogle Scholar
  51. Tachibana M, Sugimoto K, Nozaki M, Ueda J, Ohta T, Ohki M, Fukuda M, Takeda N, Niida H, Kato H, Shinkai Y (2002) G9a histone methyltransferase plays a dominant role in euchromatic histone H3 lysine 9 methylation and is essential for early embryogenesis. Genes Dev 16:1779–1791View ArticleGoogle Scholar


© The Author(s) 2019