Skip to main content

Differential expression of microRNA between normally developed and underdeveloped female worms of Schistosoma japonicum


Eggs produced by bisexual infected mature female worms (MF) of Schistosoma japonicum are important in the transmission of the parasite and responsible for the pathogenesis of schistosomiasis. The single-sex infected female worms (SF) cannot mature and do not produce normal eggs; also they do not induce severe damage to the host. In this study, the microRNA (miRNA) expression profiles of 25d MF and 25d SF were investigated through Solexa deep-sequencing technology to explore the developmental mechanisms of schistosome female worms. There were 36 differentially expressed miRNA, 20 up-regulated and 16 down-regulated found in MF/SF worms, including some development related miRNA such as bantam (ban), let-7, miR-124, miR-8, miR-1, miR-7. There were 166 target genes of up-regulated miRNA and 201 target genes of down-regulated miRNA after comparing the target gene prediction software results with RNA-Seq transcriptome results. Analysis of the target genes shows that different ones are involved in MF and SF worms in Gene Ontology terms, with a similar situation in KEGG. This observation indicates that different genes regulated by differentially expressed miRNA take part in MF and SF and lead to differential sexual status. This means that the sexual status of female worms is regulated by miRNA.


Schistosomiasis is an important helminth infection that mainly occurs in developing countries. It is endemic in 74 countries, with roughly 120 million individuals being severely affected [1]. No successful vaccine is currently available for this disease. Praziquantel is the only and most extensively used drug for treatment, but it is ineffective against young worms [2] and there are concerns that drug resistance to praziquantel may be developing [3]. Consequently, there is an urgent need to identify novel drug targets and/or develop alternative therapeutic strategies.

In schistosomiasis, the severely pathological injuries to the final host are caused by the amounts of eggs produced by mature normally developed female worms (MF). These eggs are also responsible for the spread of the disease [4], but the underdeveloped female worms (SF) which come from single-sex infection cannot develop normally or produce normal eggs. So these single-sex infected worms (female or male) will not seriously damage the host, and will not induce the prevalence of schistosomiasis. Based on the important role of eggs produced by mature females, the developmental mechanisms of female worms is the research target for many schistosome researchers.

MicroRNA (miRNA) are a class of small non-coding RNA (18–25 nucleotides in length) that are generated from endogenous transcripts [5]. They function as guide molecules for post-transcriptional gene regulation by base-pairing to their messenger RNA (mRNA) targets [6], they bind to the target and induce mRNA translation repression and/or exonucleolytic decay. Through these actions, miRNA regulate gene expression during development, differentiation, proliferation, death, and metabolism in many organisms [7]. To date, several groups have carried out studies to identify miRNA in schistosomes to provide some insight into the role of miRNA in schistosome development [8,9,10].

In this study, the expression profiles of miRNA of 25d MF and SF Schistosoma japonicum (S. japonicum) were investigated. Based on the analysis of the differentially expressed miRNA and their predicted target genes, the correlation between miRNA and the mature status of female worms was studied.

Materials and methods

Ethics statement

All animal experiments were conducted in accordance with the guidelines of the Committee for the Care and Use of Laboratory Animals of Shanghai Veterinary Research Institute, Chinese Academy of Agricultural Sciences (permit no. SYXK-2016-0010). The study protocol was approved by the Ethics and Animal Welfare Committee of the Shanghai Veterinary Research Institute, Chinese Academy of Agricultural Sciences.

Host animals and parasites

BALB/C male mice, 4- to 6-week-old, were used for worm collection. S. japonicum (Chinese strain) worms were maintained in mice and Oncomelania hupensis hupensis snails at the Shanghai Veterinary Research Institute, Chinese Academy of Agricultural Sciences (SHVRI, CAAS). Cercariae of S. japonicum were shed from infected snails after exposure to artificial light for 5 h at 25–28 °C. The 25d MF were collected from bi-sexual infected mice at 25 days post-infection by perfusion. Worms were washed five times in phosphate buffered saline (PBS, pH 7.4) at 4 °C to remove excreted and secreted proteins and other interfering substances and were then incubated at 4 °C for 30 min. The paired worms spontaneously separated into females and males, then the female worms were collected. The 25d SF were obtained after perfusion from mice infected by female cercariae. The sexes of cercariae shed from individual snails were identified by polymerase chain reaction (PCR) based on the patent: CN 101597645. Then the female worms were washed five times in PBS at 4 °C to remove excreted and secreted proteins and other interfering substances, then collected.

