Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

Identification and characterization of microRNAs and endogenous siRNAs in Schistosoma japonicum

Lili Hao1, Pengfei Cai1, Ning Jiang2, Heng Wang1 and Qijun Chen123*

Author Affiliations

1 Laboratory of Parasitology, Institute of Pathogen Biology/Institute of Basic Medicine, Chinese Academy of Medical Sciences (CAMS), Beijing 100730, PR China

2 Key Laboratory of Zoonosis, Ministry of Education, Institute of Zoonosis, Jilin University, Changchun 130062, PR China

3 Swedish Institute for Infectious Disease Control and Karolinska Institutet, 171 82, Stockholm, Sweden

For all author emails, please log on.

BMC Genomics 2010, 11:55  doi:10.1186/1471-2164-11-55


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/11/55


Received:3 September 2009
Accepted:21 January 2010
Published:21 January 2010

© 2010 Hao et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Small endogenous non-coding RNAs (sncRNAs) such as small interfering RNA (siRNA), microRNA and other small RNA transcripts are derived from distinct loci in the genome and play critical roles in RNA-mediated gene silencing mechanisms in plants and metazoa. They are approximately 22 nucleotides long; regulate mRNA stability through perfect or imperfect match to the targets. The biological activities of sncRNAs have been related to many biological events, from resistance to microbe infections to cellular differentiation. The development of the zoonotic parasite Schistosoma japonicum parasite includes multiple steps of morphological alterations and biological differentiations, which provide a unique model for studies on the functions of small RNAs. Characterization of the genome-wide transcription of the sncRNAs will be a major step in understanding of the parasite biology. The objective of this study is to investigate the transcriptional profile and potential function of the small non-coding RNAs in the development of S. japanicum.

Results

The endogenous siRNAs were found mainly derived from transposable elements (TE) or transposons and the natural antisense transcripts (NAT). In contrast to other organisms, the TE-derived siRNAs in S. japonicum were more predominant than other sncRNAs including microRNAs (miRNAs). Further, there were distinct length and 3'end variations in the sncRNAs, which were associated with the developmental differentiation of the parasite. Among the identified miRNA transcripts, there were 38 unique to S. japonicum and 16 that belonged to 13 miRNA families are common to other metazoan lineages. These miRNAs were either ubiquitously expressed, or they exhibited specific expression patterns related to the developmental stages or sex. Genes that encoded miRNAs are mainly located in clusters within the genome of S. japonicum. However, genes within one cluster could be differentially transcribed, which suggested that individual genes might be regulated by distinct mechanisms during parasite development.

Conclusions

Many miRNA and endogenous siRNA transcripts were identified in S. japonicum and the amount of siRNA was at least 4.4 and 1.6 times more than that of miRNA in both schistosomulum and adult worm stages respectively. SiRNAs are mainly derived from transposable elements (or transposons); while natural antisense transcripts (NAT)-derived siRNAs were much less. A majority of miRNA transcripts identified in the parasite were species-specific and the expression of certain miRNAs was found developmentally regulated. Both miRNA and siRNAs are potentially important regulators in the development of schistosomal parasites.

Background

Small endogenous non-coding RNA (sncRNA) transcripts including small interfering RNA (siRNA), microRNA (miRNA) and Piwi-interacting RNA (piRNA) are critical regulators in RNA-mediated silencing in plants and metazoa[1,2]. These small RNAs, which are approximately 22 nucleotides long, guide the RNA-induced silencing complex (RISC) to their target sites and exert regulatory functions including chromatin modelling, post-transcriptional repression and mRNA destabilization, which is usually through pairing within the 3'untranslated region of target mRNAs[1].

The biogenesis and regulation of sncRNAs vary among different organisms. SiRNAs are frequently found derived from transposable elements (also called transposons), repeated sequences and antisense strands of protein-coding mRNA templates[2], while miRNAs are encoded by genes either clustered or dispersed in the genome. They can be intergenic or intragenic. MiRNA genes within introns or downstream of the rRNA genes are co-transcribed with the host genes, while the genes independently located are likely controlled by a separated mechanism. Further, the sequences of endogenous siRNA transcripts seem to be more diverse, whereas conservation patterns do exist in certain miRNAs[1,3]. For example, homologues of the temporal miRNAs, lin-4 and let-7 originally identified in Caenorhabditis elegans, have been found in a variety of eukaryotes, though they might not exert similar functions in different host cells[2].

