Characterization of microRNAs in Taenia saginata of zoonotic significance by Solexa deep sequencing and bioinformatics analysis

microRNA silencing in a whole worm cestode model provides insight into miR-71 function

Parasites belonging to the class Cestoda include zoonotic species such as Echinococcus spp. and Taenia spp. that cause morbidity and mortality in endemic areas, mainly affecting pastoral and rural communities in low income countries but also upper middle income countries. Cestodes show remarkable developmental plasticity, implying tight regulation of gene expression throughout their complex life cycles. Despite the recent availability of genomic data for cestodes, little progress was made on postgenomic functional studies. MicroRNAs (miRNAs) are key components of gene regulatory systems that guide diverse developmental processes in multicellular organisms. miR-71 is a highly expressed miRNA in cestodes, which is absent in vertebrates and targets essential parasite genes, representing a potential key player in understanding the role of miRNAs in cestodes biology. Here we used transfection with antisense oligonucleotides to perform whole worm miRNA knockdown in tetrathyridia of Mesocestoides vogae (syn. Mesocestoides corti), a laboratory model of cestodes. We believe this is the first report of miRNA knockdown at the organism level in these parasites. Our results showed that M. vogae miR-71 is involved in the control of strobilation in vitro and in the establishment of murine infection. In addition, we identified miR-71 targets in M. vogae, several of them being de-repressed upon miR-71 knockdown. This study provides new knowledge on gene expression regulation in cestodes and suggests that miRNAs could be evaluated as new selective therapeutic targets for treating Neglected Tropical Diseases prioritised by the World Health Organization.

MicroRNAs in helminth parasites: a systematic review

BACKGROUND

MicroRNAs (miRNAs) are about 22-nucleotide, small, non-coding RNAs that control gene expression post-transcriptionally. Helminth parasites usually express a unique repertoire of genes, including miRNAs, across different developmental stages with subtle regulatory mechanisms.

OBJECTIVE

There is a necessity to investigate the involvement of miRNAs in the development of parasites, host-parasite interaction, immune evasion and their abilities to govern infection in hosts. miRNAs present in helminth parasites have been summarized in the current systematic review (SR).

METHODS

Electronic databases, including PubMed, Scopus, ProQuest, Embase, and Google Scholar search engine, were searched to identify helminth miRNA studies published from February 1993 till December 2019. Only the published articles in English were included in the study.

RESULTS

A total of 1769 articles were preliminarily recorded. Following the strict inclusion and exclusion criteria, 105 studies were included in this SR. Most of these studies focused on the identification of miRNAs in helminth parasites and/or probing of differentially expressed host miRNA profiles in specific relevant tissues, while 12 studies aimed to detect parasite-derived miRNAs in host circulating system and 15 studies characterized extracellular vesicles (EV)-derived miRNAs secreted by parasites.

CONCLUSION

In the current SR, information regarding all miRNAs expressed in helminth parasites has been comprehensively provided and the utility of helminth parasites-derived miRNAs in diagnosis and control of parasitic infections has been discussed. Furthermore, functional studies on helminth-derived miRNAs have also been presented.

Cestodes in the genomic era

The first cestode genomes were obtained by an international consortium led by the Wellcome Sanger Institute that included representative institutions from countries where the sequenced parasites have been studied for decades, in part because they are etiological agents of endemic diseases (Argentina, Uruguay, Mexico, Canada, UK, Germany, Switzerland, Ireland, USA, Japan, and China). After this, several complete genomes were obtained reaching 16 species to date. Cestode genomes have smaller relative size compared to other animals including free-living flatworms. Moreover, the features genome size and repeat content seem to differ in the two analyzed orders. Cyclophyllidean species have smaller genomes and with fewer repetitive content than Diphyllobothriidean species. On average, cestode genomes have 13,753 genes with 6 exons per gene and 41% GC content. More than 5,000 shared cestode proteins were accurately annotated by the integration of gene predictions and transcriptome evidence being more than 40% of these proteins of unknown function. Several gene losses and reduction of gene families were found and could be related to the extreme parasitic lifestyle of these species. The application of cutting-edge sequencing technology allowed the characterization of the terminal sequences of chromosomes that possess unique characteristics. Here, we review the current status of knowledge of complete cestode genomes and place it within a comparative genomics perspective. Multidisciplinary work together with the implementation of new technologies will provide valuable information that can certainly improve our chances to finally eradicate or at least control diseases caused by cestodes.

Taenia solium and Taenia crassiceps: miRNomes of the larvae and effects of miR-10-5p and let-7-5p on murine peritoneal macrophages

