Parasite-Derived MicroRNAs in Host Serum As Novel Biomarkers of Helminth Infection: e2701

Abstract

Background

MicroRNAs (miRNAs) are a class of short non-coding RNA that play important roles in disease processes in animals and are present in a highly stable cell-free form in body fluids. Here, we examine the capacity of host and parasite miRNAs to serve as tissue or serum biomarkers of Schistosoma mansoni infection.

Methods/Principal Findings

We used Exiqon miRNA microarrays to profile miRNA expression in the livers of mice infected with S. mansoni at 7 weeks post-infection. Thirty-three mouse miRNAs were differentially expressed in infected compared to naïve mice (>2 fold change, p<0.05) including miR-199a-3p, miR-199a-5p, miR-214 and miR-21, which have previously been associated with liver fibrosis in other settings. Five of the mouse miRNAs were also significantly elevated in serum by twelve weeks post-infection. Sequencing of small RNAs from serum confirmed the presence of these miRNAs and further revealed eleven parasite-derived miRNAs that were detectable by eight weeks post infection. Analysis of host and parasite miRNA abundance by qRT-PCR was extended to serum of patients from low and high infection sites in Zimbabwe and Uganda. The host-derived miRNAs failed to distinguish uninfected from infected individuals. However, analysis of three of the parasite-derived miRNAs (miR-277, miR-3479-3p and bantam) could detect infected individuals from low and high infection intensity sites with specificity/sensitivity values of 89%/80% and 80%/90%, respectively.

Conclusions

This work identifies parasite-derived miRNAs as novel markers of S. mansoni infection in both mice and humans, with the potential to be used with existing techniques to improve S. mansoni diagnosis. In contrast, although host miRNAs are differentially expressed in the liver during infection their abundance levels in serum are variable in human patients and may be useful in cases of extreme pathology but likely hold limited value for detecting prevalence of infection.

Abstract

Background

MicroRNAs (miRNAs) are a class of short non-coding RNA that play important roles in disease processes in animals and are present in a highly stable cell-free form in body fluids. Here, we examine the capacity of host and parasite miRNAs to serve as tissue or serum biomarkers of Schistosoma mansoni infection.

Methods/Principal Findings

We used Exiqon miRNA microarrays to profile miRNA expression in the livers of mice infected with S. mansoni at 7 weeks post-infection. Thirty-three mouse miRNAs were differentially expressed in infected compared to naïve mice (>2 fold change, p<0.05) including miR-199a-3p, miR-199a-5p, miR-214 and miR-21, which have previously been associated with liver fibrosis in other settings. Five of the mouse miRNAs were also significantly elevated in serum by twelve weeks post-infection. Sequencing of small RNAs from serum confirmed the presence of these miRNAs and further revealed eleven parasite-derived miRNAs that were detectable by eight weeks post infection. Analysis of host and parasite miRNA abundance by qRT-PCR was extended to serum of patients from low and high infection sites in Zimbabwe and Uganda. The host-derived miRNAs failed to distinguish uninfected from infected individuals. However, analysis of three of the parasite-derived miRNAs (miR-277, miR-3479-3p and bantam) could detect infected individuals from low and high infection intensity sites with specificity/sensitivity values of 89%/80% and 80%/90%, respectively.

Conclusions

This work identifies parasite-derived miRNAs as novel markers of S. mansoni infection in both mice and humans, with the potential to be used with existing techniques to improve S. mansoni diagnosis. In contrast, although host miRNAs are differentially expressed in the liver during infection their abundance levels in serum are variable in human patients and may be useful in cases of extreme pathology but likely hold limited value for detecting prevalence of infection.

Author Summary

Schistosomiasis is a chronic disease caused by blood flukes that affects over 200 million people worldwide, of which 90% live in Sub-Saharan Africa. In the field setting schistosomiasis caused by S. mansoni is diagnosed by detection of parasite eggs in stool samples using microscopic techniques. Here we investigate the potential of microRNAs (miRNAs), a class of short noncoding RNAs, to act as biomarkers of S. mansoni infection. We have identified a specific subset of murine miRNAs whose expression is significantly altered in the liver between 6-12 weeks post infection. However their abundance in serum is not significantly different between naïve and S. mansoni-infected mice until twelve weeks post infection and they do not display consistent differential abundance in the serum of infected versus uninfected humans. In contrast, three parasite-derived miRNAs (miR-277, bantam and miR-3479-3p) were detected in the serum of infected mice and human patients and the combined detection of these miRNAs could distinguish S. mansoni infected from uninfected individuals from low and high infection intensity areas with 89%/80% or 80%/90% specificity/sensitivity, respectively. These results demonstrate that miRNAs of parasite origin are a new class of serum biomarker for detecting S. mansoni and likely other helminth infections.

Citation: Hoy AM, Lundie RJ, Ivens A, Quintana JF, Nausch N, et al. (2014) Parasite-Derived MicroRNAs in Host Serum As Novel Biomarkers of Helminth Infection. PLoS Negl Trop Dis 8(2): e2701. doi:10.1371/journal.pntd.0002701

Editor: William Evan Secor, Centers for Disease Control and Prevention, Atlanta, United States of America

Received: July 4, 2013; Accepted: January 3, 2014; Published: February 20, 2014