While the knowledge regarding sncRNA biology is rapidly expanding, there is little known in schistosomes. These are parasites of at least seven developmental stages that can cause human or zoonotic schistosomiasis, which affect more than 200 million people worldwide[4]. Schistosomes are the lowest group of bilateria that diverged early from the metazoan lineage, and they are among the first animals to develop sexual dimorphism (dioecious lifestyle) and heteromorphic sex chromosomes[5]. They have seven pairs of autosomes and one pair of sex chromosomes, consisting of approximately 270 Mb of genome sequence per schistosomal parasite[6]. The unique developmental features and the availability of the genome sequences for both S. mansoni and S. japonicum has made it possible to conduct genome-wide transcriptomal and functional characterization of sncRNAs in these parasites.

Results and Discussion

Discovery of sncRNAs in S. japonicum

We identified and characterized microRNAs and endogenous siRNAs in S. japonicum both in the schistosomula (14 days post infection) and in the adult worms by direct sequencing using Solexa sequencing technology, as this technique achieves a more complete coverage of small RNA transcripts than a traditional cDNA cloning approach[7-9]. Sequencing of small RNAs yielded 6,400,876 and 5,323,610 unfiltered sequence reads from adult worms and schistosomula, respectively. After removal of reads containing ambiguous base calls, there were 5,349,115 (mixed-sex adult worms) and 4,273,194 (hepatic schistosomula) clean reads, which contained 1,193,825 and 1,135,641 unique clean reads, respectively. These unique clean reads were mapped to the genome of S. japonicum by forcing perfect alignments beginning at the first nucleotide and retaining the longest region of overlap for each read. After elimination of sequences of tRNA, rRNA and mitochondrial DNA-derived RNAs (Additional file 1), potential miRNAs and siRNAs were further analyzed[9,10].

Additional file 1. General information of the small RNA library. This file contains the reads of all small RNA transcripts identified and their relative portion in the library.

Format: PDF Size: 44KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Of the small sncRNAs transcripts identified in S. japonicum, a large portion was from potential transposable elements (TE) (Additional files 2, 3, 4, 5 and 6). Due to the fact that these sncRNAs showed similar structural signatures to siRNA transcripts reported in Drosophila[10], the small RNA transcript sequences mapped to the transposable elements (TE) and the natural antisense transcripts (NAT) were regarded as endogenous siRNAs[11,12]. Since TEs have not been previously identified in S. japonicum, we conducted TE prediction by scanning the genome sequence using REPET software. Long terminal repeats (LTR), long interspersed nucleotide elements (LINE), short interspersed elements (SINE), large terminal inverted repeats (TIR) and miniature inverted-repeat transposable elements (MITE) were the main classes of transposons identified in the genome. We found that TE-related sequences accounted for 13 percent of the S. japonicum genome (Additional file 7), which is lower than that observed in mammalian DNA that can range between 30 and 50% [13,14]. Further, SiRNAs derived from LTR, TIR and LINE accounted for the majority and transcripts from sense and antisense strands were identified (Fig. 1, Additional files 2, 3, 4, 5 and 6). Sequence scanning of the two transposons SACI-5 (LTR class retrotransposon) [15] and RTE-1 (LINE class) revealed that siRNAs were derived from more concentrated regions in the sequence (Fig. 2a, b), contrast to that found in Drosophila, in which siRNAs were generated uniformly from the whole transposon sequence[10].

thumbnailFigure 1. Percentage of siRNAs derived from different transposable elements (TE) in adult and schistosomulum stages. Majority of siRNAs were derived from LTR, LINE and TIR.

thumbnailFigure 2. The abundance of each siRNA sequence generated from the two transposons by scanning each siRNA sequence against the two transposons (A. SACI-5 and B. TRE-1). Data represents the percentage of all matched siRNAs identified in adult (red) and schistosomulum (blue) parasites.

Additional file 2. SiRNAs derived from LINE. This file contains the information of the identified transposon-LINE in the S. japonicum genome and the derived siRNAs.

Format: PDF Size: 1016KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 3. SiRNAs derived from LTR. This file contains the information of the identified transposon-LTR in the S. japonicum genome and the derived siRNAs.

Format: PDF Size: 407KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 4. SiRNAs derived from SINE. This file contains the information of the identified transposon-SINE in the S. japonicum genome and the derived siRNAs.

Format: PDF Size: 26KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 5. SiRNAs derived from TIR. This file contains the information of the identified transposon-TIR in the S. japonicum genome and the derived siRNAs.

Format: PDF Size: 34KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 6. SiRNAs derived from MITE. This file contains the information of the identified transposon-MITE in the S. japonicum genome and the derived siRNAs.