Neurocysticercosis (NCC), a major cause of neurological morbidity worldwide, is caused by the larvae of Taenia solium. Cestodes secrete molecules that block the Th1 response of their hosts and induce a Th2 response permissive to their establishment. Mature microRNAs (miRs) are small noncoding RNAs that regulate gene expression and participate in immunological processes. To determine the participation of Taenia miRs in the immune response against cysticercosis, we constructed small RNA (sRNA) libraries from larvae of Taenia solium and Taenia crassiceps. A total of 12074504 and 11779456 sequencing reads for T. solium and T. crassiceps, respectively, were mapped to the genomes of T. solium and other helminths. Both larvae shared similar miRNome, and miR-10-5p was the most abundant in both species, followed by let-7-5p in T. solium and miR-4989-3p in T. crassiceps, whereas among the genus-specific miRs, miR-001-3p was the most abundant in both, followed by miR-002-3p in T. solium and miR-003a-3p in T. crassiceps. The sequences of these miRs were identical in both. Structure and target prediction analyses revealed that these pre-miRs formed a hairpin and had more than one target involved in immunoregulation. Culture of macrophages, RT-PCR and ELISA assays showed that cells internalized miR-10-5p and let-7-5p into the cytoplasm and the miRs strongly decreased interleukin 16 (Il6) expression, tumor necrosis factor (TNF) and IL-12 secretion, and moderately decreased nitric oxide synthase inducible (Nos2) and Il1b expression (pro-inflammatory cytokines) in M(IFN-γ) macrophages and expression of Tgf1b, and the secretion of IL-10 (anti-inflammatory cytokines) in M(IL-4) macrophages. These findings could help us understand the role of miRs in the host–Taenia relationship.

Comparative Identification of MicroRNAs in Apis cerana cerana Workers' Midguts in Responseto Nosema ceranae Invasion

Here, the expression profiles and differentially expressed miRNAs (DEmiRNAs) in the midguts of Apis cerana cerana workers at 7 d and 10 d post-inoculation (dpi) with N. ceranae were investigated via small RNA sequencing and bioinformatics. Five hundred and twenty nine (529) known miRNAs and 25 novel miRNAs were identified in this study, and the expression of 16 predicted miRNAs was confirmed by Stem-loop RT-PCR. A total of 14 DEmiRNAs were detected in the midgut at 7 dpi, including eight up-regulated and six down-regulated miRNAs, while 12 DEmiRNAs were observed in the midgut at 10 dpi, including nine up-regulated and three down-regulated ones. Additionally, five DEmiRNAs were shared, while nine and seven DEmiRNAs were specifically expressed in midguts at 7 dpi and 10 dpi. Gene ontology analysis suggested some DEmiRNAs and corresponding target mRNAs were involved in various functions including immune system processes and response to stimulus. KEGG pathway analysis shed light on the potential functions of some DEmiRNAs in regulating target mRNAs engaged in material and energy metabolisms, cellular immunity and the humoral immune system. Further investigation demonstrated a complex regulation network between DEmiRNAs and their target mRNAs, with miR-598-y, miR-252-y, miR-92-x and miR-3654-y at the center. Our results can facilitate future exploration of the regulatory roles of miRNAs in host responses to N. ceranae, and provide potential candidates for further investigation of the molecular mechanisms underlying eastern honeybee-microsporidian interactions.

High-throughput identification of microRNAs in Taenia hydatigena, a cestode threatening livestock breeding industry

Infection of Cysticercus tenuicollis, the larval stage of Taenia hydatigena, is extensively found in sheep and pigs and jeopardizes the breeding and meat industry. miRNAs are a subclass of small noncoding regulatory RNAs and closely associated with the pathogenesis and biology of parasites. Here, using HiSeq sequencing we identified 49 known and 2 potential novel miRNAs in C. tenuicollis, of which both thy-miR-71 and -87 were predominant. Using RT-qPCR, 6 selected miRNAs were validated, and thy-miR-71 and -miR-87 were confirmed to be highly expressed, with the copy number of approximately 82,340 ± 2079 and 19,580 ± 609 per 1 ng total RNA, respectively. Similar to other cestodes, T. hydatigena was predicted to have two conserved miRNA clusters thy-miR-71/2c/2b and thy-miR-4989/277, and three members of the former were confirmed to reside sequentially within the genomic region of 253 bp by PCR. The current data provide us a valuable resource for further studies of a role of miRNAs in T. hydatigena biology and infection.

Conservation and diversification of small RNA pathways within flatworms

Background

Small non-coding RNAs, including miRNAs, and gene silencing mediated by RNA interference have been described in free-living and parasitic lineages of flatworms, but only few key factors of the small RNA pathways have been exhaustively investigated in a limited number of species. The availability of flatworm draft genomes and predicted proteomes allowed us to perform an extended survey of the genes involved in small non-coding RNA pathways in this phylum.

Results

Overall, findings show that the small non-coding RNA pathways are conserved in all the analyzed flatworm linages; however notable peculiarities were identified. While Piwi genes are amplified in free-living worms they are completely absent in all parasitic species. Remarkably all flatworms share a specific Argonaute family (FL-Ago) that has been independently amplified in different lineages. Other key factors such as Dicer are also duplicated, with Dicer-2 showing structural differences between trematodes, cestodes and free-living flatworms. Similarly, a very divergent GW182 Argonaute interacting protein was identified in all flatworm linages. Contrasting to this, genes involved in the amplification of the RNAi interfering signal were detected only in the ancestral free living species Macrostomum lignano. We here described all the putative small RNA pathways present in both free living and parasitic flatworm lineages.

Conclusion