Copyright: © 2014 Hoy et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work was supported by a NERC studentship (http://www.nerc.ac.uk/) to AMH, funding by the Thrasher Research Fund (https://www.thrasherresearch.org/) to FM and NN, a MRC SNCF to AsM (G0701437) and a Wellcome Trust RCDF to AHB (097394/Z/11/Z). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Helminths are parasitic worms that infect a third of the world's population and cause a diverse range of health consequences leading to significant social and economical burdens [1], [2]. Schistosomiasis is a chronic disease caused by blood flukes of the genus Schistosoma that affects more than 200 million people worldwide and is second only to malaria as the most important lethal human parasitic disease in tropical and subtropical regions. Schistosomiasis is predominantly caused by hepatic S. mansoni and urogenital S. haematobium [3]. It is estimated that the mortality rates due to haematemesis (S. mansoni) and renal failure (S. haematobium) are around 130,000 and 150,000 per year respectively [4]. In addition, schistosomiasis is associated with anaemia, diarrhoea, under nutrition, chronic pain and exercise intolerance, which are estimated to contribute to 0.02-0.15 disability-adjusted life-years (DALY) [5].

Diagnosis of schistosome infection is crucial for patient management, evaluation of treatment efficiency, monitoring of disease transmission and success of control strategies, as recommended by the World Health Organization [6]. In the field, S. mansoni and other intestinal schistosomes are currently diagnosed through the detection of the parasite eggs in stool specimens using microscopic techniques such as Kato-Katz or ether-concentration [7]. While these techniques are relatively simple, inexpensive and specific, their major drawbacks include poor sensitivity in detecting low-intensity infections (for example, in children), their inability to detect pre-patent or single sex infection and their failure to detect infection in individuals where eggs are trapped in tissues and not excreted [8]. Available antibody-based assays are useful for diagnosis in some cases (e.g. foreign travellers) but they cannot differentiate past and active infection and can also cross-react with antigens from other helminths [9]. These assays therefore do not offer a definitive diagnosis in schistosome-endemic areas. Recent studies have shown success with point-of-care tests for Schistosoma circulating cathodic and anodic antigens (CCA and CAA, respectively) in serum and urine, which decrease rapidly after chemotherapy [8], [10], [11], [12], [13], [14] and these are now being further developed for use in the field [15]. Detection of schistosome DNA in urine and stool samples or plasma by the polymerase chain reaction (PCR) method is another strategy for routine diagnosis of infection that has shown promising results [16], [17], [18], [19], [20]. Here we examine whether microRNAs (miRNAs), which are extremely stable in serum and detectable by quantitative reverse transcription PCR (qRT-PCR) could provide an additional diagnostic tool for S. mansoni infection.

miRNAs are a class of naturally occurring small non-coding RNA produced from animal, plant and viral genomes [21]. They are incorporated into the RNA-induced silencing complex (RISC) and function by binding to messenger RNAs (mRNAs) and inhibiting translation and/or causing mRNA destabilization [22]. Depending on the genes they target, miRNAs have diverse functions inside cells, from regulation of developmental programming to viral-host interactions [23], [24], [25]. In the last 5 years, reports have shown that miRNAs circulate in serum in a cell-free form and many cell types secrete miRNAs through encapsulation within exosomes or in association with specific proteins [26], [27]. These extracellular RNAs have been implicated in cell-to-cell communication in a range of systems [28], [29], and have been shown to be extremely stable in body fluids [27], [30]. They have received extensive attention for their promise as biomarkers of disease, including cancer [31], [32], [33]. Extracellular miRNAs have also been shown to be altered in human serum and/or plasma in infection settings, including hepatitis C and hepatitis B infections [34], [35] as well as during pulmonary tuberculosis [36]. A recent report has also identified miRNAs derived from rice in human serum [37], sparking interest in foreign RNA in body fluids and its diagnostic potential.

Here we investigate the potential of miRNAs to act as tissue and serum biomarkers of experimental mouse and natural human S. mansoni infection. We demonstrate that several host miRNAs are dysregulated in the liver of mice during S. mansoni infection, but do not serve as reliable serum biomarkers of infection in humans. In contrast, we identify at least three parasite-derived miRNAs in the serum of mice infected with S. mansoni that are also detected in human patients and can distinguish 'egg-negative' from 'egg-positive' individuals with high specificity and sensitivity. We anticipate that parasite-derived miRNAs may provide a general platform for specific and non-invasive detection of active helminth infection.

Methods

Animals and S. mansoni infection

For miRNA array analysis of liver and Illumina sequencing of serum in experiment 1, 8-10 week old C57BL/6 mice were left uninfected or infected percutaneously with ~180 S. mansoni cercariae, weighed regularly, and euthanized 7 weeks post infection. For the 12-week time course, 8-10 week old C57BL/6 mice were left uninfected or infected percutaneously with ~80 cercariae, weighed regularly, and euthanized at 4, 6, 8 and 12 weeks post infection.

Collection of blood and serum

Whole blood was drawn from mice by cardiac puncture. The needle was removed before emptying the syringe to avoid haemolysis and blood was allowed to sit for 1 h at RT to clot. Serum was separated by centrifugation at 2500 g for 15 min at 4°C, and the supernatant was collected into a new tube and spun at 10,000 g for 1 min to remove remaining cells. The resultant supernatant was transferred into a new tube and stored at -20°C prior to RNA extraction.

Human serum samples

Human serum was screened retrospectively from samples obtained from schistosome endemic areas in Zimbabwe and Uganda. The study group included in this investigation (Zimbabwe) was part of a larger study on the molecular immunoepidemiology of human schistosomiasis (carried out between November 1999 and March 2000) and had not been included in the National Schistosome Control Programme and therefore had not received treatment for schistosomiasis or other helminth infections. After collection of all samples, all participants were offered anti-helminthic treatment with the recommended dose of praziquantel (40 mg/kg of body weight). The selection of samples for screening here was based on availability of sufficient serum for miRNA analysis, and comprised five 'egg positive' individuals with an average egg per gram in stool (epg) = 108, range: 39-277 and nine 'egg-negative' individuals (Table 1). According to the WHO's classification, Zimbabwe has high to moderate levels of S. haemobium infection (prevalence ranging from 10% to greater than 70%) but low S. mansoni prevalence (less than 10%). The selection of serum for screening from Piida (Uganda) participants was based on availability of sufficient material and comprised twenty individuals infected with S. mansoni with an average epg = 1117 (range: 105-4030) and ten egg-negative individuals. The study group included in this investigation were part of a larger study carried out in Butiaba village, adjacent to Lake Albert, Masindi district, Uganda in 1996, described further in [38]. The cohort had moderate to high S. mansoni infection intensities and prevalence of 91%. After collection of stool and serum samples all the study participants, irrespective of the infection status, received 2 doses of praziquantel, 40 mg/kg of body weight, 6 weeks apart. Efficacy of chemotherapy was assessed 6 weeks after each treatment.

Table 1. Top miRNAs dysregulated in the liver at 7 weeks post S. mansoni infection as determined by microarray analysis and showing results obtained by qRT-PCR (microarray: p<0.05, Fold change ≥2, qRT-PCR: t-test).

Ethics statement

Animal experiments were conducted under a Project License granted by the Home Office (United Kingdom), reference 60/4104, in accordance with local guidelines and approved by the Ethical Review Committee of the University of Edinburgh.

For the human serum samples collected in Chiredzi, permission to conduct the work in this province was obtained from the Provincial Medical Director. Ethical approval was received from the Medical Research Council of Zimbabwe (MRCZ). Only compliant participants were recruited into the study and they were free to drop out at any point during the study. At the beginning of the study, participants and their parents/guardians (in case of children) had the aims and procedures of the project explained fully in the local language, Shona, and oral consent (as was customary) was obtained from participants and parents/guardian before parasitology and blood samples were obtained. For the samples collected in Piida, ethical clearance was obtained from the Uganda National Council of Science and Technology (ethics committee for Vector Control Division, Ugandan Ministry of Health). The aims and procedures were explained to the local community at the start of the study and oral consent was obtained from all adults and from the parents/legal guardians of all children under 15 who were willing to participate. Due to cultural reasons and low levels of literacy, oral consent is deemed acceptable by the Ugandan Ministry of Health and was approved by the Uganda National Council of Science and Technology. Upon oral consent participants were enrolled in the study with a written record of their name, age, sex and case number, this served as both the record of oral consent and enrolment record.

Parasitology and serum sample collection and processing

Participants in the Zimbabwe study were checked for both urogenital and intestinal schistosomiasis and for inclusion in this analysis had to be free of any soil-transmitted helminths and also free of S. haematobium infection to avoid cross reactivity between different helminth parasites. For the analysis of urogenital schistosome infection (S. haematobium) participants submitted three urine and three stool samples (over four consecutive days). 10 ml of each sample received was processed on the day of collection by a urine filtration method [39]. Stool samples were prepared and examined on the day of collection using the Kato-Katz faecal smear for detection of S. mansoni eggs and soil transmitted helminths; this was carried out by trained and experienced technical staff from the National Institute of Health Research [40]. A single slide for microscopic examination was prepared from each urine and stool sample. Serum was prepared from 10 ml venous blood collected from study participants, frozen at -20°C and afterwards stored at -80°C. Samples were transported frozen to Edinburgh and stored at -80°C prior to serological assays. Participants in the Ugandan study had duplicate 50 mg Kato-Katz slides prepared from 3 consecutive stool samples for detection of S. mansoni eggs, expressed as mean epg; this was carried out by trained and experienced technical staff from the Vector Control division at the Ministry of Health in Uganda. Serum was prepared from 10 ml venous blood samples, frozen at -20°C and transported to Cambridge for storage at -80°C prior to serological assays.

Mouse serum RNA extraction

For the 12-week time course experiment, total RNA was extracted from serum using the miRVana PARIS extraction kit (Ambion), according to the manufacturer's protocol. In brief, 100 µl of serum was thawed on ice, mixed with an equal volume of 2× Denaturing Solution and kept on ice for 10 min. Samples were extracted with an equal volume of acid-phenol chloroform, vortexed for 30 s and centrifuged for 10 min at 10,000 g at RT. The aqueous phase was mixed with 1.25 volumes of 100% ethanol and added to the mirVana PARIS column. The column was washed and RNA eluted in 100 µl of 0.1 mM EDTA. RNA samples were stored at -20°C prior to further analysis. Extracted RNA was quantified by Qubit (Invitrogen).

Human serum RNA extraction

Due to small volumes and low amounts of RNA in available human serum samples the extraction protocol was adjusted to result in more concentrated RNA. To do this, 50 µl of serum was thawed on ice, mixed with 50 µl of ddH2O and 100 µl of 2× Denaturing Solution (as supplied in the miRVana Paris kit) and kept on ice for 10 min. Samples were spiked with 10 fmoles of a synthetic RNA, Spike1: 5'-UGCUGAAUGCGUAGCUAUAAGC-3' (IDT) and extracted with an equal volume of acid-phenol chloroform, vortexed for 30 s and centrifuged for 10 min at 10,000 g at RT. The aqueous phase was mixed with 1/10 volume of 3M sodium acetate, 10 µg of GlycoBlue (Ambion) and an equal volume of isopropanol. Samples were allowed to precipitate overnight at -20°C and were then centrifuged at >10,000 g at 4°C. Pellets were washed twice with 75% ethanol, air-dried and then resuspended in 25 µl of 0.1 mM EDTA. The total RNA concentration was below the limit of detection based on Qubit (Invitrogen).

Liver RNA extraction

Liver tissue was immersed in RNA Later Solution (Ambion) overnight at 4°C prior to extraction using TRIzol Reagent (Invitrogen) according to the manufacturer's protocol. RNA was quantified by NanoDrop and integrity assessed by 10% PAGE or Bioanalyzer 2100. All RNA samples used in the microarray analysis had RIN >8.

Reverse transcription & PCR

For reverse transcription of mouse serum samples, a fixed amount of extracted RNA (1.5 ng) was used as an input and 0.1 fmoles of a synthetic RNA, Spike2: 5'-CGUAUCGAGUGAUGUCACGUA-3', was added at the RT step for normalization. For human serum samples, where the total RNA concentration was below the detection limit, a fixed volume of RNA (5 µl), corresponding to 10 µl of extracted serum, was used as the input and Spike1 (added at the time of purification) was used for normalization. For reverse transcription of RNA extracted from liver, 200 ng of total RNA was used in each reaction. Reverse transcription reactions were performed using the miScript System (Qiagen) according to the manufacturer's protocol. PCR was carried out with SYBR green real-time PCR assays (Qiagen) and miScript primers to detect mouse and human miRNAs, according to the manufacturer's protocol (Qiagen). Primers for S. mansoni-specific miRNAs and the synthetic spikes were used at 200 nM final concentration and were purchased from Invitrogen: miR-277, 5'-TAAATGCATTTTCTGGCCCG-3', miR-2162-3p, 5'-TATTATGCAACGTTTCACTCT-3', miR-3479-3p, 5'-TATTGCACTAACCTTCGCCTTG-3', bantam, 5'- TGAGATCGCGATTAAAGCTGGT-3', miR-2a-3p, 5'-TCACAGCCAGTATTGATGAAC-3', miR-71a-3p, 5'- TGAAAGACGATGGTAGTGAGAT-3', sma-miR-n1, 5'-AACTCAGTGGCCTATCGGT-3', sma-miR-n2, 5'-TCAGCTGTGTTCATGTCTTCGA-3', sma-miR-n3, 5'- TGGCGCTTAGTAGAATGTCACCG-3', Spike1, 5'-TGCTGAATGCGTAGCTATAAGC-3', Spike2, 5'-CGTATCGAGTGATGTCACGTA-3'. Data were collected on a Light Cycler 480 System (Roche) with the following temperature profile: pre-denaturation 15 min at 95°C followed by 50 cycles of denaturation for 15 s at 95°C, annealing for 30 s at 55°C, elongation for 30 s at 70°C. The efficiencies of pre-optimized miScript host miRNA probes (Qiagen) were measured from standard curves and ranged between 93-100% and the efficiency of custom probes ranged from 87-97%; both displayed homogenous melting curves and amplification products of the expected size, as described in [41], which were examined here by 6% TBE PAGE (data not shown). Two technical replicates were carried out for each biological replicate. Nuclease free water was used as a non-template control.

miRNA array

For the array analysis, 1 µg of total RNA was labelled using the Hy3 power labeling kit (Exiqon) and hybridized to codelink slides printed with the miRCury 8.1 probe set as described elsewhere [42]. Hybridization and washing were carried out following the manufacturer's protocols (Exiqon). Background signal was subtracted from foreground signal and data were transformed to log-2 scale. Between-array normalization was carried out using global array percentiles (matching median of each array); triplicate probes were represented by the median for each array. Empirical Bayes moderated t statistic (eBayes) was used to test the null hypothesis of "no differential expression" between uninfected (n = 3) and infected (n = 3) liver samples. Since this array was only used as a filter for further validation, the p values were not adjusted for multiple testing. The Exiqon 8.1 arrays contained 384 probes specific for mouse miRNAs.

qRT-PCR data analysis

For analysis of miRNAs in liver samples, the relative fold change between naïve and infected samples was calculated using the 2-ΔΔCt method [43], normalized to miR-16; values for infected mice were compared to values for age-matched naïve controls and the median value for naïve mice was set to 1 for the purpose of calculating fold change. For serum miRNA data analysis, Ct values were "median-normalized" to synthetic RNA spike oligos as described previously [27]: relative change was calculated as 2-Ctn, where Ctn stands for normalized Ct values. The spike-in sequences did not match any known miRNA in miRBase and the primers for detecting these did not yield signals in serum by qRT-PCR, indicating that they do not cross-hybridize with mouse or human small RNAs (data not shown). Fold change was calculated as the ratio of the relative change value of the sample compared to an average of the relative change values of uninfected samples. For the cumulative analysis of miRNAs, the arithmetic mean of fold changes for miR-277, miR-3479-3p and bantam were used.

Statistical analysis of qRT-PCR data and receiver operator characteristic (ROC) curve analysis (95% confidence intervals) was performed with GraphPad Prism (Version 6) software. Two-way ANOVA followed by a Sidak multiple comparison test was used to calculate statistical differences for the mouse miRNA time course data from serum and liver. For the parasite miRNA serum time course, one-way ANOVA followed by Holm-Sidak multiple comparison was used. For analysis of the human serum samples, the Man-Whitney test was used and p-values of <0.05 were considered statistically significant. The measurement of miRNA levels by qRT-PCR was carried out by a trained doctoral student who was not blinded to the results of the infection status of each sample. Of the samples available for screening, none were excluded from the analysis.

Illumina sequencing

For small RNA sequencing in experiment 1, total RNA was extracted from serum of 8 pooled S. mansoni-infected mice (Wk 7, 180 cercariae) and 8 uninfected age matched controls (400 µl total volume) according to the miRVana PARIS protocol (as described above). The small RNAs were size selected by 15% PAGE and prepared according to the Illumina small RNA Sample Preparation Kit version 1.5 and sequenced on the GAIIX. For small RNA sequencing in experiment 2, total RNA was extracted from serum of 3 pooled S. mansoni infected mice (Wk 8, 80 cercariae) and 3 uninfected age matched controls (300 µl total volume) according to the miRVana PARIS protocol. The library was prepared according to the TruSeq Small RNA protocol (without size-selecting small RNA) and sequenced on the HiSeq2.

Raw reads were obtained in fastq format and 3' adapters trimmed using cutadapt, requiring at least a 6 bp match to the adapter sequence and a quality threshold of 20. Only reads that contained the adapter were retained; reads were subsequently collapsed on primary fasta sequence and only reads present at ≥2 copies were analyzed. Trimmed, collapsed reads ≥17 bp were then aligned to mouse (MM9) or S. mansoni genomes (V5.0) using BOWTIE version 0.12.5, requiring a perfect match to the full length of the sequence. Reads that mapped to either genome were then BLASTN aligned against the RFAM database [BLASTN parameters: -max_target_seqs 1 -outfmt '6 std qseq sseq' -task blastn -word_size 6 -dust no] and categorized according to matches to Rfam class (e.g. rRNA, tRNA, etc.). Mouse reads without rfam similarities (other than miRNAs) were aligned to mature miRNAs in miRBase version 19. Some trimmed miRNA reads aligned to more than one family member: the assignment of these ambiguous reads is designated with "X" in Table S3. RNAs that aligned to the S. mansoni genome and did not show RFAM similarities were passed to mirDeep2.0.0.5 using platyhelminth miRNAs from miRBase 19 as guides (Table S4); reads mapping to known miRNAs from miRBase were identified regardless of miRdeep score, for prediction of novel miRNAs a cut-off value of 0 was used for reporting in Table 1.

Results

Specific host miRNAs are significantly dysregulated in the liver of mice at 7 weeks post S. mansoni infection

miRNAs are dysregulated in most disease contexts and play important roles in mediating how cells respond to insult and infection (reviewed in [44]). Approximately 4-6 weeks after S. mansoni infection, mature female parasites produce eggs, some of which are carried by the blood-flow to the liver where they become trapped [45]. The host immune response induced by the presence of the egg antigens leads to the formation of granulomatous lesions, which are composed of immune cells and collagen fibres [46] and result in fibrosis and associated pathology. In order to identify miRNAs associated with S. mansoni-induced liver pathology, and to prioritise candidates for further screening as biomarkers, we first compared expression profiles of miRNAs in livers of naïve mice or mice that were infected with S. mansoni at a high dose (~180 cercariae). Tissues were collected at 7 weeks post infection, at which time substantial granulomas were observed (data not shown). A total of 33 mouse miRNAs were differentially expressed: 26 miRNAs were up-regulated and 7 miRNAs were down-regulated in infected mice, based on a fold change cut-off of ≥2 and p value cut off <0.05 (Table S1, Fig. S1). The miRNAs that displayed the largest differential expression included miR-199a and miR-214, which are known to be altered in liver fibrosis caused by hepatitis C infection or induced by carbon tetrachloride [47], [48]. Among the down regulated miRNAs was the liver-enriched miR-122, which is dysregulated during hepatitis C infection, acetaminophen overdose and hepatocellular carcinoma and is involved in lipid metabolism [49], [50].

For validation of the microarray results, the miRNAs that displayed the largest fold change were quantified by qRT-PCR and normalized to miR-16 (a total of 6 up-regulated miRNAs and 6 down-regulated miRNAs were examined). Consistent with the array results, there was an increase in miR-199-5p, miR-199-3p, miR-214, miR-21, miR-210, and a reduction of miR-192, miR-194, miR-365, miR-122 and miR-151 in the liver tissue of S. mansoni infected mice as compared to naïve mice; miR-9 and miR-744 did not display differential expression and were not analysed further (Table 1). All of these miRNAs are perfectly conserved in mouse and human.

Temporal expression analysis of miR-199, miR-214, miR-21, miR-210, miR-122, miR-192 and miR-194 in the liver during S. mansoni infection

Between weeks 6 and 12, female parasites continue to produce ~300 eggs per day [51], resulting in an increase in the number of granulomas in the liver and the development of fibrosis [45]. For the 10 miRNAs validated to be differentially expressed at 7 weeks post infection, we next examined their temporal expression between 4-12 weeks post infection using a lower parasite dose (80 cercariae). All data from infected mice were compared to age-matched naive mice. To account for differences in RNA extraction or qRT-PCR efficiency, the data were normalised to miR-16, which displayed stable expression in the liver during infection (Fig. S2). Of the 10 miRNAs examined, all except miR-365 and miR-151 were differentially expressed between naïve and infected mice by 6-8 weeks post infection (Fig. 1, Table S2). This timing correlated with the deposition of eggs in the liver, which were detected by 6 weeks post infection and increased by 8 and 12 weeks post infection (Fig. S3). Our results suggest that these cellular miRNAs represent tissue biomarkers of infection that may play a role in the development and progression of liver fibrosis induced by S. mansoni egg deposition.

[Figure omitted, see PDF]

Figure 1. Differential expression of host miRNAs in the livers of mice at 4-12 weeks post infection with S. mansoni.

miRNAs were quantified by qRT-PCR, normalized to miR-16 and fold changes calculated as the ratio of values from infected versus naïve mice (*p<0.05, **p<0.01, ***p<0.001, ****p<0.0001).

Temporal host miRNA dysregulation in the liver during infection is not reflected in serum

Several studies in non-helminth systems have shown that liver-derived miRNAs are detectable in serum and can be used as biomarkers in disease states [52], [53], [54], [55], [56]. We therefore examined whether the murine miRNAs that are altered in liver tissue are similarly altered in serum during S. mansoni infection. qRT-PCR was used to measure miRNA levels in the serum of mice infected with S. mansoni over the 12 week time course. The methodology of miRNA analysis in body fluids is not well standardized, and it is still not clear which small RNAs are appropriate for normalization [57], [58], [59], [60]. Here we used a constant amount of total RNA in the reverse transcription reaction. To account for any variation in qRT-PCR efficiency due to contaminants, miRNA levels were normalized to a synthetic RNA oligo ("spike-in") that was included at the reverse-transcription step. The ratio of miRNA levels in infected versus age-matched naïve mice was quantified at each time point and plotted as fold change (Fig. 2). Compared to the analysis of liver samples, there is more variation in the serum miRNA levels between biological replicates. As shown in Fig. 2, the levels of miR-192, miR-194 and miR-122 in serum do not change between 4-12 weeks post infection, whereas five of the miRNAs that are up-regulated in the liver are also significantly elevated in serum at 12 weeks post infection (p<0.05), ranging from 2.6 fold (miR-21) to 4.7 fold (miR-214) (Table S2). These five host miRNAs represent potential serum biomarkers of S. mansoni infection. However, since their levels do not change until 12 weeks post infection in mice, they may only be useful in cases of more advanced pathology.

[Figure omitted, see PDF]

Figure 2. Differential abundance of host miRNAs in mouse serum during S. mansoni infection.

miRNAs were quantified by qRT-PCR and normalized to a synthetic RNA spike-in. Each symbol represents data from one individual mouse. Fold changes are defined as the ratio of abundance in infected versus naïve serum; the signal from naïve was set as 1. (***p<0.001, ****p<0.0001).

Small RNA sequencing reveals the presence of S. mansoni-derived miRNAs in the serum of infected mice

Multiple cell types release or secrete miRNAs into the circulation [59], [61] and it is possible that other miRNAs, beyond those differentially expressed in the liver, might represent biomarkers of infection. To examine this in an unbiased fashion, small RNAs from serum of naïve mice or mice infected with S. mansoni were sequenced using the Illumina platform. Small RNA libraries were prepared in two independent experiments, using two different preparation methods (Methods). In both experiments, the majority of small RNA reads aligned to the mouse genome (ranging from 63-71%, Table 2). A small percentage of reads unambiguously aligned to the S. mansoni genome in infected mice (0.04-0.14%); <0.01% of reads mapped to the S. mansoni genome in naïve samples. Interestingly, the majority of reads in all samples were derived from tRNAs: 92-98% of the reads that mapped to the mouse genome and 42-100% of reads that mapped to the S. mansoni draft genome (Table 2). It is important to note, however, that due to the sequence similarity of the tRNAs between species and the scope for post-transcriptional editing, we cannot at this point definitively determine the organism from which these derive. Only 1-6% of the reads that mapped to the mouse genome were miRNAs (Table 2, listed in Table S3). There was no correlation between differences in host miRNA levels in naive and infected mice from the two experiments (data not shown), this could relate to the different doses, although it is not possible to derive statistical conclusions from these data due to low read numbers. A total of 78 and 29 mature miRNAs are known to be encoded by the trematodes S. japonicum and S. mansoni, respectively [62], [63], [64], [65], and further miRNAs have been predicted [62]. Analysis of the reads that mapped to the S. mansoni draft genome using the miRdeep2 program [66] identified at least 11 miRNAs from infected samples: 8 of these were positively identified based on identity to known miRNAs in S. mansoni and/or S. japonicum (Table 3) [64], [65] and 3 are predicted by miRdeep (here named sma-miR-n1, sma-miR-n2 and sma-miR-n3). The predicted stem-loop structures for the putative pre-miRNAs are shown in Fig. 3; since the depth of coverage of S. mansoni reads in this study was very low and reads did not map to both arms of the hairpin, these cannot yet be considered prototypical miRNAs. However, sma-miR-n3 shares a seed site with Schmidtea mediterranea miR-2160 and was annotated as sma-miR-8437 in a study of S. mansoni miRNAs published after submission of this manuscript [67]. Additional analysis of the datasets allowing 1 mismatch to the S. mansoni draft genome identified one further miRNA from infected but not naïve serum: miR-277, which is identical to miR-277 in S. japonicum and Echinococcus granulosus but has a C->T mutation at position 17 in relation to the S. mansoni draft genome (Table 3). Other reads from infected samples that aligned to S. mansoni but did not make the miRdeep2 cut-offs are provided in Table S4. It is possible that some of these less abundant sequences derive from real miRNAs but the coverage in these studies is not sufficient to determine this.

[Figure omitted, see PDF]

Figure 3. Predicted stem-loop structures of the S. mansoni miRNAs identified in this study.

The mature sequence that was identified in serum is shown in bold.

Table 2. Small RNA classification in serum of naïve and infected mice.

Table 3. Parasite-miRNAs in serum of mice infected with S. mansoni identified by deep-sequencing.

S. mansoni miRNAs are present in mouse serum as early as 8 weeks post infection

To determine the kinetic profile of parasite miRNAs in serum during S. mansoni infection, qRT-PCR analysis was carried out as described above, using primers specific for 9 of the 11 parasite miRNAs. The other 2 miRNAs, sma-miR-10-5p and sma-let-7-3p, were excluded from analysis because they are highly similar to homologous mouse miRNAs that are present at >100 fold higher read frequencies (Table S3). Importantly, most miRNAs in helminth parasites have evolved after the last common ancestor with their vertebrate hosts and are therefore distinguishable in sequence, however several miRNAs are perfectly conserved across animals or highly similar in sequence [63]. Several of the parasite-specific probes showed a signal in the serum of naïve mice, which is presumably due to cross-hybridization with endogenous small RNAs. In cases where no signal was observed in naïve mice, the maximum cycle value of 50 was set as background for the purpose of calculating signal over noise (which we interchange here with "fold change"). Six of the nine parasite miRNA probes tested (miR-277, bantam, miR-3479-3p, miR-2a-3p, miR-n1, miR-n2) showed a statistically significant signal over noise at 8 or 12 weeks post infection (Fig. 4, p<0.05); the three miRNAs that were not reliably detected (miR-n3, miR-71a-3p, miR-2162-3p) were not analysed further. The average signal over noise ratios for each probe during the time course of infection are provided in Table S5 and range from 4.2 to >3,000.

[Figure omitted, see PDF]

Figure 4. Detection of parasite-derived miRNAs in mouse serum during S. mansoni infection.

miRNAs were quantified by qRT-PCR, normalized to a synthetic RNA spike-in and fold change calculated as the ratio of abundance in infected serum compared to the background abundance level detected in naïve serum, which represents the noise in the assay likely derived from cross-hybridization with endogenous small RNAs (*p<0.05, **p<0.01, ***p<0.001, ****p<0.0001).

Three S. mansoni-derived miRNAs are detected in human serum and distinguish egg-negative from egg-positive individuals

Based on the results described above, we extended our analyses to human patients, using the 6 parasite miRNAs and 5 mouse miRNAs that displayed differential abundance in serum of infected compared to naïve mice. Serum samples from two field sites were examined: an area of high infection in the Piida community of Uganda and an area of low infection in the Chiredzi community of Zimbabwe. The high infection samples were collected from mixed-age participants from Piida diagnosed as S. mansoni infected, termed 'egg positive' and compared to volunteers from the same community with undetectable parasite eggs in the stool, termed 'egg-negative'. The low infection samples were collected from children in Chiredzi diagnosed with S. mansoni and compared to age matched participants with undetectable eggs by standard stool examination methods. Demographic data of individuals are provided in Tables 4 and 5. The signal over noise was calculated as described above, using synthetic spike-ins for normalization. miR-n1, miR-n2 and miR-2a-3p were below the detection limit (Ct = 50) in both 'egg-positive' and 'egg-negative' samples. The 5 host miRNAs were detectable in serum (miR-21, miR-199-3p, miR-199-5p, miR-210, miR-214) but showed variable abundance and failed to differentiate 'egg-positive' and 'egg-negative' participants (Fig. S4). In contrast, three out of the six parasite miRNAs (bantam, miR-277 and miR-3479-3p) displayed a significant signal over noise level in the serum of S. mansoni infected individuals 'egg-positive' from both high (Fig. 5A) and low (Fig. 5B) infection endemic areas, compared to the 'egg-negative' participants from the same communities (p<0.05, Mann - Whitney test). Data presented using ROC curves show that the single parasite miRNAs discriminated between S. mansoni 'egg-negative' and 'egg-positive' with an area under the curve (AUC) of 0.785, 0.790, 0.768 for bantam, miR-277 and miR-3479-3p, respectively, in the individuals from Uganda (Fig. 5A) and 0.889, 0.933, 0.911 in the individuals from Zimbabwe (Fig. 5B). Using optimal cut-off points, this translates to detection of S. mansoni infected individuals with specificity/sensitivity of 80%/60%, 80%/70% and 80%/60%, in the patients from Uganda (Fig. 5A) and specificity/sensitivity of 100%/60%, 89%/80% and 89%/80% respectively in the patients from Zimbabwe (Fig. 5B). A repeated measurement of the parasite miRNA levels in the same samples displayed Pearson correlation values between 0.86 to 0.95 and comparable specificity/sensitivity values (Table S6). When combining the data for all three of the miRNAs into a cumulative value, the AUC increased to 0.845 and 0.933 for each cohort of participants. This resulted in improved specificity/sensitivity in detection for samples from Uganda (80%/90%) using a fold change cut-off of 1.189 (Fig. 6). This approach is very similar to a recent report that uses a "miRNA score" based on cumulative normalized signals of miRNAs [68]; analysis of our data based on cumulative fold change or cumulative normalized signals yields very similar results (Fig. S5). These results show that combining data for bantam, miR-277 and miR-3479-3p may improve sensitivity of S. mansoni diagnosis compared to analysis of individual miRNAs.

[Figure omitted, see PDF]

Figure 5. Discrimination between the S. mansoni infected and uninfected individuals using parasite-specific miRNA detection.

miRNAs were quantified by qRT-PCR, normalized to a synthetic RNA spike-in and fold changes calculated as the ratio of infected to uninfected (median with interquartile range indicated). Piida- panel A, Chiredzi- panel B. Specificity and sensitivity and ROC curves with AUC are indicated. (*p<0.05, **p<0.01).

[Figure omitted, see PDF]

Figure 6. Discrimination between the S. mansoni infected and uninfected individuals using the combined data for 3 parasite miRNAs.

Fold-changes (median with interquartile range indicated) and ROC curves for bantam, miR-277 and miR-3479-3p (Piida- panel A, Chiredzi-panel B) (**p<0.01).

Table 4. Piida study participants.

Table 5. Chiredzi study participants.

Discussion

The recent evidence that miRNAs can be released into circulation from mammalian cells and tissues has stimulated extensive interest in the potential use of these molecules as non-invasive biomarkers [31], [32], [57]. MiRNA-based diagnostics are being developed for a number of diseases and although qRT-PCR is the most common detection method at present, there is extensive interest in improving and diversifying detection technologies, which may provide more field-friendly tools. Since miRNAs have been shown to be extremely stable in body fluids [27], [30], [69], we anticipate that these nucleic acids could be particularly useful as diagnostics in field settings where collection and storage conditions can be difficult to control. Here we find that miRNAs derived from the helminth parasite S. mansoni are present in infected mouse and human serum and offer advantages over endogenous miRNAs as biomarkers of infection. Specifically, we find 9 known miRNAs and at least 2 novel putative miRNAs derived from S. mansoni that are present in the serum of infected mice. The read counts of some of these are very low (between 2-66 reads per 930,209 mouse miRNAs reads) and will require further validation with better coverage. However, we show that three of these miRNAs (bantam, miR-277, mirR-3479-3p) can be detected in human serum from schistosome endemic areas. These represent a direct marker for infection and may also provide an indirect marker for the pathology induced by infection. Notably, a study published after submission of this manuscript identified 5 miRNAs derived from S. japonicum in the plasma of infected rabbits and 3 of these are identical or homologous to those identified here: bantam, miR-3479-3p and miR-10-5p [70], providing independent validation for the presence of trematode miRNAs in the serum of infected animals. Here we demonstrate small RNAs derived from a helminth parasite are also present in patient serum, and provides a starting point for developing more field-friendly methods for their detection. This work also extends a burgeoning area of research detailing "foreign" small RNAs in body fluids. Zhang et al., (2012) recently reported that the plant miRNA, miR-168a, is present in human and animal serum and demonstrated that this derives from a rice diet [37]. Wang et al., (2012) recently reported that a wide range of exogenous RNAs can be found in human plasma, including small RNAs derived from bacteria and fungi [71]. From our results is not yet clear whether the miRNAs identified in serum are actively secreted from the parasite and how they are stabilized, since serum contains high amounts of RNases [72], [73]. From the time course analysis presented here, schistosome-derived small RNAs were reliably detected in serum by 8 weeks post infection, after deposition of the eggs in the liver (Fig. 4 and Fig. S4). At present, we cannot determine whether these RNAs derive from the adult worms or the eggs; further in vitro and in vivo studies will shed light on this issue, which is relevant to diagnostic applications, for example the capacity to detect pre-patent infection.

It is intriguing to think that existence of these miRNAs outside the parasite has a function, but this is beyond the scope of the current analysis. Given the short size of miRNAs and the fact that they do not require perfect complementarity with their targets, one could predict hundreds of possible targets in the host. Interestingly, Xue et al. (2008) showed that three of the miRNAs that we find in serum (sja-bantam, sja-miR-71 and sja-let-7) are expressed during all the stages of parasite development but are enriched in the cercariae, suggesting that they may be important during the initial stages of schistosome infection [65]. The bantam miRNA has been implicated in regulating organ growth in response to environmental conditions in Drosophila as well as C. elegans [74], [75] but functional homologues of bantam do not exist in mammals.

An initial objective in this study was to identify host miRNAs that may be involved in liver pathology associated with S. mansoni infection and to then determine whether these hold any diagnostic value. A number of reports have demonstrated an increase in miR-122 and miR-192 in plasma or serum upon viral infection as well as chemically induced liver disease [54], [56]. However, according to our analysis, although miR-192, miR-122 and miR-194 were down-regulated in the liver during infection, their levels in serum did not change significantly (Fig. 1-2). In contrast, the miRNAs up-regulated in the liver (miR-199-3p, miR-199-5p, miR-21, miR-214 and miR-210) showed significantly higher levels in mouse serum at 12 weeks post infection (Fig. 2), however these failed to differentiate S. mansoni infected from uninfected humans (Fig. S4). It should be noted that in addition to significant liver disease, pulmonary and intestinal complications can also occur during S. mansoni infection that could also contribute to serum miRNA levels [45], [76]. Related to this, a possible limitation in the use of endogenous miRNAs as biomarkers is the fact that their differential abundance in serum can derive from multiple cell types and can also be attributed to unrelated conditions [59]. Although the diagnostic value of these host miRNAs in serum is therefore not obvious, our work provides a foundation for further research into the functional role of these miRNAs in S. mansoni pathogenesis. Interestingly, our results do not overlap with those reported by Han et al., (2013) who examined changes in host miRNA levels in the liver of BALB/c mice during S. japonicum infection. We assume this may be due to the very early time point (10 days post infection) used in their study. Indeed, it is likely that the host miRNA changes we observe are primarily related to the liver pathology and/or the immune response initiated by the parasite eggs trapped in the liver, rather than the initial host immune response to the schistosomula. Several of the miRNAs we identify as differentially expressed are already known to be associated with liver fibrosis or disease in other, non-helminth, settings [49], [50], [77], [78]. The functional role of the miRNAs in the pathology induced by S. mansoni infection remains to be determined. Future work in this area will shed light on the molecular basis of pathology and may offer innovative new therapeutic strategies.

Importantly, we report here that parasite-derived miRNAs can be detected in human serum and can distinguish 'egg-negative' from 'egg-positive' individuals in areas of both low and high infection intensity. By combining data for miR-277, miR-3479-3p and bantam we detected infection with a sensitivity of 80-90% and specificity of 80-89%. We anticipate that this may be improved further by optimizing isolation protocols, probe design and more robust methods for normalizing the data, for example to identify appropriate endogenous small RNAs that could be used as controls, rather than synthetic spike-ins [79]. The main limitation of miRNA detection in serum appears to be cross-hybridization of the probes with endogenous small RNAs, which influences the signal to noise ratios. In particular, we report the finding that the majority of small RNAs in mouse serum are derived from tRNAs (Table 2), consistent with a recent study [80]. It is possible that depletion of these tRNAs prior to qRT-PCR may improve the specificity or sensitivity of miRNA detection; this requires further investigation. In addition, one of the parasite-derived miRNAs identified by sequencing, miR-2162-3p could not be validated by qRT-PCR, likely owing to its low GC content (33%). Optimization of probe design for the parasite-derived miRNAs may also greatly increase their diagnostic utility in larger scale studies in human patients. On this note, the work presented here was performed on a relatively small number of individuals. Following further optimization of the extraction methods and probe design to minimize cross hybridization and sensitivity, larger studies will be important to assess the full potential of the proposed miRNA biomarkers, including positive and negative predictive values in comparison to existing techniques. It will also be of interest to determine the origin of these parasite-derived miRNAs and examine their abundance levels in response to treatment. We anticipate that the parasite miRNAs in serum could complement or transform existing diagnostic strategies and may serve as a platform for detecting a range of helminth infections.

Supporting Information

Figure S1.

Volcano plot of miRNA microarray data comparing miRNA expression levels in mice infected with S. mansoni compared to naive mice. The horizontal line indicates a p-value cut-off of 0.05 and vertical lines indicate fold change thresholds of 2.

(EPS)

Figure S2.

Expression of miR-16 in liver and serum over the time course of S. mansoni infection. miRNAs were quantified by qRT-PCR, relative change was calculated as 2-Ct (for liver samples) or 2-Ctn,(for serum samples) and fold changes were calculated as infected to naïve ratio.

(EPS)

Figure S3.

S. mansoni egg counts in the livers of mice over the course of infection.

(EPS)

Figure S4.

Abundance of host miRNAs in serum of S. mansoni egg-positive and egg-negative individuals. miRNAs were quantified by qRT-PCR, normalized to a synthetic RNA spike-in and fold changes calculated as the ratio of infected to uninfected individuals. Piida- panel A, Chiredzi- panel B.

(EPS)

Figure S5.

Discrimination between the S. mansoni infected and uninfected patients using the miRNA-score. Dot plot comparing cumulative (bantam, miR-277 and miR-3479-3p) miRNA-score between S. mansoni egg positive and negative participants (median with interquartile range indicated) and ROC curves (Piida- panel A, Chiredzi-panel B) (**p<0.01).

(EPS)

Table S1.

miRNAs that are dysregulated in the liver upon S.mansoni infection as determined by microarray analysis (p<0.05, fold change ≥2).

(DOCX)

Table S2.

Relative expression of miRNAs in the liver and serum during the time course of S.mansoni infection based on qRT-PCR analysis, normalized to values in naïve mice.

(DOCX)

Table S3.

Tabulated mouse miRNA reads identified in serum of naïve and infected mice. The "X" denotes the miRNAs who were assigned reads that could have derived from one of multiple miRNA family members (e.g. in the case of shortened reads).

(XLSX)

Table S4.

Prediction of known or novel S. mansoni miRNAs in mouse serum based on miRdeep2. This table shows all predictions (only predictions with a miRdeep2 score >0 are included in the text). Read assignments are: uninfected experiment 1 (R1N), infected experiment 1 (R1I), uninfected experiment 2 (607), infected experiment 2 (608 or 609). The number of reads in that sample is listed after the "x".

(TXT)

Table S5.

Relative expression of parasite miRNAs in serum during the time course of S.mansoni infection, based on qRT-PCR analysis, normalized to values in naïve mice.

(DOCX)

Table S6.

miRNA detection values in Chiredzi and Uganda samples. Fold change is defined as the abundance value in each sample compared to the median of values obtained in uninfected individuals. Repeated miRNA measurements from Chiredzi samples are shown (there was not enough material to repeat measurement in Uganda samples).

(XLSX)

Table S7.

Ct values obtained for qRT-PCR.

(XLSX)

Checklist S1.

STARD Checklist.

(DOC)

Acknowledgments

We greatly acknowledge the help of Rinku Rajan with the mouse S. mansoni infections and Fabio Simbari, Nouf Laqtom, Gillian Coakley, Marissa Lear with organ harvests. We thank Sujai Kumar for processing initial sequencing data and Dan Frank and Margo Chase-Topping for helpful discussions on analysis. Schistosome life cycle stages used to infect mice for this research were supplied by the N.I.A.I.D. Schistosomiasis Resource Centre at the Biomedical Research Institute (Rockville, MD) through Contract N01-AI-30026.

The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE 49668.

Author Contributions

Conceived and designed the experiments: AMH RJL ASM AHB. Performed the experiments: AMH RJL JFQ AHB. Analyzed the data: AMH TF AI AHB. Contributed reagents/materials/analysis tools: NN FM NBK FJ DWD. Wrote the paper: AMH AHB.

References

Hotez PJ, Brindley PJ, Bethony JM, King CH, Pearce EJ, et al. (2008) Helminth infections: the great neglected tropical diseases. The Journal of clinical investigation 118: 1311-1321. doi: 10.1172/jci34261

Awasthi S, Bundy DA, Savioli L (2003) Helminthic infections. BMJ 327: 431-433. doi: 10.1136/bmj.327.7412.431

Organisation WH (2013) Schistosomiasis-Fact sheet 115. WHO

van der Werf MJ, de Vlas SJ, Brooker S, Looman CW, Nagelkerke NJ, et al. (2003) Quantification of clinical morbidity associated with schistosome infection in sub-Saharan Africa. Acta tropica 86: 125-139

King CH, Dickman K, Tisch DJ (2005) Reassessment of the cost of chronic helmintic infection: a meta-analysis of disability-related outcomes in endemic schistosomiasis. Lancet 365: 1561-1569. doi: 10.1016/s0140-6736(05)66457-4

Organization WH (2002) Prevention and Control of Schistosomiasia and Soil Transmitted helminthiasia- Report of a WHO Expert Committee. Geneva

Glinz D, Silue KD, Knopp S, Lohourignon LK, Yao KP, et al. (2010) Comparing diagnostic accuracy of Kato-Katz, Koga agar plate, ether-concentration, and FLOTAC for Schistosoma mansoni and soil-transmitted helminths. PLoS neglected tropical diseases 4: e754. doi: 10.1371/journal.pntd.0000754

Stothard JR, Sousa-Figueiredo JC, Betson M, Adriko M, Arinaitwe M, et al. (2011) Schistosoma mansoni Infections in young children: when are schistosome antigens in urine, eggs in stool and antibodies to eggs first detectable? PLoS neglected tropical diseases 5: e938. doi: 10.1371/journal.pntd.0000938

Doenhoff MJ, Chiodini PL, Hamilton JV (2004) Specific and sensitive diagnosis of schistosome infection: can it be done with antibodies? Trends in Parasitology 20: 35-39. doi: 10.1016/j.pt.2003.10.019

Shane HL, Verani JR, Abudho B, Montgomery SP, Blackstock AJ, et al. (2011) Evaluation of urine CCA assays for detection of Schistosoma mansoni infection in Western Kenya. PLoS neglected tropical diseases 5: e951. doi: 10.1371/journal.pntd.0000951

van Lieshout L, De Jonge N, Mansour MM, Bassily S, Krijger FW, et al. (1993) Circulating cathodic antigen levels in serum and urine of schistosomiasis patients before and after chemotherapy with praziquantel. Transactions of the Royal Society of Tropical Medicine and Hygiene 87: 311-312. doi: 10.1016/0035-9203(93)90142-d

van Dam GJ, Bogitsh BJ, van Zeyl RJ, Rotmans JP, Deelder AM (1996) Schistosoma mansoni: in vitro and in vivo excretion of CAA and CCA by developing schistosomula and adult worms. The Journal of parasitology 82: 557-564. doi: 10.2307/3283780

Coulibaly JT, N'Gbesso YK, Knopp S, N'Guessan NA, Silue KD, et al. (2013) Accuracy of urine circulating cathodic antigen test for the diagnosis of Schistosoma mansoni in preschool-aged children before and after treatment. PLoS neglected tropical diseases 7: e2109. doi: 10.1371/journal.pntd.0002109

Sousa-Figueiredo JC, Betson M, Kabatereine NB, Stothard JR (2013) The urine circulating cathodic antigen (CCA) dipstick: a valid substitute for microscopy for mapping and point-of-care diagnosis of intestinal schistosomiasis. PLoS neglected tropical diseases 7: e2008. doi: 10.1371/journal.pntd.0002008

van Dam GJ, de Dood CJ, Lewis M, Deelder AM, van Lieshout L, et al. (2013) A robust dry reagent lateral flow assay for diagnosis of active schistosomiasis by detection of Schistosoma circulating anodic antigen. Exp Parasitol 135: 274-282. doi: 10.1016/j.exppara.2013.06.017

Enk MJ, Oliveira e Silva G, Rodrigues NB (2012) Diagnostic accuracy and applicability of a PCR system for the detection of Schistosoma mansoni DNA in human urine samples from an endemic area. PloS one 7: e38947. doi: 10.1371/journal.pone.0038947

Pontes LA, Oliveira MC, Katz N, Dias-Neto E, Rabello A (2003) Comparison of a polymerase chain reaction and the Kato-Katz technique for diagnosing infection with Schistosoma mansoni. The American journal of tropical medicine and hygiene 68: 652-656

Wichmann D, Panning M, Quack T, Kramme S, Burchard GD, et al. (2009) Diagnosing schistosomiasis by detection of cell-free parasite DNA in human plasma. PLoS neglected tropical diseases 3: e422. doi: 10.1371/journal.pntd.0000422

Lodh N, Mwansa JC, Mutengo MM, Shiff CJ (2013) Diagnosis of Schistosoma mansoni without the Stool: Detecting DNA from Filtered Urine Comparison of Three Diagnostic Tests to Detect Schistosoma mansoni Infection from Filtered Urine in Zambia. Am J Trop Med Hyg 89: 46-50. doi: 10.4269/ajtmh.13-0104

Ibironke OA, Phillips AE, Garba A, Lamine SM, Shiff C (2011) Diagnosis of Schistosoma haematobium by detection of specific DNA fragments from filtered urine samples. The American journal of tropical medicine and hygiene 84: 998-1001. doi: 10.4269/ajtmh.2011.10-0691

Bartel DP (2009) MicroRNAs: target recognition and regulatory functions. Cell 136: 215-233. doi: 10.1016/j.cell.2009.01.002

Fabian MR, Sonenberg N (2012) The mechanics of miRNA-mediated gene silencing: a look under the hood of miRISC. Nature structural & molecular biology 19: 586-593. doi: 10.1038/nsmb.2296

Sullivan CS, Ganem D (2005) MicroRNAs and viral infection. Molecular cell 20: 3-7. doi: 10.1016/j.molcel.2005.09.012

Xiao C, Rajewsky K (2009) MicroRNA control in the immune system: basic principles. Cell 136: 26-36. doi: 10.1016/j.cell.2008.12.027

Mendell JT, Olson EN (2012) MicroRNAs in stress signaling and human disease. Cell 148: 1172-1187. doi: 10.1016/j.cell.2012.02.005

Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, et al. (2007) Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nature cell biology 9: 654-659. doi: 10.1038/ncb1596

Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, et al. (2008) Circulating microRNAs as stable blood-based markers for cancer detection. Proceedings of the National Academy of Sciences of the United States of America 105: 10513-10518. doi: 10.1073/pnas.0804549105

Xu L, Yang BF, Ai J (2013) MicroRNA transport: a new way in cell communication. Journal of cellular physiology 228: 1713-1719. doi: 10.1002/jcp.24344

Creemers EE, Tijsen AJ, Pinto YM (2012) Circulating microRNAs: novel biomarkers and extracellular communicators in cardiovascular disease? Circulation research 110: 483-495. doi: 10.1161/circresaha.111.247452

Chen X, Ba Y, Ma L, Cai X, Yin Y, et al. (2008) Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell research 18: 997-1006. doi: 10.1038/cr.2008.282

Etheridge A, Lee I, Hood L, Galas D, Wang K (2011) Extracellular microRNA: A new source of biomarkers. Mutation Research-Fundamental and Molecular Mechanisms of Mutagenesis 717: 85-90. doi: 10.1016/j.mrfmmm.2011.03.004

Kosaka N, Iguchi H, Ochiya T (2010) Circulating microRNA in body fluid: a new potential biomarker for cancer diagnosis and prognosis. Cancer Science 101: 2087-2092. doi: 10.1111/j.1349-7006.2010.01650.x

Weiland M, Gao XH, Zhou L, Mi QS (2012) Small RNAs have a large impact Circulating microRNAs as biomarkers for human diseases. Rna Biology 9: 850-859. doi: 10.4161/rna.20378

Shrivastava S, Petrone J, Steele R, Lauer GM, Bisceglie AM, et al. (2013) Upregulation of circulating miR-20a is correlated with hepatitis C virus mediated liver disease progression. Hepatology 58: 863-71. doi: 10.1002/hep.26296

Li LM, Hu ZB, Zhou ZX, Chen X, Liu FY, et al. (2010) Serum microRNA profiles serve as novel biomarkers for HBV infection and diagnosis of HBV-positive hepatocarcinoma. Cancer research 70: 9798-9807. doi: 10.1158/0008-5472.can-10-1001

Qi Y, Cui L, Ge Y, Shi Z, Zhao K, et al. (2012) Altered serum microRNAs as biomarkers for the early diagnosis of pulmonary tuberculosis infection. BMC infectious diseases 12: 384. doi: 10.1186/1471-2334-12-384

Zhang L, Hou D, Chen X, Li D, Zhu L, et al. (2012) Exogenous plant MIR168a specifically targets mammalian LDLRAP1: evidence of cross-kingdom regulation by microRNA. Cell research 22: 107-126. doi: 10.1038/cr.2011.158

Kabatereine NB, Vennervald BJ, Ouma JH, Kemijumbi J, Butterworth AE, et al. (1999) Adult resistance to schistosomiasis mansoni: age-dependence of reinfection remains constant in communities with diverse exposure patterns. Parasitology 118 (Pt 1) 101-105. doi: 10.1017/s0031182098003576

Mott KE (1983) A reusable polyamide filter for diagnosis of S. haematobium infection by urine filtration. Bulletin de la Societe de pathologie exotique et de ses filiales 76: 101-104

Katz N, Chaves A, Pellegrino J (1972) A simple device for quantitative stool thick-smear technique in Schistosomiasis mansoni. Revista do Instituto de Medicina Tropical de Sao Paulo 14: 397-400

Buck AH, Santoyo-Lopez J, Robertson KA, Kumar DS, Reczko M, et al. (2007) Discrete clusters of virus-encoded micrornas are associated with complementary strands of the genome and the 7.2-kilobase stable intron in murine cytomegalovirus. J Virol 81: 13761-13770. doi: 10.1128/jvi.01290-07

Ruckerl D, Jenkins SJ, Laqtom NN, Gallagher IJ, Sutherland TE, et al. (2012) Induction of IL-4Ralpha-dependent microRNAs identifies PI3K/Akt signaling as essential for IL-4-driven murine macrophage proliferation in vivo. Blood 120: 2307-2316. doi: 10.1182/blood-2012-02-408252

Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods 25: 402-408. doi: 10.1006/meth.2001.1262

Mendell JT, Olson EN (2012) MicroRNAs in Stress Signaling and Human Disease. Cell 148: 1172-1187. doi: 10.1016/j.cell.2012.02.005

Pearce EJ, MacDonald AS (2002) The immunobiology of schistosomiasis. Nature reviews Immunology 2: 499-511. doi: 10.1038/nri843

Wynn TA, Thompson RW, Cheever AW, Mentink-Kane MM (2004) Immunopathogenesis of schistosomiasis. Immunological reviews 201: 156-167. doi: 10.1111/j.0105-2896.2004.00176.x

Iizuka M, Ogawa T, Enomoto M, Motoyama H, Yoshizato K, et al. (2012) Induction of microRNA-214-5p in human and rodent liver fibrosis. Fibrogenesis & tissue repair 5: 12. doi: 10.1186/1755-1536-5-12

Murakami Y, Toyoda H, Tanaka M, Kuroda M, Harada Y, et al. (2011) The progression of liver fibrosis is related with overexpression of the miR-199 and 200 families. PloS one 6: e16081. doi: 10.1371/journal.pone.0016081

Girard M, Jacquemin E, Munnich A, Lyonnet S, Henrion-Caude A (2008) miR-122, a paradigm for the role of microRNAs in the liver. Journal of hepatology 48: 648-656. doi: 10.1016/j.jhep.2008.01.019

Hu J, Xu Y, Hao J, Wang S, Li C, et al. (2012) MiR-122 in hepatic function and liver diseases. Protein & cell 3: 364-371. doi: 10.1007/s13238-012-2036-3

Moore DV, Sandground JH (1956) The relative egg producing capacity of Schistosoma mansoni and Schistosoma japonicum. The American journal of tropical medicine and hygiene 5: 831-840

Waidmann O, Koberle V, Brunner F, Zeuzem S, Piiper A, et al. (2012) Serum microRNA-122 predicts survival in patients with liver cirrhosis. PloS one 7: e45652. doi: 10.1371/journal.pone.0045652

Farid WR, Pan Q, van der Meer AJ, de Ruiter PE, Ramakrishnaiah V, et al. (2012) Hepatocyte-derived microRNAs as serum biomarkers of hepatic injury and rejection after liver transplantation. Liver transplantation : official publication of the American Association for the Study of Liver Diseases and the International Liver Transplantation Society 18: 290-297. doi: 10.1002/lt.22438

Starkey Lewis PJ, Merz M, Couttet P, Grenet O, Dear J, et al. (2012) Serum microRNA biomarkers for drug-induced liver injury. Clinical pharmacology and therapeutics 92: 291-293

Wang LG, Gu J (2012) Serum microRNA-29a is a promising novel marker for early detection of colorectal liver metastasis. Cancer epidemiology 36: e61-67. doi: 10.1016/j.canep.2011.05.002

Zhang Y, Jia Y, Zheng R, Guo Y, Wang Y, et al. (2010) Plasma microRNA-122 as a biomarker for viral-, alcohol-, and chemical-related hepatic diseases. Clinical chemistry 56: 1830-1838. doi: 10.1373/clinchem.2010.147850

Weiland M, Gao XH, Zhou L, Mi QS (2012) Small RNAs have a large impact: circulating microRNAs as biomarkers for human diseases. RNA biology 9: 850-859. doi: 10.4161/rna.20378

Scholer N, Langer C, Dohner H, Buske C, Kuchenbauer F (2010) Serum microRNAs as a novel class of biomarkers: a comprehensive review of the literature. Experimental hematology 38: 1126-1130. doi: 10.1016/j.exphem.2010.10.004

Hoy AM, Buck AH (2012) Extracellular small RNAs: what, where, why? Biochemical Society transactions 40: 886-890. doi: 10.1042/bst20120019

Benz F, Roderburg C, Vargas Cardenas D, Vucur M, Gautheron J, et al. (2013) U6 is unsuitable for normalization of serum miRNA levels in patients with sepsis or liver fibrosis. Experimental & molecular medicine 45: e42. doi: 10.1038/emm.2013.81

Turchinovich A, Weiz L, Burwinkel B (2012) Extracellular miRNAs: the mystery of their origin and function. Trends in biochemical sciences 37: 460-465. doi: 10.1016/j.tibs.2012.08.003

de Souza Gomes M, Muniyappa MK, Carvalho SG, Guerra-Sa R, Spillane C (2011) Genome-wide identification of novel microRNAs and their target genes in the human parasite Schistosoma mansoni. Genomics 98: 96-111. doi: 10.1016/j.ygeno.2011.05.007

Kozomara A, Griffiths-Jones S (2011) miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic acids research 39: D152-157. doi: 10.1093/nar/gkq1027

Wang Z, Xue X, Sun J, Luo R, Xu X, et al. (2010) An "in-depth" description of the small non-coding RNA population of Schistosoma japonicum schistosomulum. PLoS neglected tropical diseases 4: e596. doi: 10.1371/journal.pntd.0000596

Xue X, Sun J, Zhang Q, Wang Z, Huang Y, et al. (2008) Identification and characterization of novel microRNAs from Schistosoma japonicum. PloS one 3: e4034. doi: 10.1371/journal.pone.0004034

Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N (2012) miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic acids research 40: 37-52. doi: 10.1093/nar/gkr688

Marco A, Kozomara A, Hui JH, Emery AM, Rollinson D, et al. (2013) Sex-biased expression of microRNAs in Schistosoma mansoni. PLoS neglected tropical diseases 7: e2402. doi: 10.1371/journal.pntd.0002402

Goren Y, Kushnir M, Zafrir B, Tabak S, Lewis BS, et al. (2012) Serum levels of microRNAs in patients with heart failure. European journal of heart failure 14: 147-154. doi: 10.1093/eurjhf/hfr155

Li Y, Jiang Z, Xu L, Yao H, Guo J, et al. (2011) Stability analysis of liver cancer-related microRNAs. Acta biochimica et biophysica Sinica 43: 69-78

Cheng G, Luo R, Hu C, Cao J, Jin Y (2013) Deep sequencing-based identification of pathogen-specific microRNAs in the plasma of rabbits infected with Schistosoma japonicum. Parasitology 140: 1751-1761. doi: 10.1017/s0031182013000917

Wang K, Li H, Yuan Y, Etheridge A, Zhou Y, et al. (2012) The complex exogenous RNA spectra in human plasma: an interface with human gut biota? PloS one 7: e51009. doi: 10.1371/journal.pone.0051009

Blank A, Dekker CA (1981) Ribonucleases of human serum, urine, cerebrospinal fluid, and leukocytes. Activity staining following electrophoresis in sodium dodecyl sulfate-polyacrylamide gels. Biochemistry 20: 2261-2267. doi: 10.1021/bi00511a030

Kamm RC, Smith AG (1972) Ribonuclease activity in human plasma. Clinical biochemistry 5: 198-200. doi: 10.1016/s0009-9120(72)80033-x

Brennecke J, Hipfner DR, Stark A, Russell RB, Cohen SM (2003) bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell 113: 25-36. doi: 10.1016/s0092-8674(03)00231-9

Schaedel ON, Gerisch B, Antebi A, Sternberg PW (2012) Hormonal signal amplification mediates environmental conditions during development and controls an irreversible commitment to adulthood. PLoS biology 10: e1001306. doi: 10.1371/journal.pbio.1001306

Crosby A, Jones FM, Kolosionek E, Southwood M, Purvis I, et al. (2011) Praziquantel reverses pulmonary hypertension and vascular remodeling in murine schistosomiasis. Am J Respir Crit Care Med 184: 467-473. doi: 10.1164/rccm.201101-0146oc

Chen XM (2009) MicroRNA signatures in liver diseases. World journal of gastroenterology : WJG 15: 1665-1672. doi: 10.3748/wjg.15.1665

Jiang X, Tsitsiou E, Herrick SE, Lindsay MA (2010) MicroRNAs and the regulation of fibrosis. The FEBS journal 277: 2015-2021. doi: 10.1111/j.1742-4658.2010.07632.x

Kang K, Peng X, Luo J, Gou D (2012) Identification of circulating miRNA biomarkers based on global quantitative real-time PCR profiling. Journal of animal science and biotechnology 3: 4. doi: 10.1186/2049-1891-3-4

Dhahbi JM, Spindler SR, Atamna H, Yamakawa A, Boffelli D, et al. (2013) 5' tRNA halves are present as abundant complexes in serum, concentrated in blood cells, and modulated by aging and calorie restriction. BMC genomics 14: 298. doi: 10.1186/1471-2164-14-298