Format: PDF Size: 20KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Furthermore, siRNAs accounted for a major portion of the sncRNAs identified, while miRNAs contributed less. In adult parasites, there were 807,672 TE-derived siRNAs including 120,582 unique clean reads, and there were 495,982 miRNAs including 12,443 unique clean reads (Additional file 7). In schistosomulum, there were 1,495,593 TE-derived siRNAs including 137,439 unique clean reads, and there were 307,984 miRNAs reads including 9,231 unique clean reads (Additional file 7). Thus, more TE-derived siRNAs in schistosomulum were obtained than in the adult, while the number of miRNAs obtained from the two stages was just opposite, i.e., more miRNAs were identified in adult worms than in schistosomulum. On the other hand, miRNAs were more abundant, and siRNAs only accounted for a small portion in Drosophila somatic cells[10]. The data indicated that S. japomicum, compared to Drosophila, may possess a different regulation mechanism in sncRNA transcription.

Additional file 7. Data Summary of siRNAs and miRNAs. This file contains the summary data of siRNA and miRNA, from stage-associated differences in length distribution to the number of transcripts.

Format: PDF Size: 163KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Until recently,the endogenous siRNAs identified in Drosophila and mouse[10,16] were found to be mainly derived from transposable elements, complementary annealed transcripts and long 'fold-back' transcripts (hpRNAs). Natural antisense transcripts (NAT) generated from bidirectional gene transcription are the second source of endogenous siRNAs [17,18]. SiRNAs of approximately 21 nt derived from double-stranded RNA are associate with Ago2 in Drosophila[10]. In contrast, piRNAs of 24-32 nt, which are only functional in the germline, appear to be Dicer-independent and are associated with Piwi proteins[2]. We found Ago2 homologues in S. japonicum using a BLAST search, but we did not identify a Piwi homologue (data not shown). NAT-derived siRNAs were identified from overlapped regions of mRNAs of 153 predicted genes, and almost all of them were trans-NAT siRNAs. There were 14,427 NAT-derived siRNAs including 2,176 unique clean reads in adult worms and 5,236 NAT-derived siRNAs including 1,241 unique clean reads in schistosomulum (Additional file 7), indicating that the biogenesis of NAT-derived siRNAs in the parasite might be associated with the activation of stage-specific genes. Further analysis found that the overlapped regions of mRNA transcripts where most of the siRNAs were generated from were almost exclusively <450 bp (data not shown) and the number of siRNAs identified from sense and antisense stands were almost equal (data not shown), which indicated that they were not degraded mRNA templates. Collectively, these data revealed that the two main pathways of siRNA biogenesis in S. japonicum are processing of TE-derived transcripts and endogenous double strand RNAs (dsRNAs) and the TE-derived siRNAs accounted for the main portion of the siRNA pool, at least in the two development stages of the parasite.

Length and end variations of sncRNAs in S. japonicum

The variations in the length of sncRNAs in S. japonicum were associated with the developmental stages of the parasite. In adult worms, there was an even distribution of the amount of sncRNAs that were between 20 and 23 nt (Fig 3a,b). The TE-derived siRNAs were predominantly 20 nt in adult worms, whereas they exhibited a wider range in length in schistosomulum between 18 and 28 nt (Fig. 3a,b). The NAT-derived siRNAs were predominantly 20 nt in parasites of both schistosomulum and adult stages. Thus TE-derived siRNAs and NAT-derived siRNAs were generated under different regulatory and processing mechanisms. The size of the miRNAs also differed between the two developmental stages. In adult worms, miRNAs ranging from 19 to 23 nt were found, although the miRNAs of 22 and 23 nt accounted for the majority (Fig. 3a). However, miRNAs of 23 nt were predominant in schistosomulum (Fig. 3b) indicating that several miRNA processing pathways might be functional in different developmental stages of the parasite.

thumbnailFigure 3. Small interfering RNAs and microRNAs identified by high-throughput sequencing in S. japonicum. A and B Length and distribution of siRNAs and miRNAs in adult and schistosomulum stages. C. 3' end variations in TE and NAT-derived siRNAs and miRNAs.

In adult worm, all TE-derived siRNAs had the tendency of ending with a uracil residue, while in schistosomulum, only the 21 nt siRNAs were with this composition (Fig. 3c). Of the NAT-derived siRNAs, the 19 nt siRNAs ended with either an adenosine or a uracil in schistosomulum but not in the adult worm, while the 21 nt siRNAs from both stages had a similar sequence at the 3' ends (Fig. 3c). The sequences of both miRNA and siRNA completely match to the genome sequence suggesting that there is no sequence modification in the generation of small RNAs, which is similar to that found in Drosophila [10].

Characterization of miRNAs derived from four developmental stages of S. japonicum including cercariae, hepatic schistosomulum, adult worm and egg was performed. The miRNAs were between 22 and 23 nt in the adult and predominantly 23 nt in schistosomulum (Fig. 3a,b). Interestingly, miRNAs of all sizes predominantly ended with a uracil, and there were no observed stage-related variations (Fig. 3c3c). This was different from miRNAs derived from other metazoan species, which always begin with a uracil[10]. It remains further investigation whether this is due to the processing of species-specific pre-miRNA or due to functional restriction of the RISC.