These findings highlight innovations specifically evolved in platyhelminths presumably associated with novel mechanisms of gene expression regulation mediated by small RNA pathways that differ to what has been classically described in model organisms. Understanding these phylum-specific innovations and the differences between free living and parasitic species might provide clues to adaptations to parasitism, and would be relevant for gene-silencing technology development for parasitic flatworms that infect hundreds of million people worldwide.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1061-5) contains supplementary material, which is available to authorized users.

Identification and profiling of Cyprinus carpio microRNAs during ovary differentiation by deep sequencing

Background

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that regulate gene expression by targeting specific mRNAs. However, the possible role of miRNAs in the ovary differentiation and development of fish is not well understood. In this study, we examined the expression profiles and differential expression of miRNAs during three key stages of ovarian development and different developmental stages in common carp Cyprinus carpio.

Results

A total of 8765 miRNAs were identified, including 2155 conserved miRNAs highly conserved among various species, 145 miRNAs registered in miRBase for common carp, and 6505 novel miRNAs identified in common carp for the first time. Comparison of miRNA expression profiles among the five libraries identified 714 co-expressed and 2382 specific expressed miRNAs. Overall, 150, 628, and 431 specifically expressed miRNAs were identified in primordial gonad, juvenile ovary, and adult ovary, respectively. MiR-6758-3p, miR-3050-5p, and miR-2985-3p were highly expressed in primordial gonad, miR-3544-5p, miR-6877-3p, and miR-9086-5p were highly expressed in juvenile ovary, and miR-154-3p, miR-5307-5p, and miR-3958-3p were highly expressed in adult ovary. Predicted target genes of specific miRNAs in primordial gonad were involved in many reproductive biology signaling pathways, including transforming growth factor-β, Wnt, oocyte meiosis, mitogen-activated protein kinase, Notch, p53, and gonadotropin-releasing hormone pathways. Target-gene prediction revealed upward trends in miRNAs targeting male-bias genes, including dmrt1, atm, gsdf, and sox9, and downward trends in miRNAs targeting female-bias genes including foxl2, smad3, and smad4. Other sex-related genes such as sf1 were also predicted to be miRNA target genes.

Conclusions

This comprehensive miRNA transcriptome analysis demonstrated differential expression profiles of miRNAs during ovary development in common carp. These results could facilitate future exploitation of the sex-regulatory roles and mechanisms of miRNAs, especially in primordial gonads, while the specifically expressed miRNAs represent candidates for studying the mechanisms of ovary determination in Yellow River carp.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3701-y) contains supplementary material, which is available to authorized users.

High-throughput identification of miRNAs of Taenia ovis, a cestode threatening sheep industry

Taenia ovis is a tapeworm that is mainly transmitted between dogs and sheep or goats and has an adverse effect on sheep industry. miRNAs are short regulatory non-coding RNAs, involved in parasite development and growth as well as parasite infection. The miRNA profile of T. ovis remains to be established. Herein, 33 known miRNAs belonging to 23 different families were identified in T. ovis metacestodes using deep sequencing approach. Of them, expression of some miRNAs such as tov-miR-10 and -let-7 was absolutely predominant. Moreover, comparative analysis revealed the presence of a miR-71/2b/2c cluster in T. ovis, which was also completely conserved in other 6 cestodes. The study provides rich data for further understandings of T. ovis biology.

High-throughput sequencing of RNAs isolated by cross-linking immunoprecipitation (HITS-CLIP) reveals Argonaute-associated microRNAs and targets in Schistosoma japonicum

Background

Schistosomiasis, caused by schistosomes, is one of the most prevalent and serious parasitic diseases in tropical and subtropical countries. This pathogen has a complex life cycle and harbors a unique repertoire of genes expressed at different life-stages. Understanding the gene regulation of schistosomes will contribute to identification of novel drug targets and vaccine candidates. Some conserved and novel microRNAs (miRNAs) have been identified in schistosomes as key transcriptional and post-transcriptional regulators in the past few years; however, little is known about their specific targets.

Methods

High-throughput sequencing of RNAs isolated by cross-linking immunoprecipitation (HITS-CLIP) was used to covalently crosslink native Argonaute protein-RNA complexes in Schistosoma japonicum. An antibody against S.japonicum Argonaute proteins, was generated and used for immunoprecipitation of the crosslinked SjAgo-RNA complex from soluble adult worm extract. Small RNAs, including miRNAs and their target mRNAs associated with the native SjAgo in adult parasites, were enriched and extracted for library construction.

Results

High-throughput sequencing produced a total of ~7.4 million high-quality reads, of which approximately 45.07 % were composed of 769 miRNAs and 35.54 % were composed of 11,854 mRNAs target sites. Further bioinformatics analysis identified 43 conserved known miRNAs and 256 novel miRNAs in the SjAgo-associated small RNA population. An average of approximately 15 target sites were predicted for each miRNA. Moreover, a positive rate of 50 % has been achieved in a small-scale verification test of the putative target sites of miRNA1.

Conclusion

In this study, we isolated and identified small RNAs including miRNAs and their targets associated with the S. japonicum Argonaute proteins, by the HITS-CLIP method combined with bioinformatics and biologic experimental analysis. These data reveal a genome-wide miRNA-mRNA interaction map in S. japonicum in vivo, which will help us understand the complex gene regulatory network in this pathogen and thereby facilitate the development of novel drug approaches against schistosomiasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1203-9) contains supplementary material, which is available to authorized users.

