Genome-wide sequencing of small RNAs reveals a tissue-specific loss of conserved microRNA families in Echinococcus granulosus
© Bai et al.; licensee BioMed Central Ltd. 2014
Received: 20 March 2014
Accepted: 20 August 2014
Published: 29 August 2014
MicroRNAs (miRNAs) are important post-transcriptional regulators which control growth and development in eukaryotes. The cestode Echinococcus granulosus has a complex life-cycle involving different development stages but the mechanisms underpinning this development, including the involvement of miRNAs, remain unknown.
Using Illumina next generation sequencing technology, we sequenced at the genome-wide level three small RNA populations from the adult, protoscolex and cyst membrane of E. granulosus. A total of 94 pre-miRNA candidates (coding 91 mature miRNAs and 39 miRNA stars) were in silico predicted. Through comparison of expression profiles, we found 42 mature miRNAs and 23 miRNA stars expressed with different patterns in the three life stages examined. Furthermore, considering both the previously reported and newly predicted miRNAs, 25 conserved miRNAs families were identified in the E. granulosus genome. Comparing the presence or absence of these miRNA families with the free-living Schmidtea mediterranea, we found 13 conserved miRNAs are lost in E. granulosus, most of which are tissue-specific and involved in the development of ciliated cells, the gut and sensory organs. Finally, GO enrichment analysis of the differentially expressed miRNAs and their potential targets indicated that they may be involved in bi-directional development, nutrient metabolism and nervous system development in E. granulosus.
This study has, for the first time, provided a comprehensive description of the different expression patterns of miRNAs in three distinct life cycle stages of E. granulosus. The analysis supports earlier suggestions that the loss of miRNAs in the Platyhelminths might be related to morphological simplification. These results may help in the exploration of the mechanism of interaction between this parasitic worm and its definitive and intermediate hosts, providing information that can be used to develop new interventions and therapeutics for the control of cystic echinococcosis.
KeywordsEchinococcus granulosus microRNA Deep sequencing Differential expression Life cycle stage development
MicroRNAs (miRNAs) are a class of endogenous non-coding RNAs of around 22 nucleotides (nts) in length. They play a role as post-transcriptional regulators, partially or completely complementary binding to messenger RNA transcripts (mRNAs), usually resulting in direct degradation or translational repression of target genes [1, 2]. In eukaryotes, miRNAs are involved in a broad variety of biological processes , such as embryonic development, cell proliferation, cell differentiation, and apoptosis. To date, more than twenty thousand miRNAs from 223 species have been registered in the miRBase (http://www.mirbase.org/, release 21.0, June 2014) . miRNAs are usually highly conserved throughout the animal kingdom . They have been continuously added to Metazoan genomes. The emergence of vertebrates is characterized by a strong increase in miRNA families, and correlates with the increase in vertebrate morphological complexity [5–7]. Therefore, miRNAs may have significantly contributed to phenotypic evolution in animals. Recently, a high rate of the loss of conserved miRNA loci has been found in the nematodes and the flatworms [8, 9]. This result suggested that the miRNA losses in Platyhelminths are not random events. However, to date, there is no report about the mechanisms of miRNA losses in metazoan phylogenetics.
Echinococcus granulosus is a flatworm and member of the Platyhelminthes, and is the causative agent of cystic echinococcosis (CE), a disease that is distributed throughout most areas of the world [10, 11]. Currently, up to 3 million people are infected with E. granulosus , and, in some areas, 10% of the population has detectable hydatid cysts by abdominal ultrasound and chest X-ray [13, 14]. This tapeworm requires two mammalian hosts to complete its life cycle. The mature adult worm resides in the small intestine of a carnivore (definitive host) and releases worm segments or proglottides containing hundreds of eggs which contaminate vegetation and water. When the eggs are swallowed by an intermediate host such as a sheep, the eggs hatch to release larval oncospheres into the digestive tract which are activated by bile and gastrointestinal enzymes. The activated oncospheres penetrate the intestinal wall and migrate via the circulatory system to various organs (mainly the liver and lungs). In these organs, the oncospheres develop into hydatid cysts over many months and the cysts generate brood capsules within which protoscoleces (PSC) are produced asexually. The cycle is completed when a canine (dog, wolf, fox) definitive host swallows PSC present in the hydatid cysts in infected offal, and the PSC develop into adult worms in the small intestine .
Although the life cycle of E. granulosus is well known, the mechanisms underlying the main developmental events throughout remain largely unclear. Our previous study of the transcriptome of E. granulosus revealed that 1,452 genes were up- or down-regulated in adult, PSC and cyst stages . Moreover, a global proteomic analysis of the expression characteristics of E. granulosus in larval and adult stages identified 22 adult-specific and 263 PSC-specific proteins . These studies suggested that transcriptional regulatory mechanisms are pivotal in the control of E. granulosus development.
In regards to flatworms, miRNAs have been experimentally identified in Schmidtea mediterranea [18, 19], Schistosoma japonicum [20–22] and Schistosoma mansoni [23, 24]. A PCR-based cloning study identified 26 mature miRNAs in PSC and the cystic stage of E. granulosus . However, due to the hitherto restricted experimental methods and the limited genomic information available, the numbers and expression profile of E. granulosus miRNAs are still unclear. Here, we used next generation sequencing technology (NGS) to further explore the diversity of E. granulosus miRNAs and their expression patterns in different life stages. We expand the miRNA repertoire of E. granulosus and identify new miRNA encoding loci. Through comparing miRNA families in the Platyhelminths, we found that the losses of E. granulosus miRNAs may be associated with the loss of ciliated cells, the gut and sensory organs. The results significantly enhance our knowledge of miRNA species in E. granulosus and provide insights into miRNA evolution, biogenesis, and expression in parasites generally.
Deep sequencing of three small RNA libraries from E. granulosus
Identification of known and novel miRNAs from E. granulosus
To date, 23 mature miRNAs of E. granulosus have been identified  and are included in the miRBase database 20.0 (http://www.mirbase.org/). By deep sequencing, we found that all of the known mature miRNAs were present in our data sets [Additional file 1: Table S3], the majority being abundant in all three libraries. Furthermore, we also identified 23 miRNA stars from the known miRNA precursors [Additional file 1: Table S3].
In addition to known miRNAs, we also used miRDeep2 to predict and score novel miRNA precursors  and identified 94 miRNA candidates encoding 91 mature miRNAs and 39 miRNA stars [Additional file 1: Table S4 and S5]. All these miRNAs can be folded into characteristic miRNA stem-loop secondary hairpin structures and have a 1-2 nt 3′overhang pattern generated by Dicer cleavage during mature miRNA generation [Additional file 2]. We evaluated evolutionary conservations by homologous searches to known metazoan miRNAs and found 11 pre-miRNAs were classified into known families based on their precursor sequences, whereas 83 did not show homology with other miRNAs. We further matched these predicted pre-miRNA candidates to the E. multilocularis genome (http://www.sanger.ac.uk/cgi-bin/blast/submitblast/Echinococcus) and found 82 of 94 miRNA candidates were evolutionarily conserved (identity ≥87%) in both species [Additional file 1: Table S6]. To validate the novel miRNAs, we randomly selected 22 mature miRNAs and 5 miRNA stars, and conducted stem-loop RT-PCR . All the selected miRNAs were expressed in E. granulosus [Additional file 3: Figure S1], suggesting that the filter criteria were sufficiently strict for predicting novel miRNAs.
Comparison of miRNA families of E. granulosuswith other flatworms
We then compared the presence or absence of the conserved miRNA families of E. granulosus (Cestoda) with three other flatworms, S. mediterranea (Turbellaria) , Gyrodactylus salaris (Monogenea)  and S. japonicum (Trematoda) . Similar to the previous study , 8 conserved bilaterian, protostomian, and lophotrochozoan miRNA families were not detected in the four flatworms (Figure 3b and c); these were mir-33, mir-34, mir-137, mir-193, mir-242, mir-252, mir-137 and mir-375. Furthermore, 8 miRNAs families, mir-29, mir-92, mir-216, mir-210, mir-315, mir-12, mir-750 and mir-1989 were only detected in the free living flatworm S. mediterranea, but not in G. salaris, E. granulosus and S. japonicum. Finally, there are 6 miRNA families (mir-76, mir-96, mir-184, mir-278, mir-2001 and mir-1175) that are present in S. mediterranea and G. salaris, but are absent in S. japonicum and E. granulosus.
Expression profiles of miRNAs in the three life stages of E. granulosus
The relative abundance of the top ten mature miRNAs in three life-stages of E. granulosus
Since the larval PSC is an important transition life cycle stage, capable of developing either into an adult worm in the dog gastrointestinal tract or a secondary hydatid cyst in the intermediate host, we compared the expression levels of miRNAs in adult worms or cyst membrane with those in the PSC. A total of 43 miRNAs (22 up-regulated and 21 down-regulated) in adult tissue and 42 miRNAs (4 up-regulated and 38 down-regulated) in the cyst membrane were identified (Figure 4b). Further removing similar trends of miRNA expression between adults vs PSC and cyst vs PSC, we found that 25 miRNAs (22 up-regulated and 3 down-regulated) and 24 miRNAs (4 up-regulated and 20 down-regulated) were specifically changed in the adult and cyst stage, respectively (Figure 4c). Amongst these, 8 miRNAs appeared to show reverse trends in expression in the adult and cyst stages compared with PSC; miR-125, miR-277, miR-4989*, new-17*, new-87 and new-87* were up-regulated in adult worms but down-regulated in cyst membrane. In contrast, miR-124a and miR-124b were down-regulated in the former stage but up-regulated in the latter [Additional file 1: Table S11]. These results imply that miRNA expression variation in the different life stages may be associated with the direction of development of PSC into either an adult worm or a secondary hydatid cyst.
Target gene prediction and functional analysis
To assign biological functions to the 65 differentially expressed miRNAs, we predicted putative target genes using the Miranda program . Based on previously published EST sequences , we found 3,622 genes possibly targeted by 114 mature miRNAs and 62 miRNA stars [Additional file 1: Table S12]. By enrichment analyses on the predicted targets, we determined that 182 genes could be categorized into 24 significant GO terms (adjust P < 0.05) [Additional file 1: Table S13]. The major target genes in the biological process categories were involved in the regulation of cell differentiation (GO:0045595), determination of adult lifespan (GO:0008340) and response to nutrient (GO:0007584). These results support our premise that the differentially expressed miRNAs may have important functions in the three different life cycle stages of E. granulosus.
In the current study, we performed deep sequencing of small RNAs in three distinct developmental stages of E. granulosus. By Bowtie mapping, there are 11,680,028, 12,966,593 and 9,375,095 18–30 nt reads were located the E. granulosus genome, respectively. The remaining 40% unmapped reads could be ascribed to the following reasons: (1) some genetic polymorphisms may occur in the E. granulosus genome which interferes with Bowtie mapping; (2) unsatisfactory available genome data; some sequences are probably absent in the current draft genome of E. granulosus ; (3) miRNAs sequencing errors. NGS may randomly produce some base calling errors and influence the accuracy of short reads; (4) post-transcriptional editing of miRNAs may also have contributed to the mismatches observed.
It has been previously reported that the gain of miRNA genes in Metazoan species is associated with an increase in morphological complexity . However, a high rate of the loss of conserved miRNA loci has been found in many subgroups of the Platyhelminthes . Basing on evolutionary conservation analysis, we also found a significant loss of conserved miRNAs in the Cestoda. A total of 22 conserved miRNA families were confirmed as not detectable in E. granulosus, suggesting that the loss of conserved miRNA loci in flatworms is not a random event in metazoan phylogenetics. It might relate to loss of targets in the corresponding species  or reduced morphological complexity . A previous study revealed that miRNA invention is closely related to the evolution of tissue identities in bilaterian species . Through comparing these lost miRNA families in E. granulosus, we found 8 are not detected in the parasitic Platyhelminthes compared with S. mediterranea and 6 are absent from S. japonicum and E. granulosus. On further analysis of their evolutionary origin, we observed that most of the lost miRNAs may be derived from some specific tissues, which have been reduced or lost in tapeworms. For example, there are 8 conserved miRNA families expressed in the free-living S. mediterranea that are absent in the Neodermata. Among these, miR-12, miR-216 and miR-283 have been found mainly present in locomotor ciliated cells, while miR-29 and miR-92 are expressed in the gut . These organs have gradually disappeared in E. granulosus and some of the other parasitic Platyhelminthes groups because of a significant decrease in locomotor and digestive activities associated with the parasitic mode of life. Furthermore, in the Neodermata, there is a significant degeneration of peripheral sensory nervous system elements, such as, the loss of eyes in tapeworms. In bilaterians, miR-183 and miR-263 have been reported to have a conserved affiliation with sensory organ differentiation [7, 39] and miR-184 and miR-278 play essential roles in the development of eyes [40, 41]. These sensory tissue-specific miRNAs are not detected in S. japonicum (Trematoda) and E. granulosus (Cestoda). Thus, it seems that the loss of miRNA loci in E. granulosus may be related to the reduction or loss of some specific tissues in the Cestoda.
Additionally, the biological evolution of one species to another is usually accompanied by the acquisition of new miRNA families in the newly evolved species . Although previous studies did not report there was any miRNA family acquired in the Platyhelminthes [8, 9], we observed that mir-3479 was present in E. granulosus, S. mansoni and S. japonicum. This result may be due to the lack of information in other taxa or due to evolutionary innovations to adapt to the new living environment. For tapeworms, in order to escape from host immune attack, E. granulosus produces many components to moderate the host immune response, such as EgAgB (antigen B), which interferes with human dendritic cell differentiation and maturation . Further, one recent study has shown that three schistosome-specific miRNAs, Bantam, miR-10 and miR-3479, were detected in the plasma of S. japonicum-infected rabbits . These results imply that the gain of miR-3479 in E. granulosus or S. japonicum may have resulted as an evolved adaptation in the parasitic Platyhelminthes to moderate the host immune response. It is important to note that there are no reports to date on the plasma levels of E. granulosus miRNAs in patients with CE. However, some Dicrocoelium-specific miRNAs have recently been shown to be present in exosome-like vesicles  and similar secretory vesicles have been described previously in the early stages of E. granulosus development .
Notably, two novel close miRNA clusters (new-15/ new-24/new-61/new-7 and new-12/new-22) were identified in E. granulosus and E. multilocularis. Beside the same genomic location, they display similar patterns of miRNA expression, suggesting that the expression of these novel miRNAs may be co-regulated by the same promoter elements. Furthermore, two clusters of miRNAs share the ‘GGCGCUU’ motif with the known sme-miR-2160. Recently, it has been reported that sme-miR-2160 is enriched in planarian stem cells and is down-regulated in individuals in which stem cells had been abrogated by irradiation . These observations imply that these two clusters of miRNAs may be associated in E. granulosus with cell proliferation and differentiation. Moreover, since new-7/12/15/22/24/61 of E. granulosus contains six copies of the same seed sequences in the 5′ hairpin arms, these multiple copies may likely be functional, resulting in the dose-dependent regulation of their downstream targets [47, 48]. Similar multiple copies of mir-2160 have been identified in S. mediterranea .
Development in E. granulosus involves a number of different life stages including the eggs, activated oncospheres, hydatid cyst, PSC and adult worms. It is a complicated and dynamic process and may be associated with miRNA regulation. By comparing the differential expression pattern of miRNAs, our results showed that let-7 exhibited 3.2-fold and 2.6-fold increased expression in PSC and in the cyst, respectively, compared with the adult. let-7 has been shown to regulate developmental timing in C.elegans through the direct target genes lin-41 and hbl-1 . Thus, the observation of substantially increased expression of let-7 in the PSC and cyst might be associated with the capability of E. granulosus for bi-directional development, differentiating into an adult worm in the definitive host dog gastrointestinal tract or into a secondary hydatid cyst in the intermediate (or human) host .
Indeed, bi-directional development is a unique biological characteristic of E. granulosus and our previous study showed that some protein-coding genes, with high homology to the bile acid receptors and transporters, were located in the E. granulosus genome, and that nuclear hormone receptors, including the farnesoid X receptor (FXR) and VDR, may play a key role in stimulating the PSC to develop into adult worms . Recently, the hormone receptors have been reported to directly activate transcription of the evolutionarily conserved let-7 and control the development of larvae to adults in C. elegans . Here, we observed a positive correlation between the expression of let-7 and VDR, suggesting that VDR may induce crucial miRNA transcriptional events and promote adult development in E. granulosus. Furthermore, through miRNA target prediction analysis, we found let-7 could complementary bind to the VDR 3′ UTR sequences, suggesting that let-7 may mediate a negative feedback loop controlling the bile acid signal pathway in the development of PSC. Recently, a feedback circuit between let-7-family miRNAs and DAF-12, a homolog of vertebrate FXR and VDR, has been found in C. elegans . Additional experimental validations are required to explore this relationship in E. granulosus in future studies.
In summary, this study describes the first large scale identification and dynamic characterisation of miRNAs in three developmental stages of E. granulosus. A total of 114 mature miRNAs and 62 miRNA stars were identified. Evolutionary conservation analysis suggested that the losses of miRNA families in E. granulosus may be associated with morphological reductions. Additionally, GO analysis revealed that the differentially expressed miRNAs and their targets may be involved in diverse development processes, including bi-directional differentiation and nervous system development. Understanding the regulatory processes involving miRNAs in E. granulosus may be helpful to explore the mechanism of interaction between this parasitic worm and its definitive and intermediate hosts, and provide new information to develop new interventions and therapeutics for the control of CE.
Ethical approval for the use of dogs and sheep was given by the Animal Ethics Committee of Xinjiang Academy of Animal Science(XJAASaec2011031506). The protocols for using dogs and sheep were covered under the < The Guidelines for the Care and Use of Experimental Animals > (China) National Science and Development Act 2006–398.
Preparation of E. granulosussamples
The preparation of E. granulosus samples for analysis was undertaken as previously described [53, 54]. PSC were aspirated aseptically from several sheep liver hydatid cysts collected from Xinjiang, China. The PSCs were sedimented at 1000 × g for 15 minutes and digested by pepsin . After PSC were removed, cyst membranes were rinsed 10 times with PBS and checked microscopy to confirm no PSC attached on the membranes, and then cut into small pieces with a pair of scissors. The cut membranes were soaked in PBS in a flask and stirred for 10 min at 4°C. The flask was placed on bench at room temperature for 3 min to precipitate laminated layer and other cyst tissues. The membrane cells in supernatant were pelleted by centrifuging at 1400 g for 10 min at 4°C. The cells were resuspended with PBS and left on the bench for 3 min to precipitate remained PSC, the supernatant was centrifuged at 1400 × g for 10 min at 4°C. Adult worms were collected from experimentally infected dogs using procedures describe earlier . All parasite materials were subjected to a final wash in PBS and then they were suspended in 10 volumes of RNAlater (Ambion, Austin, TX, USA), and stored at -80°C.
Total RNA extraction
Samples of the 3 different life cycle stages of E. granulosus (adult, cyst and PSC) were homogenized in lysis buffer and total RNA was extracted using mirVana™ miRNA isolation kits (Ambion, Austin, TX, USA), according to the manufacturer’s instructions. RNA concentration and purity were monitored using a NanoDrop ND-1000 UV spectrophotometer (Nanodrop Technologies, Wilmington, DE). The RNA integrity was evaluated by an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA) [Additional file 6: Figure S4]. Finally, all RNA samples were stored at -80°C.
Small RNA library construction and deep sequencing
We constructed each small RNA library using Illumina Small RNA v1.5 Sample Preparation Kits (Illumina, San Diego, CA). Briefly, total RNA (5 μg) of each sample was resolved on denaturing polyacrylamide gel electrophoresis (PAGE) gels, and then fractions of 18 to 30 nt in size were collected. The isolated small RNAs were sequentially ligated to 3′adapters, using T4 RNA ligase 2 (New England Biolabs, Ipswich, MA, USA), and 5′adapters using T4 RNA ligase (New England Biolabs, Ipswich, MA, USA). The ligation products were reverse transcribed using SuperScript Reverse Transcriptase II (Life Technologies, Gaithersburg, MD, USA) and amplified with 12 PCR cycles. Then, 6% (w/v) PAGE was used to purify the amplification products. Finally, the libraries were used for clustering and sequencing using an Illumina Genome Analyzer II (Illumina, San Diego, CA).
Computational analysis to search for novel miRNAs and other small RNAs
After removing low quality sequence reads and trimming adaptor sequences through miRExpress , we collected small RNAs ranging from 18–30 nucleotides and count length distributions. Then, duplicated sequences were removed from the initial dataset to produce a nonredundant set of unique sequences, hereafter referred to as sequence tags. To determine whether these small RNA sequences were from E. granulosus, we mapped these small RNA sequences by Bowtie  to the draft E. granulosus genome sequences . Only those mapping perfectly onto the draft genome were further considered as candidate miRNAs. If a sequence mapped to more than 10 loci in the genome, it was not considered for further analysis. Then, we compared all sequence tags against a database of known miRNA precursors sequences (http://www.mirbase.org/, release 21.0) , and profiled every annotated miRNA in each library.
Basing on our earlier findings  and scrutiny of the Rfam RNA family database (http://rfam.sanger.ac.uk), we discarded all sequences matching with known E. granulosus rRNAs, tRNAs, snRNAs and mRNAs. With the annotation of repeat sequences in the E. granulosus genome, small RNAs positioned at repeat loci were identified and annotated as repeat-associated small RNAs. Finally, the remaining clean sequenced reads were used to search for both conserved and novel miRNAs.
In order to identify new miRNAs, we used miRDeep 2.0  to identify miRNAs without prior information. The software uses a Bayesian algorithm and provides a score based on the expected pattern of miRNA-duplex excision from the stem of the precursor hairpins by Dicer . The secondary structures of putative precursors were identified using the RNAfold program  and computed by minimum free energy (MFE). We used the following criteria to filter the candidate miRNAs: (1) The lengths of mature miRNA ranged from 18 to 25 nucleotides; (2) The miRNA precursors should have a characteristic fold-back structure and MFE should be less than -20 Kcal/mol; (3) The miRDeep score should be more than 1; (4) The GC content of the mature miRNA should be 15%-85%; (5) Small RNAs with multiple loci in the E. granulosus genome were excluded (more than five times); (6) In order to remove accidental hits, only the precursors with at least 5 sequence reads in three libraries were used. Furthermore, these predicted mature miRNAs were also compared with known miRNAs from other species. If the E. granulosus miRNAs shared >80% homology with known miRNAs and had the same seed sequences (2–8 nt), the miRNAs were considered as conserved miRNAs and were thus named after the known miRNA . We defined the mature miRNA and miRNA star basing on the total number of reads in three life stages.
Detection of differential expression
Using the quantifier module of miRDeep2, we summed up read counts for all known or novel mature miRNAs in the sequencing data. Then, relative expression levels were calculated by counting the numbers of respective miRNA reads normalized to the total number of annotated miRNAs in each library. After ignoring the simultaneously lowly expressed miRNAs in three libraries (less than 50 TPM at all stages), differential expression in the three libraries was evaluated using the statistical R package DEGseq . (http://www.bioconductor.org/packages/2.6/bioc/html/DEGseq.html). Based on this statistical model, statistical significances were calculated using Fisher’s exact test and likelihood ratio test, and then corrected for multiple testing according to Benjamini-Hochberg multiple test method . Potentially interesting miRNA candidates were chosen according to the criteria of a 2-fold expression level change and adjusted P-value < 0.001.
MicroRNA target prediction
In order to seek miRNAs targets, we first predicted the 3′UTRs by combined analysis of the predicted mRNAs and the released EST scafford sequences for E. granulosus . Only the sequences located at the 3′untranslated region of predicted genes were regarded as potential 3′UTR sequences. Then, the miRanda program  was used to predict the target genes for all mature miRNAs. The energy thresholds were set at ≤ -20 kcal/mol and other thresholds used a default value (score threshold, 120; gap-opening penalty, -9; gap-extend penalty, -4). Then, all ESTs sequences predicted to contain miRNA target sites, were annotated as earlier described .
Gene ontology and pathway analysis
Target genes regulated by the different subgroups of miRNAs were collected and subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Based on NCBI-nr annotation, we used the BLAST2GO program (http://www.BLAST2go.org/) to obtain GO annotation of every target . Then, the hypergeometric test  was used to classify the GO category, and the false discovery rate (FDR) was calculated to correct the P-value. KEGG pathway analysis was carried out using the KEGG Automatic Annotation Server for ortholog assignment and pathway mapping (http://www.genome.jp/tools/kaas) . The hypergeometric test was used to assess the significant pathway of enrichment.
Quantitative RT-PCR of miRNAs and their targets
A stem-loop quantitative RT-PCR was performed to validate the miRNA expression levels in each developmental stage. Total RNA was extracted using mirVanaTM miRNA isolation kits (Ambion, Austin, TX, USA) according to the manufacturer’s instructions. The total RNA was reverse transcribed in a 20 μl reverse transcription (RT) reaction using 500 ng total RNA, 1 U PrimeScript RT Enzyme Mix I (Takara Bio, Otsu, Japan), 4 μl of 5 × RT buffer, and 1 μl of 2 μM miRNA-specific stem-loop primers [Additional file 1: Table S14a] designed according to Chen et al. . 5.8S rRNA was chosen as a reference and was reverse-transcribed with a specific reverse primer [Additional file 1: Table S14a]. The RT reaction mixture was incubated at 37°C for 30 minutes, then at 85°C for 5 seconds. A control was set up at the same time with no RNA input. With the cDNA products as a template, quantitative PCR was carried out using SYBR® Premix Ex TaqTM (Takara Bio, Otsu, Japan) in a StepOne Plus real-time system (Applied Biosystems, Carlsbad, CA, USA). The primer sequences are listed in [Additional file 1: Table S14a and S14b]. A 10 μl reaction mixture including 1 μl of diluted cDNA (1:5), 0.4 μl of each primer (10 mM), 5 μl SYBR Premix Ex Taq II and 3.6 ml H2O was placed in 0.2 ml eight-strip PCR tubes (Applied Biosystems, Carlsbad, CA, USA). Cycling conditions were: 95°C for 30s, followed by 40 cycles of 95°C for 5 s and 60°C for 34 s, 72°C for 60s. For each PCR, dissociation curve analysis was carried out to discriminate the specific products from the primer dimers. The CT-values are the average of three technical and three biological replicates and fold changes of miRNAs and their targets in different samples were calculated by the 2-ΔΔCt method.
The miRNA expression levels in each developmental stage were compared by using a two-tailed t-test. Correlation analysis was made using the Pearson correlation coefficient. P < 0.05 was considered statistically significant.
- E. granulosus :
- E. multilocularis :
- ATP2B3 :
Plasma membrane calcium-transporting atpase 3-like
- LCFACS :
Long-chain acyl-CoA synthetase
Minimum free energy
False discovery rate
Transcripts per million
Kyoto Encyclopedia of Genes and Genomes.
This work was supported by the National Natural Science Foundation of China (81271868, 31260272, and 31260272), National Basic Research Program (973) of China (2010CB534906), Shanghai Municipal Commission for Science and Technology (10XD1403200, 11DZ2292600, and 13DZ2291800), and Shanghai Natural Science Foundation (13ZR1429500).
High throughput sequencing data have been submitted to the NCBI Sequence Read Archive (SRA) under accession number SRP034671. The read-count files for each library have been submitted to NCBI Gene expression omnibus (GEO) under accession number GSE54180.
- Carrington JC, Ambros V: Role of microRNAs in plant and animal development. Science. 2003, 301 (5631): 336-338.PubMedView ArticleGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297.PubMedView ArticleGoogle Scholar
- Ambros V: MicroRNA pathways in flies and worms: growth, death, fat, stress, and timing. Cell. 2003, 113 (6): 673-676.PubMedView ArticleGoogle Scholar
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (Database issue): D152-157.PubMed CentralPubMedView ArticleGoogle Scholar
- Niwa R, Slack FJ: The evolution of animal microRNA function. Curr Opin Genet Dev. 2007, 17 (2): 145-150.PubMedView ArticleGoogle Scholar
- Heimberg AM, Sempere LF, Moy VN, Donoghue PC, Peterson KJ: MicroRNAs and the advent of vertebrate morphological complexity. Proc Natl Acad Sci U S A. 2008, 105 (8): 2946-2950.PubMed CentralPubMedView ArticleGoogle Scholar
- Christodoulou F, Raible F, Tomer R, Simakov O, Trachana K, Klaus S, Snyman H, Hannon GJ, Bork P, Arendt D: Ancient animal microRNAs and the evolution of tissue identity. Nature. 2010, 463 (7284): 1084-1088.PubMed CentralPubMedView ArticleGoogle Scholar
- Tarver JE, Sperling EA, Nailor A, Heimberg AM, Robinson JM, King BL, Pisani D, Donoghue PC, Peterson KJ: miRNAs: small genes with big potential in metazoan phylogenetics. Mol Biol Evol. 2013, 30 (11): 2369-2382.PubMedView ArticleGoogle Scholar
- Fromm B, Worren MM, Hahn C, Hovig E, Bachmann L: Substantial loss of conserved and gain of novel MicroRNA families in flatworms. Mol Biol Evol. 2013, 30 (12): 2619-2628.PubMed CentralPubMedView ArticleGoogle Scholar
- Eckert J, Deplazes P: Biological, epidemiological, and clinical aspects of echinococcosis, a zoonosis of increasing concern. Clin Microbiol Rev. 2004, 17 (1): 107-135.PubMed CentralPubMedView ArticleGoogle Scholar
- Craig PS, Larrieu E: Control of cystic echinococcosis/hydatidosis: 1863–2002. Adv Parasitol. 2006, 61: 443-508.PubMedView ArticleGoogle Scholar
- Craig PS, McManus DP, Lightowlers MW, Chabalgoity JA, Garcia HH, Gavidia CM, Gilman RH, Gonzalez AE, Lorca M, Naquira C, Nieto A, Schantz PM: Prevention and control of cystic echinococcosis. Lancet Infect Dis. 2007, 7 (6): 385-394.PubMedView ArticleGoogle Scholar
- Li T, Ito A, Pengcuo R, Sako Y, Chen X, Qiu D, Xiao N, Craig PS: Post-treatment follow-up study of abdominal cystic echinococcosis in tibetan communities of northwest Sichuan Province. China PLoS Negl Trop Dis. 2011, 5 (10): e1364-View ArticleGoogle Scholar
- Moro PL, Gilman RH, Verastegui M, Bern C, Silva B, Bonilla JJ: Human hydatidosis in the central Andes of Peru: evolution of the disease over 3 years. Clin Infect Dis. 1999, 29 (4): 807-812.PubMedView ArticleGoogle Scholar
- McManus DP, Zhang W, Li J, Bartley PB: Echinococcosis. Lancet. 2003, 362 (9392): 1295-1304.PubMedView ArticleGoogle Scholar
- Zheng H, Zhang W, Zhang L, Zhang Z, Li J, Lu G, Zhu Y, Wang Y, Huang Y, Liu J, Kang H, Chen J, Wang L, Chen A, Yu S, Gao Z, Jin L, Gu W, Wang Z, Zhao L, Shi B, Wen H, Lin R, Jones MK, Brejova B, Vinar T, Zhao G, McManus DP, Chen Z, Zhou Y, et al: The genome of the hydatid tapeworm Echinococcus granulosus. Nat Genet. 2013, 45 (10): 1168-1175.PubMedView ArticleGoogle Scholar
- Cui SJ, Xu LL, Zhang T, Xu M, Yao J, Fang CY, Feng Z, Yang PY, Hu W, Liu F: Proteomic characterization of larval and adult developmental stages in Echinococcus granulosus reveals novel insight into host-parasite interactions. J Proteomics. 2013, 84: 158-175.PubMedView ArticleGoogle Scholar
- Friedlander MR, Adamidi C, Han T, Lebedeva S, Isenbarger TA, Hirst M, Marra M, Nusbaum C, Lee WL, Jenkin JC, Sánchez Alvarado A, Kim JK, Rajewsky N: High-resolution profiling and discovery of planarian small RNAs. Proc Natl Acad Sci U S A. 2009, 106 (28): 11546-11551.PubMed CentralPubMedView ArticleGoogle Scholar
- Palakodeti D, Smielewska M, Graveley BR: MicroRNAs from the Planarian Schmidtea mediterranea: a model system for stem cell biology. RNA. 2006, 12 (9): 1640-1649.PubMed CentralPubMedView ArticleGoogle Scholar
- Hao L, Cai P, Jiang N, Wang H, Chen Q: Identification and characterization of microRNAs and endogenous siRNAs in Schistosoma japonicum. BMC Genomics. 2010, 11: 55-PubMed CentralPubMedView ArticleGoogle Scholar
- 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 CentralPubMedView ArticleGoogle Scholar
- Huang J, Hao P, Chen H, Hu W, Yan Q, Liu F, Han ZG: Genome-wide identification of Schistosoma japonicum microRNAs using a deep-sequencing approach. PLoS One. 2009, 4 (12): e8206-PubMed CentralPubMedView ArticleGoogle Scholar
- de Souza GM, Muniyappa MK, Carvalho SG, Guerra-Sa R, Spillane C: Genome-wide identification of novel microRNAs and their target genes in the human parasite Schistosoma mansoni. Genomics. 2011, 98 (2): 96-111.View ArticleGoogle Scholar
- Simoes MC, Lee J, Djikeng A, Cerqueira GC, Zerlotini A, da Silva-Pereira RA, Dalby AR, LoVerde P, El-Sayed NM, Oliveira G: Identification of Schistosoma mansoni microRNAs. BMC Genomics. 2011, 12: 47-PubMed CentralPubMedView ArticleGoogle Scholar
- Cucher M, Prada L, Mourglia-Ettlin G, Dematteis S, Camicia F, Asurmendi S, Rosenzvit M: Identification of Echinococcus granulosus microRNAs and their expression in different life cycle stages and parasite genotypes. Int J Parasitol. 2011, 41 (3–4): 439-448.PubMedView ArticleGoogle Scholar
- Henderson IR, Zhang X, Lu C, Johnson L, Meyers BC, Green PJ, Jacobsen SE: Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning. Nat Genet. 2006, 38 (6): 721-725.PubMedView ArticleGoogle Scholar
- Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N: miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012, 40 (1): 37-52.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33 (20): e179-PubMed CentralPubMedView ArticleGoogle Scholar
- Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel genes coding for small expressed RNAs. Science. 2001, 294 (5543): 853-858.PubMedView ArticleGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20.PubMedView ArticleGoogle Scholar
- Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115 (7): 787-798.PubMedView ArticleGoogle Scholar
- Meyers BC, Tej SS, Vu TH, Haudenschild CD, Agrawal V, Edberg SB, Ghazal H, Decola S: The use of MPSS for whole-genome transcriptional analysis in Arabidopsis. Genome Res. 2004, 14 (8): 1641-1653.PubMed CentralPubMedView ArticleGoogle Scholar
- Brenner S, Johnson M, Bridgham J, Golda G, Lloyd DH, Johnson D, Luo S, McCurdy S, Foy M, Ewan M, Roth R, George D, Eletr S, Albrecht G, Vermaas E, Williams SR, Moon K, Burcham T, Pallas M, DuBridge RB, Kirchner J, Fearon K, Mao J, Corcoran K: Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays. Nat Biotechnol. 2000, 18 (6): 630-634.PubMedView ArticleGoogle Scholar
- Wang Z, Xue X, Sun J, Luo R, Xu X, Jiang Y, Zhang Q, Pan W: An “in-depth” description of the small non-coding RNA population of Schistosoma japonicum schistosomulum. PLoS Negl Trop Dis. 2010, 4 (2): e596-PubMed CentralPubMedView ArticleGoogle Scholar
- Kato M, de Lencastre A, Pincus Z, Slack F: Dynamic expression of small non-coding RNAs, including novel microRNAs and piRNAs/21U-RNAs, during Caenorhabditis elegans development. Genome Biol. 2009, 10 (5): R54-PubMed CentralPubMedView ArticleGoogle Scholar
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol. 2003, 5 (1): R1-PubMed CentralPubMedView ArticleGoogle Scholar
- Sperling EA, Peterson KJ: microRNAs and metazoan phylogeny: big trees from little genes. Animal evolution-genomes, trees and fossils. Edited by: Telford MJ, Littlewood DTJ. 2009, Oxford: Oxford University Press, 157-170.View ArticleGoogle Scholar
- Erwin DH, Laflamme M, Tweedt SM, Sperling EA, Pisani D, Peterson KJ: The Cambrian conundrum: early divergence and later ecological success in the early history of animals. Science. 2011, 334 (6059): 1091-1097.PubMedView ArticleGoogle Scholar
- Aboobaker AA, Tomancak P, Patel N, Rubin GM, Lai EC: Drosophila microRNAs exhibit diverse spatial expression patterns during embryonic development. Proc Natl Acad Sci U S A. 2005, 102 (50): 18017-18022.PubMed CentralPubMedView ArticleGoogle Scholar
- Shalom-Feuerstein R, Serror L, De La Forest DS, Petit I, Aberdam E, Camargo L, Damour O, Vigouroux C, Solomon A, Gaggioli C, Itskovitz-Eldor J, Ahmad S, Aberdam D: Pluripotent stem cell model reveals essential roles for miR-450b-5p and miR-184 in embryonic corneal lineage specification. Stem Cells. 2012, 30 (5): 898-909.PubMedView ArticleGoogle Scholar
- Nairz K, Rottig C, Rintelen F, Zdobnov E, Moser M, Hafen E: Overgrowth caused by misexpression of a microRNA with dispensable wild-type function. Dev Biol. 2006, 291 (2): 314-324.PubMedView ArticleGoogle Scholar
- Wheeler BM, Heimberg AM, Moy VN, Sperling EA, Holstein TW, Heber S, Peterson KJ: The deep evolution of metazoan microRNAs. Evol Dev. 2009, 11 (1): 50-68.PubMedView ArticleGoogle Scholar
- Rigano R, Buttari B, Profumo E, Ortona E, Delunardo F, Margutti P, Mattei V, Teggi A, Sorice M, Siracusano A: Echinococcus granulosus antigen B impairs human dendritic cell differentiation and polarizes immature dendritic cell maturation towards a Th2 cell response. Infect Immun. 2007, 75 (4): 1667-1678.PubMed CentralPubMedView ArticleGoogle Scholar
- Cheng G, Luo R, Hu C, Cao J, Jin Y: Deep sequencing-based identification of pathogen-specific microRNAs in the plasma of rabbits infected with Schistosoma japonicum. Parasitology. 2013, 140 (14): 1751-1761.PubMedView ArticleGoogle Scholar
- Bernal D, Trelis M, Montaner S, Cantalapiedra F, Galiano A, Hackenberg M, Marcilla A: Surface analysis of Dicrocoelium dendriticum. The molecular characterization of exosomes reveals the presence of miRNAs. J Proteomics. 2014, 105: 232-241.PubMedView ArticleGoogle Scholar
- Holcman B, Heath DD: The early stages of Echinococcus granulosus development. Acta Trop. 1997, 64 (1–2): 5-17.PubMedView ArticleGoogle Scholar
- Giraldez AJ, Cinalli RM, Glasner ME, Enright AJ, Thomson JM, Baskerville S, Hammond SM, Bartel DP, Schier AF: MicroRNAs regulate brain morphogenesis in zebrafish. Science. 2005, 308 (5723): 833-838.PubMedView ArticleGoogle Scholar
- Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, Van Dongen S, Inoue K, Enright AJ, Schier AF: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006, 312 (5770): 75-79.PubMedView ArticleGoogle Scholar
- Hunter SE, Finnegan EF, Zisoulis DG, Lovci MT, Melnik-Martinez KV, Yeo GW, Pasquinelli AE: Functional genomic analysis of the let-7 regulatory network in Caenorhabditis elegans. PLoS Genet. 2013, 9 (3): e1003353-PubMed CentralPubMedView ArticleGoogle Scholar
- Smyth JD: In Vitro Cultivation Of Parasitic Helminths. 1990, Boca Raton: CRC PressGoogle Scholar
- Bethke A, Fielenbach N, Wang Z, Mangelsdorf DJ, Antebi A: Nuclear hormone receptor regulation of microRNAs controls developmental progression. Science. 2009, 324 (5923): 95-98.PubMed CentralPubMedView ArticleGoogle Scholar
- Hammell CM, Karp X, Ambros V: A feedback circuit involving let-7-family miRNAs and DAF-12 integrates environmental signals and developmental timing in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 2009, 106 (44): 18668-18673.PubMed CentralPubMedView ArticleGoogle Scholar
- Li J, Zhang CS, Lu GD, Wang JH, Wen H, Yan GQ, Wei XF, Lin RY: Molecular characterization of a signal-regulated kinase homolog from Echinococcus granulosus. Chin Med J (Engl). 2011, 124 (18): 2838-2844.Google Scholar
- Zhang W, Li J, Jones MK, Zhang Z, Zhao L, Blair D, McManus DP: The Echinococcus granulosus antigen B gene family comprises at least 10 unique genes in five subclasses which are differentially expressed. PLoS Negl Trop Dis. 2010, 4 (8): e784-PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang W, Zhang Z, Shi B, Li J, You H, Tulson G, Dang X, Song Y, Yimiti T, Wang J, Jones MK, McManus DP: Vaccination of dogs against Echinococcus granulosus, the cause of cystic hydatid disease in humans. J Infect Dis. 2006, 194 (7): 966-974.PubMedView ArticleGoogle Scholar
- Wang WC, Lin FM, Chang WC, Lin KY, Huang HD, Lin NS: miRExpress: analyzing high-throughput sequencing data for profiling microRNA expression. BMC Bioinformatics. 2009, 10: 328-PubMed CentralPubMedView ArticleGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-PubMed CentralPubMedView ArticleGoogle Scholar
- Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, Rajewsky N: Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol. 2008, 26 (4): 407-415.PubMedView ArticleGoogle Scholar
- Gruber AR, Lorenz R, Bernhart SH, Neubock R, Hofacker IL: The Vienna RNA websuite. Nucleic Acids Res. 2008, 36 (Web Server issue): W70-74.PubMed CentralPubMedView ArticleGoogle Scholar
- Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones S, Marshall M, Matzke M, Ruvkun G, Tuschl T: A uniform system for microRNA annotation. RNA. 2003, 9 (3): 277-279.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang L, Feng Z, Wang X, Zhang X: DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 26 (1): 136-138.PubMedView ArticleGoogle Scholar
- Grant GR, Liu J, Stoeckert CJ: A practical false discovery rate approach to identifying patterns of differential expression in microarray data. Bioinformatics. 2005, 21 (11): 2684-2690.PubMedView ArticleGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676.PubMedView ArticleGoogle Scholar
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-PubMed CentralPubMedView ArticleGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35 (Web Server issue): W182-185.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.