MiRNAs are generated from hairpin structures formed by the complementary sequences within one transcript. After processing by the enzyme Dicer, two complementary sequences (strands) will be released. In this study, there was a diverse variety of miRNA isomers observed. In general, if a dominant sequence was observed, the complementary sequence was also more frequently observed but with significantly different frequency. For example, sja-mir-36a was one of the dominant miRNAs observed (Fig. 4a, Additional files 8,9), and its complementary strand sequence (miR*) was also commonly observed (highlighted in blue), but it was significantly less stable than its counterpart. Furthermore, the number of reads for multiple mature miRNA strands derived from the same precursor differed dramatically, which indicated that the hydrolysis process mediated by Dicer and its associated proteins might not be very precise and the miRNA stability might be determined by the terminal nucleotide residues.

thumbnailFigure 4. Repertoire of Sja-mir-36a isomers and a typical miRNA cluster. A. Sequences and the number of reads of the mature miRNA of Sja-mir-36a and the complementary miR* are represented in red and blue respectively. The predicted structure of the pre-miRNA is represented on the right side. B. The predicted secondary structure of a pre-miRNA containing a cluster of Sja-mir-2a, Sja-mir-2b, and Sja-novel-36, which form three hairpin stem-loop structures. The sequences of the mature miRNAs are shown in red and the miRNA* strands are shown in blue.

Additional file 8. Hairpin prediction of common miRNAs in S. japonicum. This file contains the predicted hairpin structures of common miRNAs.

Format: PDF Size: 48KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 9. Hairpin prediction of novel miRNAs in S. japonicum. This file contains the predicted hairpin structures of novel miRNAs.

Format: PDF Size: 51KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Earlier study has found that expression of miRNA genes can be both monocistronic and polycistronic[19]; however, little is known regarding the genetic regulation of miRNA expression. Certain miRNA-encoding genes are located within the introns of host genes; therefore, they might be transcriptionally regulated through their host-gene promoters[20]. In this study, we found that miRNA genes were predominantly intergenic in S. japonicum, thus most miRNA genes have their own control elements ( or promoters) in the genome [21]. Surprisingly, we found that genes within the same cluster exhibited an asymmetrical transcription pattern. For instance, Sja-mir-2a, 2b and Sja-novel-36 were situated within one cluster (Fig. 4b, Table 1 and 2); however, the number of reads of Sja-mir-2a and Sja-novel-36 could be 100 times less than that of Sja-mir-2b. Even though it cannot be ruled out that the differences in reads among the miRNAs encoded in the same cluster were due to experimental factors, it is also possible that several layers of regulation may control miRNA transcription within the same gene cluster.

Table 1. Common miRNAs identified in S. japonicum

Identification and characterization of developmental stage-associated miRNAs in S. japonicum

For further categorization of the miRNAs, alignments were performed with miRNA sequences derived from S. japonicum and miRNAs from other organisms using the ClustalW 1.8 program[22]. Similar to the study on the planarian Schmidtea mediterranea [23-25], many miRNAs discovered in the parasite were common to other metazoan lineages as well as miRNAs that were unique to S. japonicum. There were 16 miRNAs classified as common miRNAs with evolutionarily conserved characteristics related to 13 metazoan miRNA families including let-7, miR-71, miR-2, and miR-36 (Fig. 5a, Table 1and Additional file 10). Of the commen miRNAs, 13 showed high similirity to that reported in S. mediterranea [23-25]. MiRNAs unique to S. japonicum (here regarded as novel miRNAs) were also identified (Table 2and Additional files 11,12).

thumbnailFigure 5. Sequence and transcriptional analysis of miRNAs dominantly expressed in S. japonicum. A. Alignment of sja-mir-7 sequences with homologues from other organisms. The seed sequences are shadowed in a dark colour. B. Northern-blot hybridization using probes complementary to the dominant miRNAs expressed during the adult worm stage. The 70 nt pre-miRNAs were indicated by arrows. C. Quantification analysis, relative to that in the egg stage, of miRNAs dominantly expressed in four developmental stages as indicated in different colours.

Table 2. Top 20 novel abundantly-expressed miRNAs in S. japonicum

Additional file 10. Common miRNAs identified in S. japonicum. This file contains the information of common miRNAs, from contig identity to stage-associated variation in transcription.

Format: PDF Size: 29KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 11. Hairpin prediction of novel miRNAs in S. japonicum. This file contains the predicted hairpin structures of the novel miRNAs, from contig identity to stage-associated variation in transcription.