The miRnome of Fasciola hepatica juveniles endorses the existence of a reduced set of highly divergent micro RNAs in parasitic flatworms

The liver fluke Fasciola hepatica is a foodborne zoonotic parasite affecting livestock worldwide, with increasing relevance in human health. The first developmental stage that the host meets after ingestion of the parasite is the newly excysted juvenile, that actively transverses the gut wall and migrates to its final location in the liver. The regulation of the early developmental events in newly excysted juveniles is still poorly understood and a relevant target for control strategies. Here we investigated the putative involvement of small regulatory RNAs in the invasion process. The small RNA population of the newly excysted juvenile fall into two classes, one represented by micro (mi)RNAs and a secondary group of larger (32-33 nucleotides) tRNA-derived sequences. We identified 40 different miRNAs, most of those belonging to ancient miRNAs conserved in protostomes and metazoans, notably with a highly predominant miR-125b variant. Remarkably, several protostomian and metazoan conserved families were not detected in consonance with previous reports of drastic miRnome reduction in parasitic flatworms. Additionally, a set of five novel miRNAs was identified, probably associated with specific gene regulation expression needs in F. hepatica. While sequence conservation in mature miRNA is high across the metazoan tree, we observed that flatworm miRNAs are more divergent, suggesting that mutation rates in parasitic flatworms could be high. Finally, the distinctive presence of tRNA-derived sequences, mostly 5' tRNA halves of selected tRNAs in the small RNA population of newly excysted juveniles, raises the possibility that both miRNA and tRNA fragments participate in the regulation of gene expression in this parasite.

MicroRNA-Sequence Profiling Reveals Novel Osmoregulatory MicroRNA Expression Patterns in Catadromous Eel Anguilla marmorata

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate gene expression by post-transcriptional repression of mRNAs. Recently, several miRNAs have been confirmed to execute directly or indirectly osmoregulatory functions in fish via translational control. In order to clarify whether miRNAs play relevant roles in the osmoregulation of Anguilla marmorata, three sRNA libraries of A. marmorata during adjusting to three various salinities were sequenced by Illumina sRNA deep sequencing methods. Totally 11,339,168, 11,958,406 and 12,568,964 clear reads were obtained from 3 different libraries, respectively. Meanwhile, 34 conserved miRNAs and 613 novel miRNAs were identified using the sequence data. MiR-10b-5p, miR-181a, miR-26a-5p, miR-30d and miR-99a-5p were dominantly expressed in eels at three salinities. Totally 29 mature miRNAs were significantly up-regulated, while 72 mature miRNAs were significantly down-regulated in brackish water (10‰ salinity) compared with fresh water (0‰ salinity); 24 mature miRNAs were significantly up-regulated, while 54 mature miRNAs were significantly down-regulated in sea water (25‰ salinity) compared with fresh water. Similarly, 24 mature miRNAs were significantly up-regulated, while 45 mature miRNAs were significantly down-regulated in sea water compared with brackish water. The expression patterns of 12 dominantly expressed miRNAs were analyzed at different time points when the eels transferred from fresh water to brackish water or to sea water. These miRNAs showed differential expression patterns in eels at distinct salinities. Interestingly, miR-122, miR-140-3p and miR-10b-5p demonstrated osmoregulatory effects in certain salinities. In addition, the identification and characterization of differentially expressed miRNAs at different salinities can clarify the osmoregulatory roles of miRNAs, which will shed lights for future studies on osmoregulation in fish.

Application of small RNA technology for improved control of parasitic helminths

•

MicroRNAs and siRNAs in helminth post-transcriptional gene regulation are reviewed.

•

Many parasitic helminth miRNAs are unique and developmentally expressed.

•

miRNAs released by parasites have diagnostic potential, particularly for filarial and schistosome spp.

•

Parasite and host miRNAs may regulate immune responses.

•

Improvements to siRNA-mediated gene silencing are important for functional genomics.

Over the last decade microRNAs (miRNAs) and small interfering RNAs (siRNAs) have emerged as important regulators of post-transcriptional gene expression. miRNAs are short, non-coding RNAs that regulate a variety of processes including cancer, organ development and immune function. This class of small RNAs bind with partial complementarity to their target mRNA sequences, most often in the 3′UTR, to negatively regulate gene expression. In parasitic helminths, miRNAs are being increasingly studied for their potential roles in development and host-parasite interactions. The availability of genome data, combined with small RNA sequencing, has paved the way to profile miRNAs expressed at particular developmental stages for many parasitic helminths. While some miRNAs are conserved across species, others appear to be unique to specific parasites, suggesting important roles in adaptation and survival in the host environment. Some miRNAs are released from parasites, in exosomes or in protein complexes, and the potential effects of these on host immune function are being increasingly studied. In addition, release of miRNAs from schistosome and filarial parasites into host plasma can be exploited for the development of specific and sensitive diagnostic biomarkers of infection. Interfering with miRNA function, as well as silencing key components of the pathways they regulate, will progress our understanding of parasite development and provide a novel approach to therapeutic control. RNA interference (RNAi) by siRNAs has proven to be inconsistent in parasitic nematodes. However, the recent successes reported for schistosome and liver fluke RNAi, encourage further efforts to enhance delivery of RNA and improve in vitro culture systems and assays to monitor phenotypic effects in nematodes. These improvements are important for the establishment of reliable functional genomic platforms for novel drug and vaccine development. In this review we focus on the important roles of miRNAs and siRNAs in post-transcriptional gene regulation in veterinary parasitic helminths and the potential value of these in parasite diagnosis and control.