RNA extraction, small RNA library construction and sequencing

Total RNA was extracted using TRIzol reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) according to the manufacturer’s instructions. The quality of RNA was measured using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The RNA samples were stored at -80 °C.

The construction of small RNA libraries was carried out using the TruSeq Small RNA Sample Prep Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s protocol. Briefly, total RNA was subjected to polyacrylamide gel electrophoresis (PAGE) and then the 18–30 nt size range of RNA was enriched for purification. The purified small RNA was subjected to sequencing with an Illumina Genome Analyzer (Illumina) according to the manufacturer’s instructions. Briefly, proprietary adapters were ligated to the 5ʹ and 3ʹ termini of the small RNA, and the ligated small RNA were then used as templates for cDNA synthesis, and further converted into 62–75 nt single-stranded cDNA. The cDNA were amplified using illumina’s 3′ adaptor reverse primer and 5′ adaptor forward primers. The purified PCR products were sequenced using Solexa’s proprietary sequencing-by-synthesis method. DNA sequencing was performed with an Illumina Genome Analyzer at the BGI (Beijing Genomics Institute, Shenzhen, China).

Data processing and mapping sequence reads to the reference genome

All raw datasets produced by deep sequencing from each library (25d MF, 25d SF) were subjected to further analysis. Clean reads were obtained after removal of low quality reads, insert null reads, adaptor null reads, 5′ adaptor contaminants, reads with poly (A) tails and reads shorter than 18 nt. Adapter sequences were then trimmed from both ends of the clean reads. The small RNA with sizes ranging from 18–30 nt were used for further analyses. The clean reads in each library were mapped onto the draft S.japonicum genome sequences ( and ( using the Short Oligonucleotide Alignment Program (SOAP). The small RNA that matched with known rRNA, tRNA, small nuclear RNA (snRNA), and small nucleolar RNA (snoRNA) deposited in the Rfam database ( and NCBI GenBank ( were excluded. Then, we aligned the predicted miRNA to known miRNA or homologous miRNA in the miRBase (version 21, to identify conserved miRNA Other small RNA were further analyzed to predict new miRNA by exploring the secondary structure, the Dicer cleavage site and minimum free energy of the unannotated small RNA tags that could be mapped to the genome using MIREAP (

Bioinformatics analysis of S. japonicum small RNA

Small RNA tags were aligned to the miRNA precursor/mature miRNA of corresponding species in miRBase20. Detailed information of the alignment, including the structure of known miRNA precursors, and the length and count of tags from the sample, among others, were collected. To make every unique small RNA mapped to only one annotation, we followed the following priority Ruler RNA etc.: rRNA etc. (in which Genbank > Rfam) > known miRNA > repeat > exon > intron3. The total rRNA proportion should be less than 40%. The expression of the miRNA in two samples were visualized by plotting the Log2-ratio figure and the scatter plot. The procedures adopted were as follows: (1) the expression of miRNA in two samples (25d MF and 25d SF) were normalized to get the expression in transcripts per million (TPM) with the normalization formula (normalized expression = actual miRNA count/total count of clean reads × 1 000 000); and (2) the fold-change and P-value were calculated from the normalized expression, and log2-ratio and scatter plots were generated using the fold-change formula (fold change = log2 (25d MF/25d SF)).

MiRNA quantification by quantitative RT-PCR analysis

Differentially expressed miRNA were examined using qPCR with SYBR green. The first cDNA strand was synthesized from 0.2 μg total RNA using PrimeScript® RT reagent Kit (Perfect Real Time, TaKaRa, Otsu, Japan, code: DRR037A). Stem-loop qRT-PCR was performed to quantify the sex-biased expressed miRNA. Stem-loop RT primers (Table 1) were used to reverse-transcribe mature miRNA to cDNA. The RNA templates for the qPCR were performed on the same samples used for Solexa deep sequencing. The 20 μL reaction RT system contained 2 μL of total RNA (0.2 μg), 0.4 μL (10 μM) of each individual stem-loop RT primer, 0.4 μL ROX Reference Dye II (50 ×), 10 μL SYBR Premix Ex Taq (2 ×), and 6.8 μL Easy Dilution. The RT primers of the miRNA, the forward primers of miRNA and the common reverse primer were designed as shown in Table 1. The cDNA were synthesized by incubation for 30 min at 16 °C, 45 min at 37 °C, and 15 s at 85 °C. The products were amplified using SYBR® Premix ExTaq (TliRNaseH Plus, TaKaRa, code: RR420A) in an ABI Prism 7500 sequence detection system (Applied Biosystems, Waltham, MA, USA) according to the manufacturer’s instruction. Our system contained 2 μL of cDNA from the RT reaction product, 5 μL of 2 × SYBR Premix ExTaq (TaKaRa), and 0.4 μL of 10 μM forward and reverse primers. The NADPH gene was used in the qRT-PCR reaction as an endogenous reference [8]. Relative levels of gene expression were calculated using the 2 − Ct method. All assays were performed in triplicate.

Table 1 Primers of reverse transcription and RT-PCR.

Prediction and bioinformatics analysis of differentially expressed miRNA target genes

The regulatory effect of differentially expressed miRNA in 25d MF and 25d SF were predicted by analyzing the functions of their predicted target genes and by comparing the proportion of their predicted target genes in differentially expressed genes. In this study, RNA hybrid ( [11] were used to predict the target genes of differentially expressed miRNA. The predicted target genes were filtered with minimum free energy (MFE <  = -20 kcal/mol). These target genes were re-filtered by comparison with the RNA-Seq (Quantification) transcriptome results of SF and MF. The worms were the same sets used in miRNA deep sequencing (Data not published); the genes that were not consistent in the gene expression profile and miRNA target gene profile were excluded. Then those corresponding to target genes were mapped to the terms in the Gene Ontology (GO) database, to identify significantly enriched GO terms in DEG (differentially expressed genes) compared with the genome background. Pathway enrichment analysis was used to identify significantly enriched metabolic pathways or signal transduction pathways in DEG compared with the whole genome background. Detailed pathway information was determined with the KEGG database.


Host animals and parasites

Cercariae were collected from individual snails and their sex was identified by PCR. There is a 153 bp specific band from female S. japonicum cercariae, and no band from male S. japonicum cercariae (Figure 1). For those snails that contain two sex cercariae, the SF female worms were guaranteed by sacrificing the mice that were infected by cercariae from the same snail when the worms were collected.

Figure 1

Identification of the sex of cercariae shed from different snails by PCR. Lane M, DNA Marker, Lanes 1–6, Cercariae shed from different snails.

Preliminary analysis of raw reads of miRNA in S. japonicum 25d MF and 25d SF

The cDNA libraries derived from 18–30 nt long RNA isolated from 25d MF and 25d SF were constructed and sequenced using an Illumina (Solexa) DNA sequencer, yielding a total of 9 358 376 and 9 314 275 clean reads corresponding to 97.74% and 97.29% of high-quality sequence for S. japonicum 25d MF and 25d SF, respectively (Table 2). The length distribution of the miRNA from the two libraries mainly ranged from 20–24 nt (Figure 2).

Table 2 Proportion of clean reads and contaminant reads that were removed.
Figure 2

Length, distribution and abundance of small miRNA tags of 25d MF and 25d SF.

Differential expression profiles of miRNA in S. japonicum 25d MF and 25d SF

All the differentially expressed miRNA profile data for S. japonicum 25d MF and 25d SF were normalized to TPM as previously described [11,12,13]. As shown in Table 3, statistical analysis revealed that 36 S. japonicum miRNA were differentially expressed in the 25d SF and 25d MF. 20 miRNA were up-regulated (fold-change (log2 (25d MF/25d SF)) > 1), and 16 miRNA were down-regulated (fold-change (log2 (25d MF/25d SF)) < -1) in 25d MF.

Table 3 Differentially expressed miRNA of 25d SF and 25d MF.

Some high-abundance miRNA, such as ban, let-7, miR-1, miR-7-5p and so on, were significantly differentially expressed between 25d MF and 25d SF. Some high-abundance miRNA, however, such as miR-125b, miR-190-5p, miR-2162-3p, miR-277, miR-2a-3p, miR-71a and miR-71b-5p, were not significantly differentially expressed between the two samples (Figure 3). Meanwhile, some low-abundance miRNA, such as miR-3490, miR-3479-3p, miR-3505, miR-2c-5p (Table 3), were significantly differentially expressed too, and there are also some low-abundance miRNA that were not differentially expressed (Data not shown). Some miRNA that were functionally characterized as developmental-related in other organisms were differentially expressed in these two samples, such as ban, let-7, miR-124, miR-8, miR-1, and miR-7.

Figure 3

Comparison of the miRNA profiles from 25d SF and 25d MF. **Fold-change (log2) > 1 or fold-change (log2) < -1, and P value < 0.01.

Confirmation of differentially expressed miRNA by quantitative RT-PCR analysis

To confirm the differential expression of miRNA in 25d MF and 25d SF, let-7, miR-1, miR-10-5p, miR-124-3p and miR-3479-3p were selected for quantitative RT-PCR analysis. The results of the RT-PCR showed that let-7, miR-1, miR-10-5p, miR-124-3p and miR-3479-3p were expressed higher in 25d SF than in 25d MF (Figure 4). These results were consistent with the Solexa analysis.

Figure 4

Result of qPCR confirmation of miRNA expression.

Predictions and bioinformatics analysis of the predicted target genes for differentially expressed miRNA

To analyze the effect of the differential expression of miRNA on 25d MF and 25d SF, the target genes of the differentially expressed miRNA between the two samples were predicted and analyzed. The target genes of up- or down-regulated miRNA were predicted and analyzed separately. Based on the target gene prediction software, 4728 and 3906 target genes of up- and down-regulated miRNA were acquired respectively. After comparison with the RNA-Seq transcriptome results, there were 166 target genes of up-regulated miRNA and 201 target genes of down-regulated miRNA left respectively. The bioinformatics analysis for target genes of either up- or down-regulated miRNA was applied separately. Not all the target genes can be applied to GO terms, for example, only 41 out of 166 target genes of up-regulated miRNA and 76 out of 201 target genes of down-regulated miRNA were logged in the GO terms of molecular function (Additional file 1). The target genes of up- or down-regulated miRNA were mainly concentrated on the same GO terms (Additional file 2 and Additional file 3). For molecular-function, they were both concentrated on binding, catalytic activity, molecular transducer activity, transporter activity (Figure 5). But the specific target genes of each GO term were different for up- or down-regulated miRNA and there were also target genes of up- and down-regulated miRNA concentrated on different GO terms (Figure 5). For biological process and cellular component GO terms, the same situation was found (Additional file 4 and Additional file 5).

Figure 5

GO terms of molecular-function in the target genes of up- and down-regulated miRNA.

Some target genes of up- or down-regulated miRNA were selected to analyze their functions. The results show that target genes were not only regulated by only one miRNA, such as the target gene Sjc_0005980, but that a target gene can be involved in more than one GO term (Additional file 6 and Additional file 7). This means they participate in multiple important biological processes.

Pathway enrichment analysis also testified to this situation. KEGG results show that some target genes of up- and down-regulated miRNA take part in the same pathway, but the individual target genes in up- and down-regulated miRNA are different. In addition also some target genes of up- or down-regulated miRNA are involved in various pathways, however, only target genes of up-regulated miRNA are involved in these pathways, such as starch and sucrose metabolism, carbohydrate digestion and absorption, TGF-beta signaling pathway, chemokine signaling pathway, and so on. For the pathways, protein processing in the endoplasmic reticulum, N-glycan biosynthesis, ribosome biogenesis in eukaryotes, etc., only target genes of down-regulated miRNA take part (Table 4).

Table 4 Comparison of predicted target genes participating in regulating metabolic processes between up-regulated miRNAs and down-regulated miRNAs in 25d MF.


The MF worms of S. japonicum can develop to sexual maturation and maintain this mature state by pairing with male worms [14]. The SF worms were considerably different in size and structure of inner tissues or systems from those MF worms [15]. Although the long term infected SF can lay very few eggs, these non-viable eggs do not contain developing or live miracidia so there are no soluble miracidial antigens that can induce the host’s immune response, therefore the pathological reaction in infected final host and the schistosomiasis spreading in the environment will be relieved or stopped [16]. During the development of S. japonicum, males and females begin to pair at 15 or 16 d post-infection, and females begin to lay eggs from 24 d post-infection [14]. So S. japonicum will develop totally mature worms at 25 d post-infection. In this study, the miRNA profiles of 25d MF and 25d SF were compared, some differentially expressed miRNA in the two samples were found, and some clues about the correlation between differentially expressed miRNA and the mature status of schistosome female worms appeared by analyzing the differentially expressed miRNA and their target genes.

The adult female and male worms have obvious gender dimorphism, which can be identified by naked eyes, but the sex of miracidium, sporocyst and cercariae cannot be identified in this way, as they do not have obvious gender dimorphism. One snail can be infected by several miracidia, but in the majority of snails (65.7%), only one miracidium can develop into a mother sporocyst [17]. The proportion of single- or mixed-sex infections of parasites within snail individuals in the field has been reported for S. japonicum, most cercariae shed from an individual snail were single-sex [18]. Based on this situation, in this study, the sex of cercariae overflow from an individual snail was confirmed by PCR and infected mice were sacrificed, to ensure that SF worms were recovered.

There are significant morphological differences between normal developed and underdeveloped female S. japonicum. In our previous research, the underdeveloped female S. japonicum was significantly shorter and thinner than normal developed, and the mature ootid in the indistinct ovary of the underdeveloped female could not be observed [19]. These distinct morphological changes may be associated with miRNA regulation, in this study; there are 20 up-regulated and 16 down-regulated miRNA in 25d MF. The differentially expressed miRNA: miR-1, miR-124, miR-8, miR-7, let-7, and ban, have been functionally characterized as developmental regulators in other organisms [20].

Bantam is a miRNA that is functionally well characterized in Drosophila. In Drosophila, it is a developmentally regulated miRNA that controls cell proliferation and regulates the proapoptotic gene hid [21]. Bantam was also identified with genes promoting organ growth and systemic body growth by connecting insulin signaling, and ecdysone production [22]. Boulan et al. found that lower ban activity induced growth cessation of late ecdysone-producing cells that led to the cells entry into a metamorphosis stage, whereas ban’s high activity promotes the larvae’s systemic growth in Drosophila [22]. In our study, ban was also up-regulated in developed females and down-regulated in underdeveloped females (Table 3). These phenomena are also observed in Caenorhabditis elegans, the deletion of ban orthologs in Caenorhabditis elegans cause severe growth defects and prevent entry into the dauer state [23].

In this study, two miRNA, miR-7 and miR-1, which have been characterized as tumor suppressors in many kinds of human cancers, were up-regulated in 25d SF. They were found to play a key role in cell proliferation and survival either by being involved in the insulin-like growth factor 1 (IGF1) receptor signaling pathway or epidermal growth factor (EGF) receptor signaling pathway. EGF is an early growth response gene that is induced by mitogenic stimuli, including serum and purified factors such as IGF, EGF, and TGF-beta [24]. Let-7 and miR-124 were also up-regulated in 25d SF in our study. The Let-7 miRNA family has an important function in cell differentiation and cancer [25]. Worringer et al. suggested that let-7 based pathways are an inhibitory influence on the reprogramming process through a regulatory pathway involving prodifferentiation factors, including EGF receptor 1 [26]. MiR-124 also plays a key role in cell differentiation, especially in the nervous system [27]; it was also thought to be a growth suppressor [28]. MiR-8, recognized as a critical mediator for several well characterized developmental regulatory networks such as Wnt, Notch, and TGF-beta in Drosophila or mammalian cells [29, 30], is up-regulated in 25d SF too. In addition, the Wnt, Notch, and TGF-beta pathways are shown to be involved in the parasite development, oogenesis and embryogenesis in schistosoma mansoni [31,32,33,34]. The up-regulation of these differentiation-associated miRNA in 25d SF may be related to the underdevelopmental status of SF. Sexual immaturity may mean more potential for differentiation when a single female pairs with her mate.

The bioinformatics analysis of the predicted target genes for differentially expressed miRNA shows that the target genes distributed in different GO terms, which means these differential functional genes, take part in many important biological processes. The KEGG results show that more target genes in 25d MF are involved in the regulation of actin cytoskeleton, focal adhesion, protein digestion and absorption, tight junction, and calcium signaling pathway metabolic processes whereas more target genes in 25d SF are involved in lysosomes, oocyte meiosis, long-term potentiation and salivary secretion pathways. Only target genes of 25d MF were detected in starch and sucrose metabolism, galactose metabolism, carbohydrate digestion and absorption, axon guidance, TGF-beta signaling pathway, and chemokine signaling pathways whereas only target genes in 25d SF were found in pathways involved in protein processing in the endoplasmic reticulum, N-glycan biosynthesis, ribosome biogenesis eukaryotes, gap junction, dopaminergic synapse, spliceosome, RNA transport, GnRH signaling pathway, and long-term depression. All these are associated with the developmental stages of female worms.

Some receptors in the Wnt, Notch, IGF, EGF, TGF-beta signaling pathways which interact with miRNA have been characterized in both S. mansoni and S. japonicum [35, 36], these signaling pathways have also been implicated in critical roles in schistosome development by the S. japonicum genome study [37], but direct evidence for miRNA regulated signaling pathways remains to be experimentally clarified. This study shows that some miRNA could regulate the sexual development of female worms and that it analyzed target genes. This will shed some clues on the understanding of the development mechanism of schistosome female worms.

Availability of data and materials

The microRNA sequence data were deposited in the SRA database at NCBI with SRA, Project accession number: PRJNA648823.


S. japonicum :

Schistosoma japonicum


ribonucleic acid


the single-sex infected female worms


deoxyribonucleic acid


the bisexual infected mature female worms


polymerase chain reaction


polyacrylamide gel electrophoresis


quantitative real-time PCR


nicotinamide adenine dinucleotide phosphate


differentially expressed genes


short Oligonucleotide Alignment Program


transcripts per million


minimum free energy


gene Ontology


Kyoto Encyclopedia of Genes and Genomes


phosphate-bufered saline




insulin-like growth factor


epidermal growth factor


transforming growth factor


  1. 1.

    Steinmann P, Keiser J, Bos R et al (2009) Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis 6:411–425

    Article  Google Scholar 

  2. 2.

    Doenhoff MJ, Pica ML (2006) Praziquantel for the treatment of schistosomiasis: its use for control in areas with endemic disease and prospects for drug resistance. Expert Rev Anti Infect Ther 4:199–210

    CAS  Article  Google Scholar 

  3. 3.

    Ismail M, Botros S, Metwally A et al (1999) Resistance to praziquantel: direct evidence from Schistosoma mansoni isolated from Egyptian villagers. Am J Trop Med Hyg 60:932–935

    CAS  Article  Google Scholar 

  4. 4.

    Andrade ZA (2009) Schistosomiasis and liver fibrosis. Parasite Immunol 31:656–663

    CAS  Article  Google Scholar 

  5. 5.

    Kim VN (2005) MicroRNA biogenesis: coordinated cropping and dicing. Nat Rev Mol Cell Biol 6:376–385

    CAS  Article  Google Scholar 

  6. 6.

    Bartel DP (2004) MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116:281–297

    CAS  Article  Google Scholar 

  7. 7.

    Bushati N, Cohen SM (2007) microRNA functions. Annu Rev Cell Dev Biol 23:175–205

    CAS  Article  Google Scholar 

  8. 8.

    Cheng G, Luo R, Hu C et al (2013) Deep sequencing-based identification of pathogen-specific microRNAs in the plasma of rabbits infected with Schistosoma japonicum. Parasitology 140:1751–1761

    CAS  Article  Google Scholar 

  9. 9.

    Sun J, Wang S, Li C et al (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:177

    Article  Google Scholar 

  10. 10.

    Han H, Peng J, Hong Y et al (2015) Comparative characterization of microRNAs in Schistosoma japonicum schistosomula from Wistar rats and BALB/c mice. Parasitol Res 114:2639–2647

    Article  Google Scholar 

  11. 11.

    Huang J, Hao P, Zhang YL et al (2007) Discovering multiple transcripts of human hepatocytes using massively parallel signature sequencing (MPSS). BMC Genomics 8:207

    Article  Google Scholar 

  12. 12.

    Meyers BC, Tej SS, Vu TH et al (2004) The use of MPSS for whole-genome transcriptional analysis in Arabidopsis. Genome Res 14:1641–1653

    CAS  Article  Google Scholar 

  13. 13.

    Brenner S, Johnson M, Bridgham J et al (2000) Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays. Nat Biotechnol 18:630–634

    CAS  Article  Google Scholar 

  14. 14.

    He YX, Yang HZ (1980) Physiological studies on the post-cercarial development of Schistosoma japonicum. Chin J Dong Wu Xue Bao 26:32–41

    Google Scholar 

  15. 15.

    Han Y, Qiao HB, Cao YF et al (2015) Morphology and apoptosis of female Schistosoma japonicum from single-sex and mixed-sex infecteions. Chin J Anim Infect Dis 23:30–35

    Google Scholar 

  16. 16.

    Shaw MK (1987) Schistosoma mansoni: vitelline gland development in females from single sex infections. J Helminthol 61(3):253–259

    CAS  Article  Google Scholar 

  17. 17.

    Zhou SL (1958) A study on the development of the larval stage of Schistosoma japonicum in oncomelania hupensis. Chin J Wei Sheng Wu Xue Bao 6:110–126

    Google Scholar 

  18. 18.

    Shi HP, Lu DB, Shen L et al (2014) Single- or mixed-sex Schistosoma japonicum infections of intermediate host snails in hilly areas of Anhui, China. Parasitol Res 113:717–721

    Article  Google Scholar 

  19. 19.

    Li X, Qiao H, Qin F et al (2020) Comparative analysis of iTRAQ-based proteome profiles of Schistosoma japonicum female worms coming from single-sex infections and bisexual infections. J Proteomics 213:103597

    CAS  Article  Google Scholar 

  20. 20.

    Cheng G, Jin Y et al (2012) MicroRNAs: potentially important regulators for schistosome development and therapeutic targets against schistosomiasis. Parasitology 139:669–679

    CAS  Article  Google Scholar 

  21. 21.

    Brennecke J, Hipfner DR, Stark A et al (2003) Bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell 113:25–26

    CAS  Article  Google Scholar 

  22. 22.

    Boulan L, Martín D, Milán M (2013) Bantam miRNA promotes systemic growth by connecting insulin signaling and ecdysone production. Curr Biol 23:473–478

    CAS  Article  Google Scholar 

  23. 23.

    Alvarez-Saavedra E, Horvitz HR (2010) Many families of C.elegans microRNAs are not essential for development or viability. Curr Biol 20:367–373

    CAS  Article  Google Scholar 

  24. 24.

    Gashler A, Sukhatme VP (1995) Early growth response protein 1(Egr-1): prototype of a zinc-finger family of transcription factors. Prog Nucleic Acid Res Mol Biol 50:191–224

    CAS  Article  Google Scholar 

  25. 25.

    Boyerinas B, Park SM, Hau A et al (2010) The role of let-7 in cell differentiation and cancer. Endocr Relat Cancer 17:F19–36

    CAS  Article  Google Scholar 

  26. 26.

    Worringer KA, Rand TA, Hayashi Y et al (2014) The let-7/LIN-41 pathway regulates reprogramming to human induced pluripotent stem cells by controlling expression of prodifferentiation genes. Cell Stem Cell 14:1–13

    Article  Google Scholar 

  27. 27.

    Makeyev EV, Zhang J, Carrasco MA et al (2007) The MicroRNA miR-124 promotes neuronal differentiation by triggering brain-specific alternative pre-mRNA splicing. Mol Cell 27:435–448

    CAS  Article  Google Scholar 

  28. 28.

    Li KK, Pang JC, Ching AK et al (2009) miR-124 is frequently down-regulated in medulloblastoma and is a negative regulator of SLC16A1. Hum Pathol 40:1234–1243

    CAS  Article  Google Scholar 

  29. 29.

    Karres JS, Hilgers V, Carrera I et al (2007) The conserved microRNA miR-8 tunes atrophin levels to prevent neurodegeneration in Drosophila. Cell 131:136–145

    CAS  Article  Google Scholar 

  30. 30.

    Kennell JA, Gerin I, MacDougald OA et al (2008) The microRNA miR-8 is a conserved negative regulator of Wnt signaling. Proc Natl Acad Sci USA 105:15417–15422

    CAS  Article  Google Scholar 

  31. 31.

    Liang S, Knight M, Jolly ER (2013) Polyethyleneimine mediated DNA transfection in schistosome parasites and regulation of the WNT signaling pathway by a dominant-negative SmMef2. PLoS Negl Trop Dis 7:e2332

    CAS  Article  Google Scholar 

  32. 32.

    Magalhães LG, Morais ER, Machado CB (2016) Uncovering Notch pathway in the parasitic flatworm Schistosoma mansoni. Parasitol Res 115(10):3951–3961

    Article  Google Scholar 

  33. 33.

    Freitas TC, Jung E, Pearce EJ (2007) TGF-beta signaling controls embryo development in the parasitic flatworm Schistosoma mansoni. PLoS Pathog 3:e52

    Article  Google Scholar 

  34. 34.

    Buro C, Oliveira KC, Lu Z et al (2013) Transcriptome analyses of inhibitor treated schistosome females provide evidence for cooperating Src-kinase and TGFbeta receptor pathways controlling mitosis and eggshell formation. PLoS Pathog 9:e1003448

    CAS  Article  Google Scholar 

  35. 35.

    Khayath N, Vicogne J, Ahier A et al (2007) Diversification of the insulin receptor family in the helminth parasite Schistosoma mansoni. FEBS J 274:659–676

    CAS  Article  Google Scholar 

  36. 36.

    You H, Zhang W, Jones MK et al (2010) Cloning and characterisation of Schistosoma japonicum insulin receptors. PLoS One 5:e9868

    Article  Google Scholar 

  37. 37.

    Zhou Y, Zheng H, Chen Y et al (2009) The Schistosoma japonicum genome reveals features of host-parasite interplay. Nature 460:345–351

    CAS  Article  Google Scholar 

Download references


The authors thank International Science Editing for revising the manuscript and their valuable comments.


This work was supported by grants from the Chinese National Natural Science Foundation (no. 31672245) and the Natural Science Foundation of Shanghai (no. 16ZR1444100; no. 20ZR1469300).

Author information




YH: overall experimental design and execution of all laboratory experiments, writing of the manuscript, interpretation of all data; JF, YR, LW, HL: participated in the qRT-PCR analyses; JL: assisted with the data analysis; YJ: obtained financial support, conceived and designed the study and analyzed the data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yamei Jin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

The author confirms: the work described has not been published before; it is not under consideration for publication elsewhere; its publication has been approved by all co-authors; its publication has been approved by the responsible authorities at the institution where the work is carried out.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Han, Y., Feng, J., Ren, Y. et al. Differential expression of microRNA between normally developed and underdeveloped female worms of Schistosoma japonicum. Vet Res 51, 126 (2020).

Download citation


  • Schistosoma japonicum
  • microRNA
  • deep-sequencing
  • differential expression
  • predicted target genes