Format: PDF Size: 49KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 12. MicroRNAs in S. japonicum with high homology to that of S. mediterranea. This file contains information of the comparison between miRNAs identified in S. japonicum and that found in S. mediterranea.

Format: XLS Size: 33KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Real-time PCR was used to relatively quantify the dominantly expressed miRNAs, and the results were confirmed by Northern blotting (Fig. 5b, c). In general, more miRNAs were expressed at stages other than the egg, except Sja-mir-2b and Sja-mir-71, which were constantly expressed during all stages with less than 10-fold variations (Fig. 5c). The expression of Sja-mir-7 and 36 was predominant in the cercarial stage. MiR-7 has been reported to regulate segmentation and sensory organ development in Drosophila, and miR-36 has been associated with the transition from the embryo to L1 stage in C. elegans[24,25]. Thus it is reasonable to observe significantly increased expression of the homologous miRNAs (Sja-mir-7 and 36) in the cercariae stage of S. japonicum, since this is the stage the parasite quickly transforms into schistosomulum after invasion into a mammalian host. The miRNA let-7, which was originally identified to be dominantly expressed during the L3/L4 transition stage in C. elegans[24], played a critical role in stage-specific differentiation (timing of cell fate determination). Its homologous miRNA Sja-let-7 was also identified in S. japonicum. In contrast to its counterpart in C. elegans, Sja-let-7 was expressed in all developmental stages with a similar magnitude except during the egg stage, indicating that this miRNA has evolved to other functions than developmental timing in the schistosomal parasite (Fig. 5c). Mir-71 was previously found to be expressed in the Drosophila embryo, and it was implicated in controlling cell differentiation[26]. However, the homolog to mir-71, Sja-mir-71, was predominantly expressed in male parasites. There were also miRNAs identified that were predominantly expressed in female parasites, such as Sja-Novel-70; its expression was approximately 1000 times higher in females, as compared to cercariae, male worms and eggs. The amount of Sja-Novel-70 transcripts expressed during the female stage was one hundred thousand times higher than that in the schistosomulum stage (Fig. 5c). Even though the expression of Sja-Novel-110 was more predominant in female worms, its expression was also detected in males. The discovery of these sex-related miRNAs indicated that miRNAs might regulate parasite sexual differentiation or reproduction machinery. Since schistosomal parasites are among the first animals to develop sexual dimorphism[5], further study on this aspect may result in increased understanding of the molecular control of female and male development. Five miRNAs in S. japonicum have been previously identified via a conventional cDNA cloning and sequencing approach[27]. Compared to the previous report, apart from the reported miRNA transcripts, more miRNAs common to other organisms as well as novel miRNA transcripts were identified in this study due to the powerful coverage of the large scale sequencing capacity of Solexa technology. Further, we identified and characterized endogenous siRNAs which were dominantly derived from transposable elements in the genome of the parasite.

Conclusions

Schistosomiasis is one of the most important human helminth infections[28] that affects more than 200 million people worldwide. Currently, praziquantel is the only drug available for treating schistosomiasis. Investigations focused on parasite biology and identification of novel drug targets have become of great importance. Here we have identified a panel of common as well unique miRNAs and siRNAs in S. japonicum. The identification and characterization of siRNAs and miRNAs in the parasite and their possible biological functions have opened a new avenue towards a final dissection of parasite biology and will hopefully facilitate the discovery of potential anti-parasitic drug targets.

Methods

Parasites

Adult worms and eggs of S. japonicum were isolated from the infected rabbits at 7 weeks post-infection. Male and female adult worms were manually separated under a light microscope. Hepatic schistosomula were isolated from the infected rabbits at 2 weeks post-infection. Cercariae were harvested from the intermediate host, the snail O. hupensis hupensis. Parasites of all stages (adult worms, hepatic schistosomula, cercariae and eggs) were stored in RNAlater Solution (Invitrogen, CA, USA) before RNA preparation according the manufacture's protocol.

RNA isolation

Total RNA of S. japonicum (egg, cercariae, schistosomulum and adult worm) was extracted using Trizol reagent (Invitrogen, CA, USA). For complete precipitation of small RNA, after addition of isopropanol, the RNA extract was incubated overnight at -20°C.

Small RNA library construction

Approximately 20 μg of total RNA was size-fractionated on a 15% TBE urea polyacrylamide gel and a 15-30 base-pair fraction was excised and ligated with proprietary adapters to the 5' and 3' termini of the RNA. The RNA was converted to single-stranded cDNA using Superscript II reverse transcriptase (Invitrogen, CA, USA) and Illumina's small RNA RT-Primer following the manufacturer's instructions. The cDNA was PCR-amplified with Hotstart Phusion DNA Polymerase (New England Lab, USA) in 15 cycles using Illumina's small RNA primer set. The PCR products resolved in a 12% TBE urea polyacrylamide gel were eluted with elution buffer, quantified and subjected to Solexa's proprietary sequencing-by-synthesis method.