Identification of microRNA genes in three opisthorchiids

Background

Opisthorchis felineus, O. viverrini, and Clonorchis sinensis (family Opisthorchiidae) are parasitic flatworms that pose a serious threat to humans in some countries and cause opisthorchiasis/clonorchiasis. Chronic disease may lead to a risk of carcinogenesis in the biliary ducts. MicroRNAs (miRNAs) are small noncoding RNAs that control gene expression at post-transcriptional level and are implicated in the regulation of various cellular processes during the parasite- host interplay. However, to date, the miRNAs of opisthorchiid flukes, in particular those essential for maintaining their complex biology and parasitic mode of existence, have not been satisfactorily described.

Methodology/Principal Findings

Using a SOLiD deep sequencing-bioinformatic approach, we identified 43 novel and 18 conserved miRNAs for O. felineus (miracidia, metacercariae and adult worms), 20 novel and 16 conserved miRNAs for O. viverrini (adult worms), and 33 novel and 18 conserved miRNAs for C. sinensis (adult worms). The analysis of the data revealed differences in the expression level of conserved miRNAs among the three species and among three the developmental stages of O. felineus. Analysis of miRNA genes revealed two gene clusters, one cluster-like region and one intronic miRNA in the genome. The presence and structure of the two gene clusters were validated using a PCR-based approach in the three flukes.

Conclusions

This study represents a comprehensive description of miRNAs in three members of the family Opistorchiidae, significantly expands our knowledge of miRNAs in multicellular parasites and provides a basis for understanding the structural and functional evolution of miRNAs in these metazoan parasites. Results of this study also provides novel resources for deeper understanding the complex parasite biology, for further research on the pathogenesis and molecular events of disease induced by the liver flukes. The present data may also facilitate the development of novel approaches for the prevention and treatment of opisthorchiasis/clonorchiasis.

Liver flukes of the family Opisthorchiidae cause diseases of the hepatobiliary system, known as opisthorchiasis/clonorchiasis. The chronic forms of these diseases greatly increase the risk of cancer developing in the biliary ducts. Much has been elucidated regarding the developmental biology of opisthorchiid flukes and the molecular pathological effects on the definitive host; however, the role of microRNAs (short non-coding RNAs) capable of influencing the pathogenic process and host-parasite interactions have not yet been comprehensively studied. The aim of the present work was to identify the miRNA genes of the liver flukes and provide a basis for further investigating the roles of these miRNAs in the complex opisthorchiidae life cycle and the pathogenesis of disease.

microRNA profiling in the zoonotic parasite Echinococcus canadensis using a high-throughput approach

Background

microRNAs (miRNAs), a class of small non-coding RNAs, are key regulators of gene expression at post-transcriptional level and play essential roles in fundamental biological processes such as development and metabolism. The particular developmental and metabolic characteristics of cestode parasites highlight the importance of studying miRNA gene regulation in these organisms. Here, we perform a comprehensive analysis of miRNAs in the parasitic cestode Echinococcus canadensis G7, one of the causative agents of the neglected zoonotic disease cystic echinococcosis.

Methods

Small RNA libraries from protoscoleces and cyst walls of E. canadensis G7 and protoscoleces of E. granulosus sensu stricto G1 were sequenced using Illumina technology. For miRNA prediction, miRDeep2 core algorithm was used. The output list of candidate precursors was manually curated to generate a high confidence set of miRNAs. Differential expression analysis of miRNAs between stages or species was estimated with DESeq. Expression levels of selected miRNAs were validated using poly-A RT-qPCR.

Results

In this study we used a high-throughput approach and found transcriptional evidence of 37 miRNAs thus expanding the miRNA repertoire of E. canadensis G7. Differential expression analysis showed highly regulated miRNAs between life cycle stages, suggesting a role in maintaining the features of each developmental stage or in the regulation of developmental timing. In this work we characterize conserved and novel Echinococcus miRNAs which represent 30 unique miRNA families. Here we confirmed the remarkable loss of conserved miRNA families in E. canadensis, reflecting their low morphological complexity and high adaptation to parasitism.

Conclusions

We performed the first in-depth study profiling of small RNAs in the zoonotic parasite E. canadensis G7. We found that miRNAs are the preponderant small RNA silencing molecules, suggesting that these small RNAs could be an essential mechanism of gene regulation in this species. We also identified both parasite specific and divergent miRNAs which are potential biomarkers of infection. This study will provide valuable information for better understanding of the complex biology of this parasite and could help to find new potential targets for therapy and/or diagnosis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0686-8) contains supplementary material, which is available to authorized users.

Detection of circulating parasite-derived microRNAs in filarial infections

