BACKGROUND
human schistosomiasis, which is second only to malaria in terms of morbidity and mortality, is a chronic debilitating disease caused by infections of schistosoma species that vary depending on the endemic region of the parasites
<cit> .
three principal schistosoma species can infect humans and cause severe diseases: schistosoma japonicum, schistosoma mansoni, and schistosoma haematobium.
s. japonicum is the causative agent of zoonotic schistosomiasis, affecting millions of people in several east and southeast asian countries.
despite the availability of a highly effective chemotherapeutic drug , the high re-infection rates in humans and animals plus the requirement of frequent administration of the agent still limits the overall success of chemotherapy and disease control efforts.
novel targets for drug and vaccine development remain to be defined for optimal treatment and disease prevention; however, the lack of knowledge about this parasite’s biology remains a hurdle.
schistosoma parasites can persist in a mammalian host for decades in the presence of the host immune system, and current knowledge about the mechanism of parasitization is still fragmented.
what is known is that the successful host-evasion mechanisms of the parasite involve the inert tegument that covers the surface in most developmental stages, the recruitment of host components to the surface, and the expression of various antigens and immune-regulating factors
.
schistosoma parasites have a complicated developmental and biological cycle.
they are among the few platyhelminth parasites to adopt a dioecious lifestyle and possess heteromorphic sex chromosomes.
the genome of s. japonicum contains eight pairs of chromosomes comprising seven pairs of autosomes and one pair of sexual chromosomes, with an estimated 397 mb containing primarily  <dig>  protein-coding sequences
<cit>  that account for 4% of the genome.
the decoding and availability of the genome sequences of the three most pathogenic parasites, s. mansoni, s. japonicum, and s. haematobium, has proved pivotal for the systematic dissection of the parasite biology
.
deep transcriptome sequencing  with next-generation sequencing technologies has provided unprecedented opportunities to investigate the genome-wide transcriptional property of many species
.
this technique allows for the survey of the entire transcriptome in a very high-throughput and quantitative manner, making it possible to identify exons and introns, map their boundaries and the 5′ and 3′ ends of genes, and understand the complexity of genome organization and activity comprehensively.
a majority of eukaryotic protein-coding genes contain intron sequences that must be removed by splicing after transcription from the dna templates.
however, some pre-mrnas can be processed alternatively by the splicing out or retention of the transcript regions of exons or introns.
this alternative splicing allows individual genes to produce two or more variant mrna templates, which in many cases encode functionally distinct proteins.
alternative splicing is an integrated process in regulation of gene transcription and expression and results in structural and functional diversity of molecules
<cit> .
because of the powerful readout of rna-seq, which can generate many sequence reads that span exon–exon junctions, rna transcripts generated from different splicing events can be identified
<cit> .
so far, five basic modes of alternative splicing are generally recognized: exon skipping , mutually exclusive exons, alternative donor sites , alternative acceptor sites , and intron retention
<cit> .
es, also called exon cassette, indicates that an exon is spliced out from the primary transcript and occurs most commonly in mammalian cells
<cit> .
in the event of mutually exclusive exons, only one of the two exons is retained in mrnas after splicing.
an ads results when an alternative 5′ splice junction  is used, leading to a change in the 3′ boundary of the upstream exon.
an aas arises when an alternative 3′ splice junction  is used, leading to a change in the 5′ boundary of the downstream exon .
ir occurs when a sequence is spliced out as an intron or simply retained and is distinguished from es because introns do not flank the retained sequence.
the retained transcript of the intron region in most cases encodes amino acids in-frame with the neighboring exons
<cit> .
recent results have suggested that schistosomes create multiple protein variants by splicing micro-exon gene transcripts, which might be involved in immune evasion mechanisms
<cit> ; however, the general feature of alternative splicing in the parasites remains understudied.
here, we investigated the alternative splicing of transcripts in both male and female s. japonicum after deep rna sequencing.
we found that the gene transcripts were diversely processed and that four types of rna splicing were identifiable after transcription of the genome.figure  <dig>
schematic illustration of alternative splicing.
a) exon skipping.
gene a forms two different transcripts; the first transcript has a new exon compared to the second transcript, the new exon is an inclusive exon, and the other two exons are constitutive.
b) intron retention.
gene b forms two different transcripts; the second transcript is a new exon formed from retained intron and exons on both sides.
c) alternative donor site.
gene c forms two different transcripts; the difference is one exon of an alternative 5′ splice site of the second transcript extended.
d) alternative acceptor site.
gene d forms two different transcripts; the difference is one exon of the alternative 3′ splice site of the second transcript extended.
e) mutually exclusive exon.
gene e forms two different transcripts; the different exon is an inclusive exon, the same exon is a constitutive exon, and two transcripts have different inclusive exons.
methods
parasites and rna purification
schistosoma japonicum–infected oncomelania hupensis were purchased from jiangxi institute of parasitic disease, nanchang, china.
cercariae were freshly shed from the infected snails.
one new zealand white female rabbit was percutaneously infected with ~ <dig>  s. japonicum cercariae, as described previously
<cit> .
mature adult parasites were isolated at 6 weeks post-infection from the rabbit by flushing the blood vessels with phosphate-buffered saline, as described previously
.
male and female parasites were manually separated with the aid of a light microscope.
total rna from the parasites was purified with trizol reagent , and contaminating genomic dna was removed using the rnase-free dnase set .
rna quantification and quality were examined with a nanodrop nd- <dig> spectrophotometer  and standard agarose gel electrophoresis.
all rna samples were stored at -80°c until use.
library preparation and sequencing
polyadenylated rna samples from adult male and female s. japonicum parasites were isolated from total rna using oligo- conjugated magnetic beads .
the mrna was interrupted into short fragments by adding the fragmentation buffer provided by the manufacturer .
with these short fragments as templates, random hexamer primers were used to synthesize the first-strand cdna.
the second-strand cdna was synthesized using buffer, dntps, rnase h, and dna polymerase i, respectively.
short fragments were purified following instructions accompanying the kit , and double-stranded cdnas were end-repaired according to manufacturer-recommended protocols, followed by connection with illumina adapters .
the fragments were first amplified by pcr.
purified cdna fragments were pooled and indexed and loaded onto one lane of an illumina ga iix flow cell.
a total of  <dig> pair-end sequencing cycles were carried out.
cluster formation, primer hybridization, and pair-end sequencing were performed according to the provided protocols
<cit> .
sequence analysis
low-quality reads , reads in which unknown bases represented more than 10%, and adapter sequences were removed from the reads, and the clean reads were mapped onto the s. japonicum genome of sgst,  by tophat
<cit> , then assembled with cufflinks
<cit>  to construct unique transcript sequences using the parameter: -g –b –u –o .
the cufflinks assembler is freely available at
http://cufflinks.cbcb.umd.edu/.
cuffcompare
<cit>  was used to compare the assembled transcripts of each library to the referenced annotated genes and build a non-redundant transcript dataset among the libraries.
then, cuffdiff was used to find significant changes in gene expression level
<cit> .
we used fdr to correct p values and obtained q values; for q value ≤5%, we considered the genes to be differentially expressed .
several perl scripts were written to summarize the splicing forms in each library.
the following algorithms were used to detect alternative splicing events.
first, junction sites, which give information about boundaries and combinations of different exons in a transcript, were detected by tophat .
then, all junction sites of the same gene were used to distinguish the type of alternative splicing event
<cit>  .
functional annotation and classification
transcripts were first compared using the kyoto encyclopedia of genes and genomes database
<cit>  with blastx
<cit>  at e values ≤ 1e- <dig>  a perl script was used to retrieve ko  information from the blast result and establish pathway associations between unigene and the database.
interpro
<cit>  domains were annotated by interproscan
<cit>  , and functional assignments were mapped onto gene ontology
<cit> .
wego
<cit>  was employed to do go classification and draw the go tree.
the significance analysis of functional pathways was performed using ideg6
<cit> .
to identify pseudogenes in the s. japonicum genome, we used pseudopipe
<cit> .
the assembled transcripts that fell into or included the predicted position of pseudogenes were designated as pseudogenes.
wego was used for the go classification.
non-coding rna annotation
rfam
<cit>   databases were used to annotate the non-coding transcripts.
the assembled novel transcripts were compared to rfam by blast at e values ≤ 1e- <dig>
verification of alternative splicing transcripts by rt-pcr and sequencing
genomic dna of s. japonicum  was purified with the dneasy blood & tissue kit  according to the manufacturer’s instructions.
total rna was prepared using trizol reagent , as previously described
<cit> , and contaminating genomic dna was removed with the rnase-free dnase set .
pcr was conducted in triplicate, and each reaction involved  <dig> amplification cycles on an applied biosystems  <dig> pcr system .
the 20 μl reaction system contained 50 ng of total rna  or 80 ng dna,  <dig>  μm of each primer, and 10 μl of premix ex taq .
the reaction conditions were as follows: 94°c for 3 min; 35 cycles of 94°c, 30 s; 55°c, 30 s; and 72°c, 90 s; and then 10 min at 72°c.
an 8 μl aliquot of each pcr sample was then subjected to electrophoresis in a  <dig> % agarose gel.
the rt-pcr primer sequences are listed in additional file
3: table s <dig>
RESULTS
identification of a large number of novel transcripts from un-annotated genome loci following deep sequencing of the s. japonicumtranscriptome
in this study, we determined the transcriptomes of male and female adult worms of s japonicum by high-throughput rna-seq with poly-a–purified rna samples.
a total of  <dig> , <dig> and  <dig> , <dig> paired reads were obtained, with a total length of  <dig> , <dig>  and  <dig> , <dig>  bp from female and male worms, respectively .
the number of predicted genes of female and male worms of s. japonicum was  <dig>  and  <dig> , respectively, which was more than that predicted based on the genome sequence
<cit> .
of the  <dig>  genes predicted in the female parasite, a total of  <dig>  were known and  <dig>  were novel, while of the  <dig>  predicted genes in the male parasite, a total of  <dig>  were known and  <dig>  were novel.
the number of predicted transcripts in the two libraries of female and male parasites was  <dig>  and  <dig> , respectively, with  <dig>  known and  <dig>  novel transcripts in females and  <dig>  known and  <dig>  novel transcripts in males.
the finding of so many novel transcripts should assist with the upgrade or reassembling of the genome sequence of s. japonicum
<cit> .
however, the novel transcripts may also be generated by alternative post-transcriptional rna processing or alternative splicing
<cit> .
indeed, we found  <dig>  multi-transcription loci in female and  <dig>  in male parasites, with about  <dig>  transcripts per locus in both sexes .
thus, the general transcription activity of both sexes was diverse but in a similar manner, confirming earlier studies
<cit> .
however, the sequence data generated from this study was much more than that of earlier studies
<cit> .
this was likely due to the approach applied in this study which is technically advanced than the digital gene expression technique.
all sequence data have been deposited in the database  with an accession number of gse <dig>
the sequence reads can be classified into four types: exons, introns, intergenic, and spliced.
the proportions of transcripts in female parasites from exons, introns, and intergenic loci were 56%, 7%, and 24%, respectively, and 13% of the transcripts were generated by alternative splicing .
similarly, in male parasites, the percentages of transcripts from exons, introns, and intergenic loci were 57%, 7%, and 22% respectively; 14% of the transcripts were generated by alternative splicing .table  <dig>
summary data of the transcriptome analysis
summary statistics of the reads and the reads’ total length from the two libraries female  or male  parasite were listed.
proportions of sequence reads  generated from different genetic regions in the genomes of female and male
s. japonicum
.
more than 50% of the transcripts were generated from exons while transcripts from intron, intergenic regions, and splicing events were around 7%, 22%, and 13%, respectively.
alternative splicing in s. japonicum
four types of alternative splicing in both female and male worms of s. japonicum  were identified, including es, ir, ads, and aas.
of the alternative splicing events, aas and ads were more common than the other two types, suggesting that the gene regulation mechanism of the schistosoma parasite is diverged from that of the mammalian taxa, in which es has been more commonly observed
<cit> .
in female s. japonicum, a total of  <dig>  alternative splicing events were bioinformatically predicted while in male worms, a total of  <dig>  were predicted .
the percentage of different alternative splicing events was similar between the two sexes ; however, the genes undergoing alternative splicing were not necessarily the same between them .table  <dig>
statistics for alternative splicing events
sj-f and sj-m are the detected alternative splicing events of the female  or male  parasite.
reference* is the reference annotated genes of the two libraries for the female  or male  parasite.
to confirm experimentally the prediction of the alternative splicing events in the bioinformatic analysis, eight transcripts in which alternatively spliced fragments were more than 100 bp were chosen randomly.
transcripts generated by es skipping of five genes and transcripts generated by ir of three genes were validated by pcr and rt-pcr.
the five genes with es activity included one that encodes the protein c14orf <dig> homolog; a novel gene; s. japonicum zinc finger ccch domain-containing protein 5; s. japonicum cell division cycle and apoptosis regulator protein 1; and s. japonicum protein phosphatase  <dig> regulatory subunit sds <dig>  the three genes with ir were respectively beta-amyloid binding protein , s. japonicum ipr <dig> pdz/dhr/glgf domain-containing protein, and deoxyribodipyrimidine photolyase.
the amplicons of all rt-pcr reactions were cloned and sequenced and were correlated with the predicated alternative splicing events , suggesting that the bioinformatic prediction based on the primary sequencing data was reliable.figure  <dig>
sequence mapping and verification of  <dig> genes with exon skipping events detected by rna-seq by pcr and rt-pcr.
the expression profiles of the same gene in female  and male  parasites were placed under the line representing the chromosome position.
the black lines represent original annotated gene structures , while the active transcripts in red and blue identified from the same genes in female and male parasites are underneath.
the five genes and transcripts  with exon skipping events were confirmed by pcr and rt-pcr .
ge indicates pcr products amplified from genomic dna, and ce indicates pcr products amplified from cdna.
red arrows indicate transcripts generated by exon skipping, and green arrows indicate primer locations.
sequence mapping and verification of three genes with intron retention events detected by rna-seq by pcr and rt-pcr.
the three genes and transcripts  with intron retention events were confirmed by pcr and rt-pcr .
gi indicates pcr products amplified from genomic dna, and ci indicates pcr products amplified from cdna.
red arrows indicate transcripts generated by intron retention, and green arrows indicate primer locations.
functional category of alternatively spliced genes in s. japonicum
after mapping of the rna-seq reads to the s. japonicum reference genome, transcripts were assembled and their relative abundances were calculated.
cuffdiff was used to find significant changes in gene level expression in the two libraries
<cit> .
genes subject to alternative splicing and showing significant differences in expression in the two libraries are listed in table
<dig> and additional file
6: figure s <dig>  genes related to the function of genetic information processing were found to be more biased to the female parasite while genes related to the environmental information processing were more active in the male parasites .
these data reflect the biology of the two sexes of the parasite.
the female parasites, which are kept in the cavity of the male parasites, are more active in the reproduction process while the male parasites are principally responsible in the host–parasite interaction.table  <dig>
go classification statistics of alternatively processed genes that were differentially expressed in female and male
s. japonicum
as alternative splicing.
diff differentially expressed.
total as & diff, genes with both alternative splicing and differential expression in female and male s. japonicum.
identification of novel transcripts from intergenic regions and previously determined pseudogenes
one of the advantages of transcriptomic analysis is that it allows identification of novel transcripts that may not be predicted based on genomic sequences.
the novel sequences can thus provide a powerful tool for re-annotation of the genome of s. japonicum, which has been poorly assembled .
of the novel polyadenylated sequences, two classes of transcripts have been identified: one that does not map either to regions of the genome corresponding to annotated genes or to the untranslated regions, and another that maps to previously annotated pseudogenes
<cit> .
we identified  <dig>  novel transcripts that completely matched the previously annotated intergenic regions of the genome.
the length of these transcripts was from  <dig> to  <dig>  bp, with an average length of  <dig>  bp .
it has been reported that the s. mansoni genome contains many small open reading frames .
our results indicated that the small transcripts derived from both intergenic and "pseudogenes" in s. japonicum may encode important functions, as reported for the human genome
<cit> .
further, 31%  of these transcripts had at least one complete open reading frame that could be translated into proteins; the other 69% were not annotated .
thus, of the  <dig>  novel transcripts, at least  <dig>  genomic sequences corresponding to the transcripts can be re-annotated as protein-coding genes.
a total of  <dig> transcripts were mapped to the non-coding rna database of the rfam
, which is also frequently used as a source of high-quality alignments for training and benchmarking rna sequence analysis software tools.
these transcripts were found as either micrornas or ribosomal or other non-coding rnas .
they were likely the contaminated sequences which were not completely depleted during mrna purification process.figure  <dig>
numbers of non-coding rna transcripts identified in the transcriptomes of female and male
s. japonicum
.
furthermore, among the  <dig>  novel transcripts, we detected  <dig>  that were derived from pseudogenes; of these,  <dig> were derived from annotated pseudogenes, and the rest were from unannotated pseudogenes .
pseudogenes can be transcribed from either direction, which contributed to the templates for generation of small endogenous interfering rnas in s. japonicum, which is more common in transposable elements
.
the identification of these pseudogene-derived transcripts suggested that all sequences are polyadenylated and that there is no discrimination between coding and non-coding transcripts in the rna polyadenylation process in schistosoma parasites.
on the other hand, complete reading frames were indeed identified in a number of pseudogene-derived transcripts that encoded proteins with known functions .
thus, these "pseudogene" genes can be re-annotated as protein-coding genes.figure  <dig>
go categories of transcripts with a complete open reading frame and derived from pseudogenes.
the percentages of the transcripts encoding the same proteins with similar function are indicated on the left while the numbers of the transcripts identified are indicated on the right.
CONCLUSIONS
in summary, by using rna-seq technology, we obtained the global transcriptomes of male and female s. japonicum.
approximately 80% of the total reference genes  were expressed in the adult stage of the parasite, representing the majority of the transcriptomes.
these results further provide a comprehensive view of the global transcriptome of s. japonicum.
the findings of a substantial level of alternative splicing events dynamically occurring in the parasitization in the mammalian hosts of the s. japonicum suggest complicated transcriptional and post-transcriptional regulatory mechanisms employed by the parasite.
the data should not only significantly improve the re-annotation of the genome sequences but also should provide new information about the biology of the parasite.
electronic supplementary material
additional file 1:
differentially expressed genes between females and males.
additional file 2: figure s1: junction sites.
additional file 3: table s1: primers and sequences for verification of the alternative splicing events.
additional file 4:
total transcripts identified in female and male parsites.
additional file 5:
genes with alternative splicing in female and male parsites.
additional file 6: figure s2: go category of the genes that were alternatively spliced and also differentially transcribed.
additional file 7:
genes related to the function of genetic information processing and environmental information.
additional file 8:
novel transcripts identified.
additional file 9:
noncoding transcripts identified in the female
s. japonicum
.
additional file 10:
noncoding transcripts identified in the male
s. japonicum.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
xp, nh and qc conceived and designed the experiments.
xp, nh, pc, sl, and cw performed the experiments.
xp, nh and qc analyzed the data.
xp and qc wrote the manuscript.
all authors read and approved the final manuscript.
we appreciate very much the bioinformatic support of dr. haibo sun at miningene biotechnology and the efforts of technicians at shenzhen bgi for solexa sequencing.
we also thank the schistosoma japonicum genome sequencing and functional analysis consortium for making the s. japonicum genome available on line publicly.
funding
this study was supported by the national natural science foundation of china , the intramural grant from institute of pathogen biology, cams , the national s & t major program  and the program for changjiang scholars and innovative research team in university.

 human schistosomiasis, which is second only to malaria in terms of morbidity and mortality, is a chronic debilitating disease caused by infections of schistosoma species that vary depending on the endemic region of the parasites. three principal schistosoma species can infect humans and cause severe diseases: schistosoma japonicum, schistosoma mansoni, and schistosoma haematobium. s. japonicum is the causative agent of zoonotic schistosomiasis, affecting millions of people in several east and southeast asian countries. despite the availability of a highly effective chemotherapeutic drug, the high re-infection rates in humans and animals plus the requirement of frequent administration of the agent still limits the overall success of chemotherapy and disease control efforts. novel targets for drug and vaccine development remain to be defined for optimal treatment and disease prevention; however, the lack of knowledge about this parasite’s biology remains a hurdle. schistosoma parasites can persist in a mammalian host for decades in the presence of the host immune system, and current knowledge about the mechanism of parasitization is still fragmented. what is known is that the successful host-evasion mechanisms of the parasite involve the inert tegument that covers the surface in most developmental stages, the recruitment of host components to the surface, and the expression of various antigens and immune-regulating factors. schistosoma parasites have a complicated developmental and biological cycle. they are among the few platyhelminth parasites to adopt a dioecious lifestyle and possess heteromorphic sex chromosomes. the genome of s. japonicum contains eight pairs of chromosomes comprising seven pairs of autosomes and one pair of sexual chromosomes, with an estimated 397 mb containing primarily protein-coding sequences that account for 4% of the genome. the decoding and availability of the genome sequences of the three most pathogenic parasites, s. mansoni, s. japonicum, and s. haematobium, has proved pivotal for the systematic dissection of the parasite biology. deep transcriptome sequencing with next-generation sequencing technologies has provided unprecedented opportunities to investigate the genome-wide transcriptional property of many species. this technique allows for the survey of the entire transcriptome in a very high-throughput and quantitative manner, making it possible to identify exons and introns, map their boundaries and the 5′ and 3′ ends of genes, and understand the complexity of genome organization and activity comprehensively. a majority of eukaryotic protein-coding genes contain intron sequences that must be removed by splicing after transcription from the dna templates. however, some pre-mrnas can be processed alternatively by the splicing out or retention of the transcript regions of exons or introns. this alternative splicing allows individual genes to produce two or more variant mrna templates, which in many cases encode functionally distinct proteins. alternative splicing is an integrated process in regulation of gene transcription and expression and results in structural and functional diversity of molecules. because of the powerful readout of rna-seq, which can generate many sequence reads that span exon–exon junctions, rna transcripts generated from different splicing events can be identified. so far, five basic modes of alternative splicing are generally recognized: exon skipping, mutually exclusive exons, alternative donor sites, alternative acceptor sites, and intron retention. es, also called exon cassette, indicates that an exon is spliced out from the primary transcript and occurs most commonly in mammalian cells. in the event of mutually exclusive exons, only one of the two exons is retained in mrnas after splicing. an ads results when an alternative 5′ splice junction is used, leading to a change in the 3′ boundary of the upstream exon. an aas arises when an alternative 3′ splice junction is used, leading to a change in the 5′ boundary of the downstream exon. ir occurs when a sequence is spliced out as an intron or simply retained and is distinguished from es because introns do not flank the retained sequence. the retained transcript of the intron region in most cases encodes amino acids in-frame with the neighboring exons. recent results have suggested that schistosomes create multiple protein variants by splicing micro-exon gene transcripts, which might be involved in immune evasion mechanisms; however, the general feature of alternative splicing in the parasites remains understudied. here, we investigated the alternative splicing of transcripts in both male and female s. japonicum after deep rna sequencing. we found that the gene transcripts were diversely processed and that four types of rna splicing were identifiable after transcription of the genome.figure schematic illustration of alternative splicing. a) exon skipping. gene a forms two different transcripts; the first transcript has a new exon compared to the second transcript, the new exon is an inclusive exon, and the other two exons are constitutive. b) intron retention. gene b forms two different transcripts; the second transcript is a new exon formed from retained intron and exons on both sides. c) alternative donor site. gene c forms two different transcripts; the difference is one exon of an alternative 5′ splice site of the second transcript extended. d) alternative acceptor site. gene d forms two different transcripts; the difference is one exon of the alternative 3′ splice site of the second transcript extended. e) mutually exclusive exon. gene e forms two different transcripts; the different exon is an inclusive exon, the same exon is a constitutive exon, and two transcripts have different inclusive exons. methods parasites and rna purification schistosoma japonicum–infected oncomelania hupensis were purchased from jiangxi institute of parasitic disease, nanchang, china. cercariae were freshly shed from the infected snails. one new zealand white female rabbit was percutaneously infected with ~  s. japonicum cercariae, as described previously. mature adult parasites were isolated at 6 weeks post-infection from the rabbit by flushing the blood vessels with phosphate-buffered saline, as described previously. male and female parasites were manually separated with the aid of a light microscope. total rna from the parasites was purified with trizol reagent, and contaminating genomic dna was removed using the rnase-free dnase set. rna quantification and quality were examined with a nanodrop nd- spectrophotometer and standard agarose gel electrophoresis. all rna samples were stored at -80°c until use. library preparation and sequencing polyadenylated rna samples from adult male and female s. japonicum parasites were isolated from total rna using oligo- conjugated magnetic beads. the mrna was interrupted into short fragments by adding the fragmentation buffer provided by the manufacturer. with these short fragments as templates, random hexamer primers were used to synthesize the first-strand cdna. the second-strand cdna was synthesized using buffer, dntps, rnase h, and dna polymerase i, respectively. short fragments were purified following instructions accompanying the kit, and double-stranded cdnas were end-repaired according to manufacturer-recommended protocols, followed by connection with illumina adapters. the fragments were first amplified by pcr. purified cdna fragments were pooled and indexed and loaded onto one lane of an illumina ga iix flow cell. a total of pair-end sequencing cycles were carried out. cluster formation, primer hybridization, and pair-end sequencing were performed according to the provided protocols. sequence analysis low-quality reads, reads in which unknown bases represented more than 10%, and adapter sequences were removed from the reads, and the clean reads were mapped onto the s. japonicum genome of sgst, by tophat, then assembled with cufflinks to construct unique transcript sequences using the parameter: -g –b –u –o. the cufflinks assembler is freely available at http://cufflinks.cbcb.umd.edu/. cuffcompare was used to compare the assembled transcripts of each library to the referenced annotated genes and build a non-redundant transcript dataset among the libraries. then, cuffdiff was used to find significant changes in gene expression level. we used fdr to correct p values and obtained q values; for q value ≤5%, we considered the genes to be differentially expressed. several perl scripts were written to summarize the splicing forms in each library. the following algorithms were used to detect alternative splicing events. first, junction sites, which give information about boundaries and combinations of different exons in a transcript, were detected by tophat. then, all junction sites of the same gene were used to distinguish the type of alternative splicing event. functional annotation and classification transcripts were first compared using the kyoto encyclopedia of genes and genomes database with blastx at e values ≤ 1e- a perl script was used to retrieve ko information from the blast result and establish pathway associations between unigene and the database. interpro domains were annotated by interproscan, and functional assignments were mapped onto gene ontology. wego was employed to do go classification and draw the go tree. the significance analysis of functional pathways was performed using ideg6. to identify pseudogenes in the s. japonicum genome, we used pseudopipe. the assembled transcripts that fell into or included the predicted position of pseudogenes were designated as pseudogenes. wego was used for the go classification. non-coding rna annotation rfam databases were used to annotate the non-coding transcripts. the assembled novel transcripts were compared to rfam by blast at e values ≤ 1e- verification of alternative splicing transcripts by rt-pcr and sequencing genomic dna of s. japonicum was purified with the dneasy blood & tissue kit according to the manufacturer’s instructions. total rna was prepared using trizol reagent, as previously described, and contaminating genomic dna was removed with the rnase-free dnase set. pcr was conducted in triplicate, and each reaction involved amplification cycles on an applied biosystems pcr system. the 20 μl reaction system contained 50 ng of total rna or 80 ng dna,  μm of each primer, and 10 μl of premix ex taq. the reaction conditions were as follows: 94°c for 3 min; 35 cycles of 94°c, 30 s; 55°c, 30 s; and 72°c, 90 s; and then 10 min at 72°c. an 8 μl aliquot of each pcr sample was then subjected to electrophoresis in a % agarose gel. the rt-pcr primer sequences are listed in additional file 3: table s  identification of a large number of novel transcripts from un-annotated genome loci following deep sequencing of the s. japonicumtranscriptome in this study, we determined the transcriptomes of male and female adult worms of s japonicum by high-throughput rna-seq with poly-a–purified rna samples. a total of, and, paired reads were obtained, with a total length of, and,  bp from female and male worms, respectively. the number of predicted genes of female and male worms of s. japonicum was and, respectively, which was more than that predicted based on the genome sequence. of the genes predicted in the female parasite, a total of were known and were novel, while of the predicted genes in the male parasite, a total of were known and were novel. the number of predicted transcripts in the two libraries of female and male parasites was and, respectively, with known and novel transcripts in females and known and novel transcripts in males. the finding of so many novel transcripts should assist with the upgrade or reassembling of the genome sequence of s. japonicum. however, the novel transcripts may also be generated by alternative post-transcriptional rna processing or alternative splicing. indeed, we found multi-transcription loci in female and in male parasites, with about transcripts per locus in both sexes. thus, the general transcription activity of both sexes was diverse but in a similar manner, confirming earlier studies. however, the sequence data generated from this study was much more than that of earlier studies. this was likely due to the approach applied in this study which is technically advanced than the digital gene expression technique. all sequence data have been deposited in the database with an accession number of gse the sequence reads can be classified into four types: exons, introns, intergenic, and spliced. the proportions of transcripts in female parasites from exons, introns, and intergenic loci were 56%, 7%, and 24%, respectively, and 13% of the transcripts were generated by alternative splicing. similarly, in male parasites, the percentages of transcripts from exons, introns, and intergenic loci were 57%, 7%, and 22% respectively; 14% of the transcripts were generated by alternative splicing.table summary data of the transcriptome analysis summary statistics of the reads and the reads’ total length from the two libraries female or male parasite were listed. proportions of sequence reads generated from different genetic regions in the genomes of female and male s. japonicum. more than 50% of the transcripts were generated from exons while transcripts from intron, intergenic regions, and splicing events were around 7%, 22%, and 13%, respectively. alternative splicing in s. japonicum four types of alternative splicing in both female and male worms of s. japonicum were identified, including es, ir, ads, and aas. of the alternative splicing events, aas and ads were more common than the other two types, suggesting that the gene regulation mechanism of the schistosoma parasite is diverged from that of the mammalian taxa, in which es has been more commonly observed. in female s. japonicum, a total of alternative splicing events were bioinformatically predicted while in male worms, a total of were predicted. the percentage of different alternative splicing events was similar between the two sexes; however, the genes undergoing alternative splicing were not necessarily the same between them.table statistics for alternative splicing events sj-f and sj-m are the detected alternative splicing events of the female or male parasite. reference* is the reference annotated genes of the two libraries for the female or male parasite. to confirm experimentally the prediction of the alternative splicing events in the bioinformatic analysis, eight transcripts in which alternatively spliced fragments were more than 100 bp were chosen randomly. transcripts generated by es skipping of five genes and transcripts generated by ir of three genes were validated by pcr and rt-pcr. the five genes with es activity included one that encodes the protein c14orf homolog; a novel gene; s. japonicum zinc finger ccch domain-containing protein 5; s. japonicum cell division cycle and apoptosis regulator protein 1; and s. japonicum protein phosphatase regulatory subunit sds the three genes with ir were respectively beta-amyloid binding protein, s. japonicum ipr pdz/dhr/glgf domain-containing protein, and deoxyribodipyrimidine photolyase. the amplicons of all rt-pcr reactions were cloned and sequenced and were correlated with the predicated alternative splicing events, suggesting that the bioinformatic prediction based on the primary sequencing data was reliable.figure sequence mapping and verification of genes with exon skipping events detected by rna-seq by pcr and rt-pcr. the expression profiles of the same gene in female and male parasites were placed under the line representing the chromosome position. the black lines represent original annotated gene structures, while the active transcripts in red and blue identified from the same genes in female and male parasites are underneath. the five genes and transcripts with exon skipping events were confirmed by pcr and rt-pcr. ge indicates pcr products amplified from genomic dna, and ce indicates pcr products amplified from cdna. red arrows indicate transcripts generated by exon skipping, and green arrows indicate primer locations. sequence mapping and verification of three genes with intron retention events detected by rna-seq by pcr and rt-pcr. the three genes and transcripts with intron retention events were confirmed by pcr and rt-pcr. gi indicates pcr products amplified from genomic dna, and ci indicates pcr products amplified from cdna. red arrows indicate transcripts generated by intron retention, and green arrows indicate primer locations. functional category of alternatively spliced genes in s. japonicum after mapping of the rna-seq reads to the s. japonicum reference genome, transcripts were assembled and their relative abundances were calculated. cuffdiff was used to find significant changes in gene level expression in the two libraries. genes subject to alternative splicing and showing significant differences in expression in the two libraries are listed in table  and additional file 6: figure s genes related to the function of genetic information processing were found to be more biased to the female parasite while genes related to the environmental information processing were more active in the male parasites. these data reflect the biology of the two sexes of the parasite. the female parasites, which are kept in the cavity of the male parasites, are more active in the reproduction process while the male parasites are principally responsible in the host–parasite interaction.table go classification statistics of alternatively processed genes that were differentially expressed in female and male s. japonicum as alternative splicing. diff differentially expressed. total as & diff, genes with both alternative splicing and differential expression in female and male s. japonicum. identification of novel transcripts from intergenic regions and previously determined pseudogenes one of the advantages of transcriptomic analysis is that it allows identification of novel transcripts that may not be predicted based on genomic sequences. the novel sequences can thus provide a powerful tool for re-annotation of the genome of s. japonicum, which has been poorly assembled. of the novel polyadenylated sequences, two classes of transcripts have been identified: one that does not map either to regions of the genome corresponding to annotated genes or to the untranslated regions, and another that maps to previously annotated pseudogenes. we identified novel transcripts that completely matched the previously annotated intergenic regions of the genome. the length of these transcripts was from to  bp, with an average length of  bp. it has been reported that the s. mansoni genome contains many small open reading frames. our results indicated that the small transcripts derived from both intergenic and "pseudogenes" in s. japonicum may encode important functions, as reported for the human genome. further, 31% of these transcripts had at least one complete open reading frame that could be translated into proteins; the other 69% were not annotated. thus, of the novel transcripts, at least genomic sequences corresponding to the transcripts can be re-annotated as protein-coding genes. a total of transcripts were mapped to the non-coding rna database of the rfam, which is also frequently used as a source of high-quality alignments for training and benchmarking rna sequence analysis software tools. these transcripts were found as either micrornas or ribosomal or other non-coding rnas. they were likely the contaminated sequences which were not completely depleted during mrna purification process.figure numbers of non-coding rna transcripts identified in the transcriptomes of female and male s. japonicum. furthermore, among the novel transcripts, we detected that were derived from pseudogenes; of these, were derived from annotated pseudogenes, and the rest were from unannotated pseudogenes. pseudogenes can be transcribed from either direction, which contributed to the templates for generation of small endogenous interfering rnas in s. japonicum, which is more common in transposable elements. the identification of these pseudogene-derived transcripts suggested that all sequences are polyadenylated and that there is no discrimination between coding and non-coding transcripts in the rna polyadenylation process in schistosoma parasites. on the other hand, complete reading frames were indeed identified in a number of pseudogene-derived transcripts that encoded proteins with known functions. thus, these "pseudogene" genes can be re-annotated as protein-coding genes.figure go categories of transcripts with a complete open reading frame and derived from pseudogenes. the percentages of the transcripts encoding the same proteins with similar function are indicated on the left while the numbers of the transcripts identified are indicated on the right.  in summary, by using rna-seq technology, we obtained the global transcriptomes of male and female s. japonicum. approximately 80% of the total reference genes were expressed in the adult stage of the parasite, representing the majority of the transcriptomes. these results further provide a comprehensive view of the global transcriptome of s. japonicum. the findings of a substantial level of alternative splicing events dynamically occurring in the parasitization in the mammalian hosts of the s. japonicum suggest complicated transcriptional and post-transcriptional regulatory mechanisms employed by the parasite. the data should not only significantly improve the re-annotation of the genome sequences but also should provide new information about the biology of the parasite. electronic supplementary material additional file 1: differentially expressed genes between females and males. additional file 2: figure s1: junction sites. additional file 3: table s1: primers and sequences for verification of the alternative splicing events. additional file 4: total transcripts identified in female and male parsites. additional file 5: genes with alternative splicing in female and male parsites. additional file 6: figure s2: go category of the genes that were alternatively spliced and also differentially transcribed. additional file 7: genes related to the function of genetic information processing and environmental information. additional file 8: novel transcripts identified. additional file 9: noncoding transcripts identified in the female s. japonicum. additional file 10: noncoding transcripts identified in the male s. japonicum. competing interests the authors declare that they have no competing interests. authors’ contributions xp, nh and qc conceived and designed the experiments. xp, nh, pc, sl, and cw performed the experiments. xp, nh and qc analyzed the data. xp and qc wrote the manuscript. all authors read and approved the final manuscript. we appreciate very much the bioinformatic support of dr. haibo sun at miningene biotechnology and the efforts of technicians at shenzhen bgi for solexa sequencing. we also thank the schistosoma japonicum genome sequencing and functional analysis consortium for making the s. japonicum genome available on line publicly. funding this study was supported by the national natural science foundation of china, the intramural grant from institute of pathogen biology, cams, the national s & t major program and the program for changjiang scholars and innovative research team in university.