Bioinformatic analysis of sequencing data

Individual sequence read with the base quality scores was produced by Solexa. Adapter sequences were removed from both ends of Solexa reads. All identical sequences were counted and eliminated from the initial data set. The resulting set of the unique sequences with associated read counts was referred as sequence tags. The unique reads were mapped onto the S. japonicum genome of SGSThttp://lifecenter.sgst.cn webcite using the program SOAP[29].

To identify potential miRNA sequences, the perfectly matched reads were Blast-searched against the metazoan mature miRNA of Sanger miRBase [30] (Release 12) using the program Patscan[31]. Conserved and novel miRNAs were predicted as described[7]. Using similar approach to credibility interval approaches reported for the analysis of SAGE data [32], we employed IDEG6[33] to identify miRNAs showing statistically significant difference in relative abundance (as reflected by total count of individual sequence reads) between the two small RNA libraries from schistolomulum and adult worm stages. We used the general Chi2 method, because it resulted the most efficient tests[33]. Finally, miRNA with a P value <= 0.01 were deemed to be significantly different between the two samples.

Secondary structure predictions of the miRNAs cluster were performed using MFOLD Version 2.38 using the default folding conditions (1 M NaCl, 37°C). To generate graphical output, the mfold predictions were imported into RnaViz9 (Version 2.0). Multiple sequence alignments were produced by ClustalW 1.84 and graphical display of the output was performed using Boxshade 3.21.5.

The transposable elements (TEs) in the S. japonicum genome were predicted by using REPET http://urgi.versailles.inra.fr/development/repet/ webcite. The reads that perfectly matched S. japonicum genome were aligned to unique consensus sequences of TEs using the program SOAP[29]. The reads that perfectly matched TEs were considered as TE-derived siRNAs.

The natural antisense transcripts (NATs) were detected by aligning the predicted sequences by SGST http://lifecenter.sgst.cn webcite to each other. If a pair of overlapping sequences was matched at opposite strand and the length of overlap ≥ 100 bp and the identity of the overlaps was ≥ 95%, this pair of overlapping sequences was defined as NATs pair. The reads that perfectly matched S. japonicum genome were aligned to overlapped sequences of NAT pairs using tSOAP[29]. The reads that perfectly matched the overlap region were considered as NAT-derived siRNAs.

Northern blot analysis of miRNA transcripts

Northern blot was performed according to the method described[34]. Thirty micrograms of total RNA were resolved on a 15% denaturing polyacrylamide gel. The RNA was transferred onto a sheet of nylon membrane (Hybond-NX, Amersham/Pharmacia, Sweden) and EDC-crosslinked. Oligonucleotide probes complementary to the miRNAs were 32P-end-labeled with T4 kinase and hybridized in 2 × SSC, 1% SDS buffer containing 100 μg/ml sheared, denatured salmon sperm DNA at temperatures between 40 and 60°C. Post-hybridization washes were carried out in 0.1 × SSC, 0.2% SDS buffer at temperatures ranging from 40 and 60°C. All isotopic images were recorded using a FLA-5000 system (Fuji, Japan) with Aida Image Analyser software. RNA markers were 32P-end-labeled Decade RNA markers (Ambion, USA) prepared according to the manufacturer's instructions.

MiRNA quantification by real-time PCR

Total RNA (1 μg) was polyadenylated with Poly(A) Tailing Kit (Ambion, California, USA) at 37°C for 1 h in a 20 μl reaction mixture following the manufacturer's instruction. After phenol-chloroform extraction and ethanol precipitation, the RNAs were dissolved in diethylpyrocarbonate (DEPC)-treated water and reverse-transcribed with 200 U SuperScript™ III Reverse Transcriptase (Invitrogen) and the FirstChoice® RLM-RACE kit (Ambion, USA) according to the manufacturer's instruction. The following primers were designed as forward primers of miRNAs:

Sja-let-7 3'-GGAGGTAGTTCGTTGTGTGGT-5'; Sja-mir-2b 5'-TCACAGCCAGTATTGATGAACG-3'; Sja-mir-36 5-'CCACCGGGTAGACATTCATTCGC-3'; Sja-mir-7 5'-TGGAAGACTGGTGATATGTTGTT-3'; Sja-mir-71 5'-TGAAAGACTTGAGTAGTGAGACG-3'; Sja-novel-16 5'-TGATATGTATGGGTTACTTGGT-3'; Sja-novel-70 5'-TCAGCTGTGTTCATGTCTTCGA-3'; Sja-novel-110 5'-TGAGATCGCGATTAAAGCT-3'; Sja-novel-137 5'-AGAGGTAGTGATTCATATGACT-3' Sja-novel-148 5'-TCCCTGAGACTGATAATTGCT-3';