Filarial nematodes cause chronic and profoundly debilitating diseases in both humans and animals. Applications of novel technology are providing unprecedented opportunities to improve diagnosis and our understanding of the molecular basis for host-parasite interactions. As a first step, we investigated the presence of circulating miRNAs released by filarial nematodes into the host bloodstream. miRNA deep-sequencing combined with bioinformatics revealed over 200 mature miRNA sequences of potential nematode origin in Dirofilaria immitis-infected dog plasma in two independent analyses, and 21 in Onchocerca volvulus-infected human serum. Total RNA obtained from D. immitis-infected dog plasma was subjected to stem-loop RT-qPCR assays targeting two detected miRNA candidates, miR-71 and miR-34. Additionally, Brugia pahangi-infected dog samples were included in the analysis, as these miRNAs were previously detected in extracts prepared from this species. The presence of miR-71 and miR-34 discriminated infected samples (both species) from uninfected samples, in which no specific miRNA amplification occurred. However, absolute miRNA copy numbers were not significantly correlated with microfilaraemia for either parasite. This may be due to the imprecision of mf counts to estimate infection intensity or to miRNA contributions from the unknown number of adult worms present. Nonetheless, parasite-derived circulating miRNAs are found in plasma or serum even for those species that do not live in the bloodstream.

Filarial parasites commonly infect humans and animals, especially in tropical settings. The strongly debilitating panel of diseases they cause in humans contributes to an entrenched cycle of poverty. For efficient treatment strategies, reliable diagnostic tests are necessary. We investigated the potential of parasite-derived microRNAs (miRNAs; short non-coding RNA molecules present in eukaryotes) as biomarkers of infection. Using deep-sequencing technologies and bioinformatics, we identified over two-hundred mature miRNA candidates of nematode origin in plasma from Dirofilaria immitis-infected dogs. Similarly, we discovered twenty-one miRNA candidates predicted to be released by Onchocerca volvulus in infected human sera. We developed two RT-qPCR assays for the detection of D. immitis miR-71 and miR-34 in dog plasma that discriminated infected from uninfected samples. We demonstrated the presence of filarial miRNAs in host blood, regardless of localization in their respective hosts, and suggest that they are suitable targets for detection by RT-qPCR.

Identification and characterization of microRNAs in normal equine tissues by Next Generation Sequencing

The role of microRNAs (miRNAs) as a post-transcriptional gene regulator has been elucidated in a broad range of organisms including domestic animals. Characterization of miRNAs in normal tissues is an important step to investigate the functions of miRNAs in various physiological and pathological conditions. Using Illumina Next Generation Sequencing (NGS) technology, we identified a total of 292 known and 329 novel miRNAs in normal horse tissues including skeletal muscle, colon and liver. Distinct sets of miRNAs were differentially expressed in a tissue-specific manner. The miRNA genes were distributed across all the chromosomes except chromosomes 29 and 31 in the horse reference genome. In some chromosomes, multiple miRNAs were clustered and considered to be polycistronic transcript. A base composition analysis showed that equine miRNAs had a higher frequency of A+U than G+C. Furthermore, U tended to be more frequent at the 5′ end of miRNA sequences. This is the first experimental study that identifies and characterizes the global miRNA expression profile in normal horse tissues. The present study enriches the horse miRNA database and provides useful information for further research dissecting biological functions of miRNAs in horse.

Identification and characterization of the expression profile of microRNAs in Anopheles anthropophagus

Background

Anopheles anthropophagus, one of the most important mosquito-borne disease vectors in Asia, mainly takes blood meals from humans and transmits both malaria and filariae. MicroRNAs (miRNAs) are small non-coding RNAs, and play a critical role in many cellular processes, including development, differentiation, apoptosis and innate immunity.

Methods

We investigated the global miRNA expression profile of male and female adults of A. anthropophagus using illumina Hiseq2000 sequencing combined with Northern blot.

Results

By using the miRNAs of the closely-related species Anopheles gambiae and Aedes aegypti as reference, we obtained 102 miRNAs candidates out of 12.43 million raw sequencing reads for male and 16.51 million reads for female, with 81 of them found as known miRNAs in An. gambiae and/or Ae. aegypti, and the remaining 21 miRNAs were considered as novel. By analyzing the revised read count of miRNAs in male and female, 29 known miRNAs show sexual difference expression: >2-fold in the read count of the same miRNAs in male and female. Especially for miR-989, which is highly expressed in the female mosquitoes, but shows almost no detected expression in male mosquitoes, indicating that miR-989 may be involved in the physiological activity of female mosquito adults. The expression of four miRNAs in different growth stages of mosquito were further identified by Northern blot. Several miRNAs show the stage-specific expression, of which miR-2943 only expressed in the egg stage, suggesting that miR-2943 may be associated with the development of mosquito eggs.

Conclusions

The present study represents the first global characterization of An. anthropophagus miRNAs in sexual differences and stage-specific functions. A better understanding of the functions of these miRNAs will offer new insights in mosquito biology and has implications for the effective control of mosquito-borne infectious diseases.

microRNAs of parasitic helminths - Identification, characterization and potential as drug targets

•

Importance of microRNAs in helminth post-transcriptional gene regulation is reviewed.

•

Increasing helminth miRNA data are available from deep sequencing.

•

Some miRNAs are helminth-specific, many are novel to each species.

•