in this study, the polyadenylated transcripts from the male and female s. japonicum were sequenced using a high-throughput rna-seq technique. by using rna-seq technology, we obtained the global transcriptomes of male and female s. japonicum. further, we identified four types of post-transcriptional processing, or alternative splicing, in both female and male worms of s. japonicum: exon skipping, intron retention, and alternative donor and acceptor sites. the findings of a substantial level of alternative splicing events dynamically occurring in s. japonicum parasitization of mammalian hosts suggest complicated transcriptional and post-transcriptional regulation mechanisms employed by the parasite. these results further provide a comprehensive view of the global transcriptome of s. japonicum. schistosoma japonicum is a pathogen of the phylum platyhelminthes that causes zoonotic schistosomiasis in china and southeast asian countries where a lack of efficient measures has hampered disease control. the numbers of protein-coding sequences identified in the transcriptomes of the female and male s. japonicum were and respectively, which is more than predicted from the annotated genome sequence. unlike mammalian organisms, in s. japonicum, the alternative donor and acceptor sites were more common than the other two types of post-transcriptional processing. in total, respectively and alternative splicing events were predicted in the transcriptomes of female and male s. japonicum. these data should not only significantly improve the re-annotation of the genome sequences but also should provide new information about the biology of the parasite. bioinformatic and experimental analyses focused on post-transcriptional rna processing, which revealed extensive alternative splicing events in the adult stage of the parasite.
 schistosoma japonicum is a pathogen of the phylum platyhelminthes that causes zoonotic schistosomiasis in china and southeast asian countries where a lack of efficient measures has hampered disease control. the development of tools for diagnosis of acute and chronic infection and for novel antiparasite reagents relies on understanding the biological mechanisms that the parasite exploits.  in this study, the polyadenylated transcripts from the male and female s. japonicum were sequenced using a high-throughput rna-seq technique. bioinformatic and experimental analyses focused on post-transcriptional rna processing, which revealed extensive alternative splicing events in the adult stage of the parasite. the numbers of protein-coding sequences identified in the transcriptomes of the female and male s. japonicum were and respectively, which is more than predicted from the annotated genome sequence. further, we identified four types of post-transcriptional processing, or alternative splicing, in both female and male worms of s. japonicum: exon skipping, intron retention, and alternative donor and acceptor sites. unlike mammalian organisms, in s. japonicum, the alternative donor and acceptor sites were more common than the other two types of post-transcriptional processing. in total, respectively and alternative splicing events were predicted in the transcriptomes of female and male s. japonicum.  by using rna-seq technology, we obtained the global transcriptomes of male and female s. japonicum. these results further provide a comprehensive view of the global transcriptome of s. japonicum. the findings of a substantial level of alternative splicing events dynamically occurring in s. japonicum parasitization of mammalian hosts suggest complicated transcriptional and post-transcriptional regulation mechanisms employed by the parasite. these data should not only significantly improve the re-annotation of the genome sequences but also should provide new information about the biology of the parasite. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. issue-copyright-statement© the author 2014 