The sequence of 5'-GCGAGCACAGAATTAATACGAC-3' complementary to the adaptor was used as the common reverse primer.

The primers, Pf α-tubulin 5'-CATGGTAGACAACGAAGCTATTTATGA-3' and Pr α-tubulin 5'-GATTAGTGTAGGTTGGACGCTCTATG-3', were used to amplify the α-tubulin transcript as the endogenous control. PCR reactions were set up by combining 0.4 nM final concentration of each primer pair, cDNA, 12.5 μl of Power SYBR Green PCR Master Mix (ABI, USA) adjusted to a final volume of 25 μl with DEPC-treated water in triplicates. Quantification was performed using the ABI PRISM 7300 sequence detection system (Applied Biosystems, Foster City, USA). The Relative expression was analyzed using the SDS 1.4 software (Applied Biosystems, Foster City, USA). PCR efficiency for the amplicon was calculated using the method described by Ruijter et al.[35].

Authors' contributions

LH and QC conceived and designed the experiments. LH, PC, NJ performed the experiments. LH, QC analyzed the data. HW contributed reagents/materials/analysis tools. LH, PC and QC wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We are very grateful to the farmers who helped to collect the snails that infected with S. japonicum.

Funding:

This study was carried out with supports to Q. Chen from the grants (2007IPB002, 2007IPB13) to invited professor from Institute of Pathogen Biology (IPB), CAMS and Jilin University, China; the National Basic Research Program of China (973 Program, No. 2007CB513100) and Young Distinguished Scientist Grant (30625029) of the Chinese Natural Scientific Foundation.