miRNAs may regulate parasite and host gene expression.

•

Uptake of miRNA inhibitors and mimics is feasible for functional analysis.

microRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional gene regulation. They were first identified in the free-living nematode Caenorhabditis elegans, where the miRNAs lin-4 and let-7 were shown to be essential for regulating correct developmental progression. The sequence of let-7 was subsequently found to be conserved in higher organisms and changes in expression of let-7, as well as other miRNAs, are associated with certain cancers, indicating important regulatory roles. Some miRNAs have been shown to have essential functions, but the roles of many are currently unknown. With the increasing availability of genome sequence data, miRNAs have now been identified from a number of parasitic helminths, by deep sequencing of small RNA libraries and bioinformatic approaches. While some miRNAs are widely conserved in a range of organisms, others are helminth-specific and many are novel to each species. Here we review the potential roles of miRNAs in regulating helminth development, in interacting with the host environment and in development of drug resistance. Use of fluorescently-labeled small RNAs demonstrates uptake by parasites, at least in vitro. Therefore delivery of miRNA inhibitors or mimics has potential to alter miRNA activity, providing a useful tool for probing the roles of miRNAs and suggesting novel routes to therapeutics for parasite control.

Identification and profiling of sex-biased microRNAs from sea urchin Strongylocentrotus nudus gonad by Solexa deep sequencing

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate gene expression by post-transcriptional repression of messenger RNA. The echinoderm, Strongylocentrotus nudus, is an excellent model organism for studying development and commercially important as a food source. However, to date no miRNAs have been reported to modulate sex gonad differentiation in S. nudus. In this study, we constructed two small RNA libraries from male and female S. nudus gonad respectively for Solexa sequencing. A total of 184 miRNAs including 60 known and 124 novel miRNAs were identified from the two libraries. Furthermore, the nucleotide bias and end variation of the known miRNAs were also analyzed. In addition, 67 differently expressed of the 86 co-expressed and 98 gender-specific (47 male-specific and 51 female-specific) miRNAs that may be involved in sexual differentiation were found by comparing the miRNA expression profiles in the two libraries. This study reveals the first miRNA profile related to the gonad differentiation of the S. nudus. This study gives a first insight into sex differences in miRNA expression of sea urchin which could facilitate studies of the reproductive organ-specific roles of miRNAs.

Substantial loss of conserved and gain of novel MicroRNA families in flatworms

Recent studies on microRNA (miRNA) evolution focused mainly on the comparison of miRNA complements between animal clades. However, evolution of miRNAs within such groups is poorly explored despite the availability of comparable data that in some cases lack only a few key taxa. For flatworms (Platyhelminthes), miRNA complements are available for some free-living flatworms and all major parasitic lineages, except for the Monogenea. We present the miRNA complement of the monogenean flatworm Gyrodactylus salaris that facilitates a comprehensive analysis of miRNA evolution in Platyhelminthes. Using the newly designed bioinformatics pipeline miRCandRef, the miRNA complement was disentangled from next-generation sequencing of small RNAs and genomic DNA without a priori genome assembly. It consists of 39 miRNA hairpin loci of conserved miRNA families, and 22 novel miRNAs. A comparison with the miRNA complements of Schmidtea mediterranea (Turbellaria), Schistosoma japonicum (Trematoda), and Echinococcus granulosus (Cestoda) reveals a substantial loss of conserved bilaterian, protostomian, and lophotrochozoan miRNAs. Eight of the 46 expected conserved miRNAs were lost in all flatworms, 16 in Neodermata and 24 conserved miRNAs could not be detected in the cestode and the trematode. Such a gradual loss of miRNAs has not been reported before for other animal phyla. Currently, little is known about miRNAs in Platyhelminthes, and for the majority of the lost miRNAs there is no prediction of function. As suggested earlier they might be related to morphological simplifications. The presence and absence of 153 conserved miRNAs was compared for platyhelminths and 32 other metazoan taxa. Phylogenetic analyses support the monophyly of Platyhelminthes (Turbellaria + Neodermata [Monogenea {Trematoda + Cestoda}]).

Identification of neglected cestode Taenia multiceps microRNAs by illumina sequencing and bioinformatic analysis

Background

Worldwide, but especially in developing countries, coenurosis of sheep and other livestock is caused by Taenia multiceps larvae, and zoonotic infections occur in humans. Infections frequently lead to host death, resulting in huge socioeconomic losses. MicroRNAs (miRNAs) have important roles in the post-transcriptional regulation of a large number of animal genes by imperfectly binding target mRNAs. To date, there have been no reports of miRNAs in T. multiceps.

Results

In this study, we obtained 12.8 million high quality raw reads from adult T. multiceps small RNA library using Illumina sequencing technology. A total of 796 conserved miRNA families (containing 1,006 miRNAs) from 170,888 unique miRNAs were characterized using miRBase (Release 17.0). Here, we selected three conserved miRNA/miRNA* (antisense strand) duplexes at random and amplified their corresponding precursors using a PCR-based method. Furthermore, 20 candidate novel miRNA precursors were verified by genomic PCR. Among these, six corresponding T. multiceps miRNAs are considered specific for Taeniidae because no homologs were found in other species annotated in miRBase. In addition, 181,077 target sites within T. multiceps transcriptome were predicted for 20 candidate newly miRNAs.

Conclusions

Our large-scale investigation of miRNAs in adult T. multiceps provides a substantial platform for improving our understanding of the molecular regulation of T. multiceps and other cestodes development.

Identification of Dirofilaria immitis miRNA using illumina deep sequencing

The heartworm Dirofilaria immitis is the causative agent of cardiopulmonary dirofilariosis in dogs and cats, which also infects a wide range of wild mammals and humans. The complex life cycle of D. immitis with several developmental stages in its invertebrate mosquito vectors and its vertebrate hosts indicates the importance of miRNA in growth and development, and their ability to regulate infection of mammalian hosts. This study identified the miRNA profiles of D. immitis of zoonotic significance by deep sequencing. A total of 1063 conserved miRNA candidates, including 68 anti-sense miRNA (miRNA*) sequences, were predicted by computational methods and could be grouped into 808 miRNA families. A significant bias towards family members, family abundance and sequence nucleotides was observed. Thirteen novel miRNA candidates were predicted by alignment with the Brugia malayi genome. Eleven out of 13 predicted miRNA candidates were verified by using a PCR-based method. Target genes of the novel miRNA candidates were predicted by using the heartworm transcriptome dataset. To our knowledge, this is the first report of miRNA profiles in D. immitis, which will contribute to a better understanding of the complex biology of this zoonotic filarial nematode and the molecular regulation roles of miRNA involved. Our findings may also become a useful resource for small RNA studies in other filarial parasitic nematodes.

Comparative profiling of microRNAs in male and female adults of Ascaris suum

Ascaris nematodes, which cause ascariasis in humans and pigs, are among the most important nematodes from both health and economic perspectives. microRNA (miRNA) is now recognized as key regulator of gene expression at posttranscription level. The public availability of the genome and transcripts of Ascaris suum provides powerful resources for the research of miRNA profiles of the parasite. Therefore, we investigated and compared the miRNA profiles of male and female adult A. suum using Solexa deep sequencing combined with bioinformatic analysis and stem-loop reverse transcription polymerase chain reaction. Deep sequencing of small RNAs yielded 11.71 and 11.72 million raw reads from male and female adults of A. suum, respectively. Analysis showed that the noncoding RNA of the two genders, including tRNA, rRNA, snRNA, and snoRNA, were similar. By mapping to the A. suum genome, we obtained 494 and 505 miRNA candidates from the female and male parasite, respectively, and 87 and 82 of miRNA candidates were consistent with A. suum miRNAs deposited in the miRBase database. Among the miRNA candidates, 154 were shared by the two genders, and 340 and 351 were female and male specific with their target numbers ranged from one to thousands, respectively. Functional prediction revealed a set of elongation factors, heat shock proteins, and growth factors from the targets of gender-specific miRNAs, which were essential for the development of the parasite. Moreover, major sperm protein and nematode sperm cell motility protein were found in targets of the male-specific miRNAs. Ovarian message protein was found in targets of the female-specific miRNAs. Enrichment analysis revealed significant differences among Gene Ontology terms of miRNA targets of the two genders, such as electron carrier and biological adhesion process. The regulating functions of gender-specific miRNAs was therefore not only related to the fundamental functions of cells but also were essential to the germ development of the parasite. The present study provides a framework for further research of Ascaris miRNAs, and consequently leads to the development of potential nucleotide vaccines against Ascaris of human and animal health significance.

Phylogenetic analysis of the Argonaute protein family in platyhelminths

Argonaute proteins (AGOs) are mediators of gene silencing via recruitment of small regulatory RNAs to induce translational regression or degradation of targeted molecules. Platyhelminths have been reported to express microRNAs but the diversity of AGOs in the phylum has not been explored. Phylogenetic relationships of members of this protein family were studied using data from six platyhelminth genomes. Phylogenetic analysis showed that all cestode and trematode AGOs, along with some triclad planarian AGOs, were grouped into the Ago subfamily and its novel sister clade, here referred to as Cluster 1. These were very distant from Piwi and Class 3 subfamilies. By contrast, a number of planarian Piwi-like AGOs formed a novel sister clade to the Piwi subfamily. Extensive sequence searching revealed the presence of an additional locus for AGO2 in the cestode Echinococcus granulosus and exon expansion in this species and E. multilocularis. The current study suggests the absence of the Piwi subfamily and Class 3 AGOs in cestodes and trematodes and the Piwi-like AGO expansion in a free-living triclad planarian and the occurrence of exon expansion prior to or during the evolution of the most-recent common ancestor of the Echinococcus species studied.

Epigenetics: a key regulator of platyhelminth developmental biology?

The Platyhelminthes (flukes/flatworms) are a large group of derived metazoans beautifully adapted for existence in diversely challenging ecosystems. As tractable examples of development and self-regeneration or as causative agents of aquacultural, veterinary and biomedically-relevant parasitic diseases, the platyhelminths are subject to intensive inter-disciplinary research. Given the complex lifestyles exhibited by individuals within this phylum, we postulate that epigenetic processes feature in many aspects of platyhelminth lifecycle diversity, development and environmentally-driven adaptations.