References

  1. Post-transcriptional processing generates a diversity of 5'-modified long and short RNAs

    Nature 2009, 457(7232):1028-1032. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Ghildiyal M, Zamore PD: Small silencing RNAs: an expanding universe.

    Nature reviews 2009, 10(2):94-108. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Bartel DP: MicroRNAs: target recognition and regulatory functions.

    Cell 2009, 136(2):215-233. PubMed Abstract | Publisher Full Text OpenURL

  4. Fenwick A, Webster JP: Schistosomiasis: challenges for control, treatment and drug resistance.

    Current opinion in infectious diseases 2006, 19(6):577-582. PubMed Abstract | Publisher Full Text OpenURL

  5. Hausdorf B: Early evolution of the bilateria.

    Systematic biology 2000, 49(1):130-142. PubMed Abstract | Publisher Full Text OpenURL

  6. Hirai H, Taguchi T, Saitoh Y, Kawanaka M, Sugiyama H, Habe S, Okamoto M, Hirata M, Shimada M, Tiu WU, et al.: Chromosomal differentiation of the Schistosoma japonicum complex.

    International journal for parasitology 2000, 30(4):441-452. PubMed Abstract | Publisher Full Text OpenURL

  7. Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, et al.: Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells.

    Genome research 2008, 18(4):610-621. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Cummins JM, He Y, Leary RJ, Pagliarini R, Diaz LA Jr, Sjoblom T, Barad O, Bentwich Z, Szafranska AE, Labourier E, et al.: The colorectal microRNAome.

    Proceedings of the National Academy of Sciences of the United States of America 2006, 103(10):3687-3692. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Pfeffer S, Sewer A, Lagos-Quintana M, Sheridan R, Sander C, Grasser FA, van Dyk LF, Ho CK, Shuman S, Chien M, et al.: Identification of microRNAs of the herpesvirus family.

    Nature methods 2005, 2(4):269-276. PubMed Abstract | Publisher Full Text OpenURL

  10. Ghildiyal M, Seitz H, Horwich MD, Li C, Du T, Lee S, Xu J, Kittler EL, Zapp ML, Weng Z, et al.: Endogenous siRNAs derived from transposons and mRNAs in Drosophila somatic cells.

    Science (New York, NY) 2008, 320(5879):1077-1081. OpenURL

  11. Czech B, Malone CD, Zhou R, Stark A, Schlingeheyde C, Dus M, Perrimon N, Kellis M, Wohlschlegel JA, Sachidanandam R, et al.: An endogenous small interfering RNA pathway in Drosophila.

    Nature 2008, 453(7196):798-802. PubMed Abstract | Publisher Full Text OpenURL

  12. Kawamura Y, Saito K, Kin T, Ono Y, Asai K, Sunohara T, Okada TN, Siomi MC, Siomi H: Drosophila endogenous small RNAs bind to Argonaute 2 in somatic cells.

    Nature 2008, 453(7196):793-797. PubMed Abstract | Publisher Full Text OpenURL

  13. Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al.: Initial sequencing and analysis of the human genome.

    Nature 2001, 409(6822):860-921. PubMed Abstract | Publisher Full Text OpenURL

  14. Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, et al.: Initial sequencing and comparative analysis of the mouse genome.

    Nature 2002, 420(6915):520-562. PubMed Abstract | Publisher Full Text OpenURL

  15. DeMarco R, Machado AA, Bisson-Filho AW, Verjovski-Almeida S: Identification of 18 new transcribed retrotransposons in Schistosoma mansoni.

    Biochemical and biophysical research communications 2005, 333(1):230-240. PubMed Abstract | Publisher Full Text OpenURL

  16. Watanabe T, Totoki Y, Toyoda A, Kaneda M, Kuramochi-Miyagawa S, Obata Y, Chiba H, Kohara Y, Kono T, Nakano T, et al.: Endogenous siRNAs from naturally formed dsRNAs regulate transcripts in mouse oocytes.

    Nature 2008, 453(7194):539-543. PubMed Abstract | Publisher Full Text OpenURL

  17. Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu JK: Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis.

    Cell 2005, 123(7):1279-1291. PubMed Abstract | Publisher Full Text OpenURL

  18. Okamura K, Balla S, Martin R, Liu N, Lai EC: Two distinct mechanisms generate endogenous siRNAs from bidirectional transcription in Drosophila melanogaster.

    Nature structural & molecular biology 2008, 15(6):581-590. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Zeng Y: Principles of micro-RNA production and maturation.

    Oncogene 2006, 25(46):6156-6162. PubMed Abstract | Publisher Full Text OpenURL

  20. Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel genes coding for small expressed RNAs.

    Science (New York, NY) 2001, 294(5543):853-858. OpenURL

  21. He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation.

    Nature reviews 2004, 5(7):522-531. PubMed Abstract | Publisher Full Text OpenURL

  22. Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs.

    Nucleic acids research 2003, 31(13):3497-3500. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Palakodeti D, Smielewska M, Graveley BR: MicroRNAs from the Planarian Schmidtea mediterranea: a model system for stem cell biology.

    RNA (New York, NY) 2006, 12(9):1640-1649. OpenURL

  24. Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G: The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans.

    Nature 2000, 403(6772):901-906. PubMed Abstract | Publisher Full Text OpenURL

  25. Lau NC, Lim LP, Weinstein EG, Bartel DP: An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans.

    Science (New York, NY) 2001, 294(5543):858-862. OpenURL

  26. Aboobaker AA, Tomancak P, Patel N, Rubin GM, Lai EC: Drosophila microRNAs exhibit diverse spatial expression patterns during embryonic development.

    Proceedings of the National Academy of Sciences of the United States of America 2005, 102(50):18017-18022. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Xue X, Sun J, Zhang Q, Wang Z, Huang Y, Pan W: Identification and characterization of novel microRNAs from Schistosoma japonicum.

    PLoS ONE 2008, 3(12):e4034. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. King CH, Dickman K, Tisch DJ: Reassessment of the cost of chronic helmintic infection: a meta-analysis of disability-related outcomes in endemic schistosomiasis.

    Lancet 2005, 365(9470):1561-1569. PubMed Abstract | Publisher Full Text OpenURL

  29. Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program.

    Bioinformatics (Oxford, England) 2008, 24(5):713-714. PubMed Abstract | Publisher Full Text OpenURL

  30. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics.

    Nucleic acids research 2008, (36 Database):D154-158. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Dsouza M, Larsen N, Overbeek R: Searching for patterns in genomic data.

    Trends Genet 1997, 13(12):497-498. PubMed Abstract | Publisher Full Text OpenURL

  32. Vencio RZ, Brentani H, Pereira CA: Using credibility intervals instead of hypothesis tests in SAGE analysis.

    Bioinformatics (Oxford, England) 2003, 19(18):2461-2464. PubMed Abstract | Publisher Full Text OpenURL

  33. Romualdi C, Bortoluzzi S, D'Alessi F, Danieli GA: IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments.

    Physiological genomics 2003, 12(2):159-162. PubMed Abstract OpenURL

  34. Pall GS, Codony-Servat C, Byrne J, Ritchie L, Hamilton A: Carbodiimide-mediated cross-linking of RNA to nylon membranes improves the detection of siRNA, miRNA and piRNA by northern blot.

    Nucleic acids research 2007, 35(8):e60. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Ruijter JM, Ramakers C, Hoogaars WM, Karlen Y, Bakker O, van den Hoff MJ, Moorman AF: Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data.

    Nucleic acids research 2009, 37(6):e45. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL