- Research article
- Open Access
A novel member of the let-7 microRNA family is associated with developmental transitions in filarial nematode parasites
© Winter et al.; licensee BioMed Central. 2015
- Received: 15 January 2015
- Accepted: 13 April 2015
- Published: 22 April 2015
Filarial nematodes are important pathogens in the tropics transmitted to humans via the bite of blood sucking arthropod vectors. The molecular mechanisms underpinning survival and differentiation of these parasites following transmission are poorly understood. microRNAs are small non-coding RNA molecules that regulate target mRNAs and we set out to investigate whether they play a role in the infection event.
microRNAs differentially expressed during the early post-infective stages of Brugia pahangi L3 were identified by microarray analysis. One of these, bpa-miR-5364, was selected for further study as it is upregulated ~12-fold at 24 hours post-infection, is specific to clade III nematodes, and is a novel member of the let-7 family, which are known to have key developmental functions in the free-living nematode Caenorhabditis elegans. Predicted mRNA targets of bpa-miR-5364 were identified using bioinformatics and comparative genomics approaches that relied on the conservation of miR-5364 binding sites in the orthologous mRNAs of other filarial nematodes. Finally, we confirmed the interaction between bpa-miR-5364 and three of its predicted targets using a dual luciferase assay.
These data provide new insight into the molecular mechanisms underpinning the transmission of third stage larvae of filarial nematodes from vector to mammal. This study is the first to identify parasitic nematode mRNAs that are verified targets of specific microRNAs and demonstrates that post-transcriptional control of gene expression via stage-specific expression of microRNAs may be important in the success of filarial infection.
- Lymphatic filariasis
The molecular mechanisms by which eukaryotic pathogens infect their hosts remain key questions in infection biology, as improved understanding of these events could lead to new methods of control. This is particularly so for vector-borne parasites, where the ability to survive and differentiate within the new environment of a mammalian host is critical for successful transmission. In this paper we explore the molecular events underlying the transmission of filarial nematodes from the mosquito to the mammalian host. We focus on microRNAs (miRNAs), which are small (~22 nucleotide) non-coding RNAs that regulate gene expression in animals, plants and viruses. miRNAs bind to complementary sequences, often in the 3′ untranslated region (UTR) of mRNAs, leading to translational repression and mRNA destabilisation . miRNAs are essential during development in the free-living nematode Caenorhabditis elegans, with functions in developmental timing [2-4], lifespan and stress responses [5,6], and embryogenesis [7,8]. C. elegans lin-4, the first miRNA to be identified , regulates developmental switching during the early larval stages (L1 and L2) by repressing the heterochronic proteins LIN-14  and LIN-28 , allowing epidermal cell fate transitions. C. elegans let-7 was the second miRNA to be identified and acts from the L3 stage onwards to regulate the transition from L4 to adult by down-regulating a number of targets, including LIN-41 [3,11]. let-7 is the founding member of a family of C. elegans miRNAs that share identity in the seed sequence (5′ nucleotides 2–7). This family includes a further six miRNAs; miR-48, miR-84, miR-241, miR-793, miR-794, and miR-795 . While let-7 is involved in specifying larval-adult cell fate, miR-48, miR-84 and miR-241 are co-activated at an earlier time point where they co-ordinate the L2-L3 transition .
Although information continues to accrue on miRNAs in the model nematode (reviewed recently in ), much less is known in parasitic nematodes. We recently identified 104 miRNAs from the filarial nematode Brugia pahangi , the sister species to the human pathogen B. malayi. In our study the miRNA identification process required that all sequences matched perfectly to the B. malayi genome and we therefore consider these miRNAs common to both species, with a recent study confirming this . B. malayi and the related Wuchereria bancrofti are the causative agents of lymphatic filariasis. Collectively, filarial worms infect around 120 million individuals in 73 countries worldwide [16,17]. These parasites cause a wide spectrum of pathology in the infected individual, including the debilitating conditions of elephantiasis and hydrocele. Lymphatic filariae have complex life cycles with no free-living forms; the parasites develop from first-stage larvae (microfilariae) to infective L3 within the mosquito intermediate host, and then from L3 to adults within the mammalian definitive host. Following transmission from an infected mosquito, the L3 migrate to the lymphatics where they develop through two moults to adults, which have a reproductive life-span of approximately 8 years . Mated females release millions of microfilariae into the bloodstream where they are available for ingestion by a mosquito taking a blood meal. Lymphatic filariae therefore have two periods of developmental arrest, as L3 in the mosquito and as microfilariae in the blood of the mammalian host. While parasite development is dependent upon transmission between hosts, little is known of the molecular mechanisms that control arrest and development within the different hosts.
In this paper we show that a novel member of the let-7 miRNA family, bpa-miR-5364, is highly up-regulated to coincide with the transition of the L3 from mosquito to mammalian host. The role of bpa-miR-5364 in the infection event was investigated by identifying potential mRNA targets using bioinformatic predictions, comparative genomics, and transcriptomic analysis, with selected targets verified using a mammalian cell transfection system. We propose that by regulating specific mRNAs, bpa-miR-5364 plays a key role in the transmission of the L3 stage from mosquito to mammalian host. This is the first example of a parasitic nematode miRNA for which function has been investigated and its mRNA targets identified and confirmed experimentally.
Expression profiling identifies developmentally regulated B. pahangi miRNAs
The complete lifecycle of B. pahangi is maintained in our laboratory allowing access to specific developmental stages. Therefore, from the miRNAs found in our previous genome-wide discovery study , we aimed to identify those with important developmental functions by detailed expression profiling at six key developmental time-points. The following life cycle stages were analysed: mosquito-derived L3, L3 isolated from the mammalian host at day 1 and day 5 post-infection (p.i.), L4 isolated at day 10 p.i., and sexually mature adult males and females isolated approximately 3 months p.i. The mosquito-derived L3 and day 1 p.i. L3 time-points were selected in order to identify changes occurring coincident with infection, day 5 p.i. L3 and day 10 p.i. L4 to identify changes prior to and after a moult, and males and females to define miRNA profiles specific to the adult stages. RNA was prepared from worms recovered at each time-point using six biological replicates, with three replicates used for microarray and three for qRT-PCR. A custom array (LC Sciences) contained probes against Brugia spp. [14,19] and C. elegans (miRBase release 15) miRNA sequences. This paper focuses on miRNAs differentially expressed in the infective stage as it transits from the mosquito to the mammalian host; our analysis of miRNAs functioning specifically in adult parasites will be presented elsewhere. The microarray data for the larval stages is provided in Additional file 1.
Significant changes in miRNA profiles are linked to invasion of the host
Mean bpa-miR-5364 microarray signal in five B. pahangi life cycle stages
Reporter name (a)
In addition to bpa-let-7 and bpa-miR-84, the three other miRNAs that met the criteria were bpa-miR-5364, cel-miR-789 and bpa-miR-5853. We chose to focus on bpa-miR-5364 for the following reasons: its expression profile is particularly striking as it is up-regulated ~12 fold (Log2 3.64) within 24 hours of infection of the mammalian host and doubles again over the period to day 10 p.i. (Table 1); bpa-miR-5364 shares sequence similarity with the let-7 miRNA family and appears to be specific to clade III parasitic nematodes (as described below); it is one of the most abundant miRNAs with only bpa-lin-4 and bpa-miR-71 found at higher levels at day 10 p.i. (Additional file 1). Abundant miRNA are often evolutionarily ancient  suggesting important function.
bpa-miR-5364 is a novel member of the let-7 miRNA family
Expression of bpa-miR-5364 is induced in culture
To examine potential mechanisms regulating expression, the genomic sequences 5′ and 3′ of bpa-mir-5364 were analysed. Brugia spp. mir-5364 is an intergenic miRNA with the closest protein coding region residing 1.2 kb upstream and is therefore likely to be regulated by its own promoter. However, analysis of this region revealed no significant transcription factor binding sites (not shown). In C. elegans, let-7 is autoregulated via a let-7 binding site situated at the 3′ end of the primary transcript approximately 500 bp downstream of the sequence encoding the mature miRNA . Autoregulation of C. elegans lin-4 also occurs, in this case by activation of lin-4 transcription via a lin-4 binding site in its promoter . The B. malayi mir-5364 genomic sequence was analysed using the Probability of Interaction by Target Accessibility (PITA) algorithm . Interestingly, predicted binding sites for mir-5364 were found 614 bases downstream (7 nucleotide seed containing one G:U pairing, ΔΔG of −6.15) and 798 bases upstream (8 nucleotide seed containing one G:U pairing, ΔΔG −11.21) of the mir-5364 mature region. This suggests that similar mechanisms of autoregulation to those described for let-7 and lin-4 in C. elegans might exist for mir-5364 in the parasite.
mir-5364 is specific to clade III parasitic nematodes
miRBase (release 21) (http://www.mirbase.org/) shows that, to date, mir-5364 has been identified only in the parasitic nematodes B. pahangi  and Ascaris suum , both of which are members of nematode clade III . mir-5364 is not present in the deep-sequencing-defined miRNA sets from the animal parasitic nematodes Haemonchus contortus (clade V)  or Strongyloides ratti (clade IVa) , nor in the free-living nematode species C. elegans, C. briggsae, C. brenneri, C. remanei, Pristionchus pacificus (all clade V) or Panagrellus redivivus (clade IV) [26-29], suggesting that mir-5364 miRNA is specific to clade III parasitic nematodes. To investigate this further, we examined a number of other nematode genomes to determine the extent of mir-5364 conservation. Using the bpa-mir-5364 precursor sequence as a query for BLASTN analysis at 959 Nematode Genomes (http://xyala.cap.ed.ac.uk/downloads/959nematodegenomes/blast/blast.php, with default settings), mir-5364 was found in all clade III parasitic nematode genomes for which sequence was available. Regions of ≥ 84% identity were identified in the filarial nematodes Dirofilaria immitis, Loa loa, Wuchereria bancrofti, Litomosoides sigmodontis and Onchocerca ochengi (identity ranging from 66/78 nucleotides in O. ochengi to 128/137 nucleotides in W. bancrofti) with complete conservation of the mature miRNA sequence in all five species (Additional file 3). When folded using mfold  (http://mfold.rit.albany.edu/?q=mfold/RNA-Folding-Form, default settings) the regions of precursor homology identified from each species formed miRNA-like hairpins (Additional file 3) which also passed the miRNA hairpin classifier algorithm miPred (http://www.bioinf.seu.edu.cn/miRNA/index.html)  (Additional file 4). In contrast, our bioinformatic analyses showed no evidence for mir-5364 orthologs in the genomic data at 959 Nematode Genomes of the non-clade III nematodes Bursaphelenchus xylophilus, Heligosomoides polygyrus, Howardula aoronymphium or Oscheius tipulae using either the precursor or mature sequences as BLAST queries. Although it is possible that divergent sequences might not have been detected using BLASTN, miR-5364 was not found in H. polygyrus by deep sequencing of small RNAs . Therefore, as miR-5364 was also recently experimentally confirmed in B. malayi , this miRNA is highly conserved in eight clade III parasitic nematodes but has not been identified in any other species to date.
Computational and phylogenetic analysis predicts target mRNAs
As miRNAs bind to partially complementary sequences in their target mRNAs, computational analysis can be used to predict targets [33,34]. Genome-wide target predictions for miR-5364 were performed utilising B. malayi and D. immitis genomic data [35,36]. Although miRNA target sites have been defined throughout genes , we focussed on the 3′UTR, as many well-characterised target sites reside in this region [1,9,11,38,39]. Also, rather than rely on annotated 3′UTRs, which are often dependent on developmental stage  or cell type, and for which the current B. malayi and D. immitis data are incomplete, we used 1 kb of genomic sequence downstream of the annotated end of all genes. The miRNA target prediction algorithms PITA  and miRanda  were used and only predictions found using both programs were retained. The genome-wide miR-5364 target predictions for both parasites are provided in Additional file 5.
Predicted mRNA targets of bpa-miR-5364
B. malayi g ene annotation
B. malayi ortholog
D. immitis ortholog
C. elegans ortholog (a)
Delta-1-pyrroline-5-carboxylate dehydrogenase, mitochondrial precursor, putative
Zinc finger in N-recognin family protein
Twik (KCNK-like) family of potassium channels, alpha subunit 5
gag protein, putative
ADP,ATP carrier protein, heart/skeletal muscle isoform T1, putative
high mobility group protein, putative
Ets-domain containing protein
LEM domain containing protein
NifU-like protein, putative
protein C33A11.1, putative
Zinc finger, C2H2 type family protein
The expression profile of predicted target mRNAs correlates inversely with that of bpa-mir-5364
bpa-miR-5364 interacts directly with predicted mRNA targets
Attempts at bpa-miR-5364 inhibition using an antisense oligonucleotide
We also tested the effect of RNAi of the single predicted miRNA-specific argonaute gene Bm1_13960  using endoribonuclease-prepared pools of complex siRNA (esiRNAs) [49-51], together with the bpa-miR-5364 ASO. This approach is similar to that used in C. elegans where miRNA mutants were analysed in the background of global reduction in miRNA levels due to alg-1 mutation . This combined approach did not produce any phenotype or de-repression of target mRNAs in B. pahangi (not shown).
Co-ordinated regulation of mRNA targets by multiple miRNAs
As individual mRNAs are often regulated by multiple miRNAs, this most likely explains the lack of de-repression with the miR-5364 ASO. Therefore predictions of additional miRNAs that might function alongside bpa-miR-5364 to target Bm1_27305, Bm1_05425 and Bm1_25620, and their D. immitis orthologs, were made using PITA and a set of miRNA which included all 104 mature Brugia spp. miRNA identified previously . Predicted miRNA target sites for the three target genes in both parasite species are shown in Additional file 9. Interestingly, conserved sites for let-7 and miR-84 are found for Bm1_27305 with almost identical positioning in each species. Given the sequence similarity between miR-5364 and the let-7 family (Figure 2A) this might be expected; however, the binding sites for miR-5364 and let-7/miR-84 are not identically positioned (see Additional file 9). In addition, the C. elegans ortholog of Bm1_27305, C33A11.4 (Table 2) also known as tag-97, has predicted binding sites for let-7 family miRNAs using PITA (results not shown) and TargetScan  (http://www.targetscan.org/worm_52/), indicating that regulation of Bm1_27305 and its orthologs by this important miRNA family may be a feature conserved throughout the nematode phylum. Similarly, using PITA the C. elegans ortholog of Bm1_05425 (T22C1.1) has a site for the let-7 family member miR-794 (results not shown). Predicted binding sites for the orthologous miRNAs bantam and mir-82  are found for Bm1_25620 (Additional file 9) and in the TargetScan predictions for the C. elegans ortholog (Y48B6A.14, hmg-1.1), suggesting conservation of miRNA-mediated regulation across species. Further work will reveal if these interactions are functional and if these miRNAs act in conjunction with miR-5364 to regulate its targets.
Infection of the mammalian host by the L3 stage of filarial worms is a key event in the life cycle: how the parasite survives and differentiates within its new environment is poorly understood. Here we show that bpa-miR-5364, which is specific to clade III parasitic nematodes and related to the let-7 family, is highly up-regulated to coincide with transmission to the mammalian host and may be involved in the infection event by targeting specific mRNAs. The L3 of many parasitic nematodes show parallels with the developmentally arrested dauer larva of C. elegans. In a recent study, four C. elegans miRNAs were found to be up-regulated in samples of mixed stages compared to dauers, two of which are the let-7 family members miR-241 and miR-795 . Similarly, in a detailed analysis of miRNA expression following dauer or continuous development, let-7 family members showed increased expression post-dauer . Brugia spp. possess three let-7 family members, let-7, miR-84 and miR-5364, all of which are up-regulated post-infection, with expression of miR-5364 increasing first, within 24 hours, and rising to a significantly greater level than either let-7 or miR-84, suggesting that in Brugia spp. miR-5364 has a dominant role. We therefore aimed to identify the possible functions of this miRNA by identifying its target genes.
Identifying mRNA targets is complex because most bioinformatics tools search primarily for complementarity between the short miRNA seed sequence (5′ nucleotides 2–7) and the mRNA, and then take other factors into account, including the free energy of the miRNA-mRNA duplex and position of binding within the 3′UTR. Thus, bioinformatics tools potentially identify hundreds of possible targets, many of which may not be verifiable. In an attempt to overcome this problem, the availability of other filarial genomes was exploited in a comparative genomics approach. We identified orthologous mRNAs in different clade III nematodes that have conserved binding sites for miR-5364, with the further requirement that the site be similarly positioned. This approach, which has been used previously in vertebrates and nematodes [53,55], significantly reduced the number of possible target mRNAs from over 850 to 13. Three of these were selected for in-depth study: Bm1_27305 encoding an ETS domain transcription factor, Bm1_05425 encoding an N-recognin-type zinc finger protein, and Bm1_25620 which encodes a putative high mobility group protein. All three genes were found to have a reciprocal expression profile to bpa-miR-5364, both from analysis of transcriptomic data and by qRT-PCR, consistent with miRNA-mediated negative regulation. To demonstrate a direct interaction between miRNA and target we used transient transfection of mammalian HEK293 cells and a dual luciferase assay system and found that the presence of bpa-miR-5364 resulted in the down-regulation of the luciferase reporter gene for all three 3′UTRs, providing further evidence that these genes are true targets of this miRNA.
All three verified targets of bpa-miR-5364 could potentially function to alter gene and/or protein levels and a decrease in the levels of these post-infection may be important for larval activation and/or developmental progression in the mammalian host. The human ortholog of Bm1_27305 (EHF, NCBI Gene ID: 26298) belongs to the ETS domain transcription factor family, which are characterised by epithelial-specific expression (reviewed in ) and by tight tissue-specific and temporal regulation. For example, in prostate cancer this transcription factor controls the balance between cell proliferation and differentiation . The C. elegans homolog of Bm1_27305 (C33A11.4) is upregulated in mutants of the miRNA-specific argonaute gene alg-1  and is enriched in AIN-1 (ALG-1 INteracting protein) pull-downs , consistent with miRNA regulation. Bm1_05425 contains an N-recognin zinc finger motif conserved in ubiquitin-ligases, which target proteins for degradation. In C. elegans, an important function of let-7 is the negative regulation of LIN-41 [3,11], which contains a ubiquitin ligase domain. In turn, LIN-41 negatively regulates levels of the transcription factor LIN-29, which is required for adult cell fates [3,11]. It is possible that in Brugia spp. L3, Bm1_05425 may perform a similar function to LIN-41 in C. elegans, as proteins with ubiquitin ligase activity could potentially target a subset of proteins for degradation. The final mir-5364 target studied, Bm1_25620, is predicted to encode an HMG protein. These bind to chromatin and can act in the regulation of transcription during development, DNA repair, apoptosis or cell senescence . In mammalian cells, HMGA2 is abundant during early development and levels decrease during later stages of development and differentiation, via let-7 regulation . A reduction in HMG Bm1_25620 following transmission to the mammalian host may allow larval development and/or exit from an arrested state.
It is important to note that whatever the function of bpa-miR-5364, it is not likely to be restricted to the L3 stage, as in the life cycle stages tested, expression peaked in day 10 p.i. L4, an actively growing stage and continued into adulthood. Whether the target mRNAs are the same between the p.i. L3 and subsequent developmental stages of Brugia spp. is currently unknown. Even in C. elegans where miRNAs are comparatively well defined, little is known of target conservation between different life cycle stages.
In an attempt to define the possible function of the mRNA targets in B. pahangi, we tested the effects of inhibiting bpa-miR-5364 in the larval stages using a specific ASO. ASOs bind with sequence complementarity to the mature miRNA to disrupt binding to mRNA targets. While we found efficient uptake of labelled RNA by B. pahangi L3 and p.i. L3s in the culture conditions tested, no de-repression of any of the three verified target genes was found using the ASO. As loss of a single miRNA is often not sufficient to reveal function, we also tested the effect of the mir-5364 ASO along with siRNAs targeting the Brugia spp. miRNA-specific argonaute, but found no observable effect. Therefore, it may be necessary to include specific inhibitors of additional miRNAs to reveal the effect of bpa-miR-5364 disruption. However, the conditions required for inhibition of miRNA function in parasitic nematodes are completely unknown, with this being the first study to attempt to inhibit miRNAs. Factors such as the concentration of ASO required, efficiency of uptake in different cells and tissues, and the ability of exogenously applied inhibitors to enter the RISC and block target binding will all influence the ability to inhibit miRNA function.
As expression of bpa-miR-5364 increased 12-fold between vector derived L3 and 24 h p.i.L3, it was of interest to determine what factors might induce the expression of this miRNA. Previous studies on the control of gene expression in vector-derived L3 versus 24 hour p.i. L3 highlighted a number of mRNAs which appeared to be regulated largely by the temperature shift associated with infection . By simply moving L3 to 37°C in an atmosphere of 5% CO2 in air, the expression of bpa-miR-5364 was induced, with serum having minimal additional effect. Recent studies in Arabidopsis have shown that increased temperature and CO2 lead to significant changes in expression of several miRNAs , but whether temperature or CO2 has any role in the induction of bpa-mir-5364 requires further study. In addition, autoregulation of let-7 has been identified in C. elegans, in which binding of let-7 to its own 3′UTR acts to increase processing of the pri-miRNA by argonaute, leading to higher levels of the mature miRNA . Similarly, binding of C. elegans lin-4 to a site in its own promoter activates lin-4 transcription . In this study, we identified potential bpa-miR-5364 binding sites in genomic sequence upstream and downstream of the mature miRNA region, suggesting that autoregulation could also be important for control of this parasite miRNA.
mir-5364 appears to be restricted to clade III nematodes as, to date, related sequences have not been found in other parasitic or free-living nematode genomes. In addition to mir-5364 in Brugia spp. [14,15], we found this miRNA in the genomes of five additional clade III filarial parasites, with a recent study confirming expression in D. immitis . mir-5364 is not restricted to filarial nematodes as it is also expressed in the clade III gastrointestinal nematode A. suum  suggesting its role is not specific to tissue-dwelling nematodes. Therefore mir-5364 may have diverged from the let-7 family during the evolution of clade III nematodes, an event that occurred 350–540 million years ago . These observations support the concept that small RNA sequence and function have evolved and diversified between nematode clades, as suggested previously following the identification of 21U (or piwi-interacting) RNAs in clade V nematodes, but not, to date, in clade III [14,24] or clade IV  nematodes.
In conclusion, we report a B. pahangi miRNA that is highly up-regulated to coincide with the transition between hosts and identify some of its mRNA targets. bpa-miR-5364 is a novel member of the let-7 family, consistent with the function of other members of that family in nematode development. Infection of the mammalian host is the key event in the filarial life cycle and our data indicate that post-transcriptional regulation of gene expression may be necessary for this event.
B. pahangi maintenance and culture
B. pahangi was maintained as described previously [14,66]. For developmental gene expression analysis, mosquito-derived L3 were frozen immediately following isolation from the mosquito, while p.i. larval stages and adult worms were collected by peritoneal lavage of infected jirds at specific time-points. B. pahangi larvae for culture were prepared as described previously  and cultured in Worm Culture Medium (WCM) (RMPI 1640 with glutamine and HEPES, supplemented with 100 units/ml of penicillin and 100 μg/ml streptomycin, 5% heat inactivated fetal calf serum (FCS), and 1% glucose, all from Invitrogen) at 37°C in 5% CO2/air. All animal protocols were carried out in accordance with the guidelines of the UK Home Office, under the Animal (Scientific Procedures) Act 1986, following approval by the University of Glasgow Ethical Review Panel. Experiments were performed under the authority of the UK Home Office, project numbers 60/4448 and 60/3792.
Total RNA extraction and miRNA microarray analysis
Total RNA was extracted using the Trizol protocol (Invitrogen) after manual disruption, with RNA for miRNA microarray prepared with an additional 75% ethanol wash. Biological replicates were defined as parasites isolated from different animals, or, for mosquito-derived L3, extracted at different times. Microarrays were performed by a service provider (LC Sciences) using a custom array containing probes to 252 Brugia spp. sequences (all the mature and star miRNA sequences identified in our previous study  with three additional sequences found by others ), and all C. elegans sequences present in miRBase release 15 (ftp://mirbase.org/pub/mirbase/15/).
miRNA and mRNA qRT-PCR
For miRNA qRT-PCR, the miRNA 1st Strand cDNA Synthesis protocol (Agilent Technologies) was followed using 1 μg of DNase I-treated (Ambion) total RNA for the polyadenylation reaction and enzyme-plus and enzyme-minus reverse transcription reactions. qPCR was performed following the miRNA qPCR Master Mix protocol (Agilent Technologies). Results were normalised to bpa-miR-100c, which was found by microarray to be present at similar levels throughout the developmental stages examined. For mRNA qRT-PCR, 1 μg of DNase I-treated total RNA was used in enzyme-plus and enzyme-minus reverse transcriptions using the SuperScript II (Invitrogen) protocol with oligo(dT) primer. mRNA qRT-PCR was performed following the Brilliant III Ultra-Fast SYBR Green QPCR Master Mix (Agilent Technologies) protocol using β-tubulin (Genbank AY705382) as a normalising gene. For both miRNA and mRNA qRT-PCR, as well as biological triplicate samples, each PCR was carried out in triplicate and run on an Agilent Mx3005P qPCR System. The sequences of the HPLC-purified oligonucleotides used are given in Additional file 10.
Induction of bpa-miR-5364 in cultured parasites
Four sets of 500 L3 were isolated from mosquitoes and prepared for culture as described above. One set was frozen immediately while the other three sets were cultured at 37°C, 5% CO2 in air in WCM adjusted to contain either no FCS, 5% FCS, or 10% FCS. After 24 hours the larvae were collected and frozen. Total RNA was extracted as described above and the entire preparation DNAase I-treated (Ambion), ethanol precipitated and used for qRT-PCR.
Bioinformatic prediction of miR-5364 targets
The B. malayi genome  was downloaded from the following URL, http://www.ncbi.nlm.nih.gov/nuccore?term=DS236884:DS264093[PACC] as genomic scaffolds. From the same site the annotated genome was downloaded in Genbank(full) format then converted to GFF using ReadSeq (http://iubio.bio.indiana.edu/soft/molbio/readseq/java/). B. malayi proteins sequences were downloaded from http://ghedinlab.csb.pitt.edu/GhedinLab/Parasite%20Genomics?action=AttachFile&do=view&target=brugia_aa_112707.fa.gz. The Dirofilaria immitis genome  version 2 reassembly files “Dirofilaria_immitis_2.2”, “nDi.2.2.2.aug.blast2go”, and “nDi.2.2.2.aug.proteins” were downloaded from http://badger.bio.ed.ac.uk/filarial/home/download. Use of this updated version of the published D. immitis genome was kindly provided pre-publication by Sujai Kumar and colleagues. Sequences 1 kb 3′ of the annotated stop codon for all genes were extracted for both species. The miRNA target prediction algorithms PITA  (http://genie.weizmann.ac.il/pubs/mir07/mir07_exe.html) and miRanda , (http://www.microrna.org/microrna/getDownloads.do) were run as standalone programs (using default settings) on the 3′ downstream sequence sets from each species searching for mir-5364 sites. PITA default settings permit seeds of the following types, 6:0:0, 7:0:1, 8:1:1 (seed length: number of G-U pairings allowed in seed: number of mismatches allowed in seed). For each species the PITA and miRanda predictions were combined and only target predictions found by both methods were retained. The predicted protein sets for B. malayi and D. immitis were analysed using InParanoid  (http://inparanoid.sbc.su.se/cgi-bin/faq.cgi) to identify orthologs. Reciprocal BLASTP analysis at NCBI (http://www.blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome) was used to identify the orthologs of B. malayi genes in L. loa and W. bancrofti, while for L. sigmodontis one-way TBLASTN was performed at 959 Nematode Genomes (http://xyala.cap.ed.ac.uk/downloads/959nematodegenomes/blast/blast.php). To test for the presence of mir-5364 binding sites, sequence 1 kb downstream from the annotated stop of each gene was analysed using PITA online (http://genie.weizmann.ac.il/pubs/mir07/mir07_prediction.html, with settings: 7 nucleotide seed, no mismatch allowed, 1 G:U allowed). Reciprocal BLAST analysis was used to identify the C. elegans orthologs.
Prediction of additional miRNAs targeting Bm1_27305, Bm1_05425 and Bm1_25620
The downstream sequences of Bm1_27305, Bm1_05425 and Bm1_25620 and their D. immitis orthologs were analysed for the presence of miRNA binding sites using the set of 252 Brugia spp. miRNA sequences (described above) with PITA online. Settings for Brugia spp. were as follows: 7 nucleotide seed, 0 mismatch allowed, 1 G-U allowed; and for D. immitis 7:1:1 to compensate for potential miRNA sequence variation. Only those sites with PITA ΔΔG energies of < −7 and located within the longest experimentally defined 3′UTR (see below) were considered.
3′RACE analysis of predicted mRNA targets of bpa-miR-5364
mRNA 3′ ends were analysed using the 3′ RACE System for Rapid Amplification of cDNA Ends (Invitrogen) following the manufacturer’s protocol. 1 μg of DNase I (Ambion) treated mosquito-derived B. pahangi L3 total RNA was used along with the gene-specific primers designed to B. malayi sequences (primers in Additional file 10). 3’ RACE products were cloned into pCR2.1-TOPO (Invitrogen) and sequenced. The longest experimentally defined 3′UTRs were 418, 182 and 504 nucleotides downstream of the annotated stop for Bm1_27305, Bm1_05425, and Bm1_25620 respectively.
RNAseq library construction and sequencing
B. pahangi RNA isolated from mosquito-derived L3, day 1, day 5 and day 10 p.i. larvae was used to generate standard Illumina RNAseq libraries with two biological replicates. Libraries were sequenced using the Illumina HiSeq 2000 platform as 100 bp pair-end sequences. Sequencing reads were submitted to the European Nucleotide Archive (experiment accession number: ERX199547; http://www.ebi.ac.uk/ena/data/view/ERX199547).
Mapping of sequencing reads and differential expression analysis
Mapping of sequencing reads was performed using Tophat version v2.0.8b  with default parameters except for “-g 1 --library-type fr-unstranded -r 200 --mate-std-dev 100 -a 6 -i 10 -I 40000 --microexon-search --min-segment-intron 10 --max-segment-intron 40000”. The B. pahangi draft genome version 1.5.4, available through Wormbase ParaSite, was used as reference (http://parasite.wormbase.org/Brugia_pahangi_prjeb497/). Alignment files resulting from the mapping step were further sorted and indexed using SAMtools version: 0.1.19-44428 cd  using the sort and index functions respectively. Counts per transcript were retrieved using the htseq-count module from the HTSeq package version 0.6.0  with default parameters except for “-f bam -a 30 -t exon -s no -m union”. The reference annotation used in this analysis is the version WBPS1 available through Wormbase ParaSite. Differential expression analysis was performed using the R language package EdgeR version 3.4.2 . Normalised counts plots were generated using ggplot2 version 1.0.0.
Analysis of miRNA/mRNA interactions in HEK293 cells
3′UTRs from Bm1_27305, Bm1_05425 and Bm1_25620 were amplified from B. pahangi genomic DNA (isolated as described previously ) using PfuUltra II Fusion HS DNA Polymerase (Agilent) (primers in Additional file 10). Products were cloned into pCR2.1-TOPO (Invitrogen) then subcloned into the Not I site downstream of firefly luciferase in pMir-Target (Origene) to generate plasmids pABR212 (Bm1_27305 3′UTR in pMir-Target), pABR210 (Bm1_05425 3′UTR in pMir-Target) and pABR211 (Bm1_25620 3′UTR in pMir-Target). A construct containing five repeats of the perfect target sequence for bpa-miR-5364 was synthesised by GeneArt (Invitrogen) (sequence in Additional file 10) and cloned into pMir-Target using Eco RI and Not I to generate plasmid pMir-T:Bp-mir-5364-5Xtarget. A 317 bp section of the bpa-mir-5364 locus was amplified from B. pahangi genomic DNA (primers in Additional file 10), cloned into pCR2.1-TOPO and subcloned into vector pEGFP-N1 (Clontech) using Kpn I to generate plasmids pABR213 (bpa-mir-5364 in vector pEGFP-N1, sense orientation) and pEGFPN1:Bp-mir-5364-Rev (bpa-mir-5364 in vector pEGFP-N1, reverse orientation).
HEK293 cells were maintained in Dulbecco's Modified Eagle's Medium (containing 4500 mg/L glucose and sodium bicarbonate, Sigma D5671) (supplemented with 2 mM L-Glutamine, 100 units/ml of penicillin, 100 μg/ml streptomycin, and 10% heat inactivated FCS, all from Invitrogen) at 37°C in 5% CO2/air. For transfections, 2 × 104 cells/well were seeded into the wells of 96-well plates in a volume of 100 μl using the above growth medium but without addition of antibiotics. Cells were transfected after 24 hours when they were ~60% confluent using Lipofectamine LTX (Invitrogen) with 50 ng of a mir-5364-containing plasmid, 6.25 ng of a pMir-Target-derived plasmid and 0.5 ng of phRG-TK (Renilla luciferase, Promega). Transfections were performed following the manufacturer’s protocol for HEK293 cells. DNA was diluted to a volume of 20 μL using Opti-MEM I Reduced Serum Medium (Invitrogen) and 0.35 μL Lipofectamine LTX Reagent added. After incubation at room temperature for 30 minutes 20 μL of the DNA-Lipofectamine complex was then added directly to the cells in each well. Cells were grown for 48 hours then analysed using a Dual Luciferase Assay kit (Promega) following the manufactures protocol with six replicates used per test condition.
Conditions for uptake of double-stranded and single-stranded RNA from medium by B. pahangi were tested using Silencer Cy3 Labeled Negative Control No. 1 siRNA (Ambion) and Cy3 dye-Labeled Anti-miR Negative Control No. 1 (Ambion) at a final concentration of 0.5 μM. ~100 mosquito-derived L3 and ~30 day 5 p.i. L3, both prepared as described above, were added to 0.5 ml WCM (5% FCS) in a 48-well culture plate, with and without the addition of Cy3-siRNA and Cy3-anti-miR, and uptake analysed after 24 and 48 hours. An antisense oligonucleotide (ASO)  to inhibit bpa-mir-5364 was purchased from Eurofins MWG Operon along with a scrambled version of the ASO (sequences in Additional file 10). The ASOs consisted of all 2’O-methyl RNA bases, phosphorothioate bonds at the 5′ and 3′ ends, as well as nucleotides added at each end. The B. pahangi ortholog of the full-length predicted Bm1_13960 (argonaute) sequence was amplified from mosquito-derived B. pahangi cDNA using PfuUltra II Fusion HS DNA Polymerase (Agilent) (primers in Additional file 10). A ~1.7 kb Pst I/Pst I fragment of this was cloned into pPD129.36 (Addgene) to create plasmid pADW0062 which was linearised with Not I and Hind III. In vitro transcription and preparation of endoribonuclease-prepared pools of complex siRNA (esiRNA) was carried out as described previously . Larvae were harvested at either 4 or 5 days p.i. and cultured in 24-well plates containing 1 ml WCM with ~150-200 larvae per well. Larvae were treated under three conditions, 1) 5 μM mir-5364 ASO plus 1 μM Bm1_13960 (argonaute) esiRNA, 2) 5 μM scrambled ASO plus 1 μM Bm1_13960 esiRNA, and 3) medium only. Each condition was replicated in three wells and the experiment repeated twice. Larvae were cultured for 48 hours, frozen at −80°C and analysed by qRT-PCR.
B. pahangi larvae were prepared for microscopy by washing three times in 1 ml of PBS/0.1% Tween in 1.5 ml tubes, using centrifugation at 2000 rpm in a benchtop microfuge for 2 minutes. After the final wash, the volume was reduced to ~20 μL and 10 μL pipetted onto microscope slides with a 2% agarose pad. Images were collected using an Axioskop 2 Plus microscope (Zeiss), ORCA-ER digital camera (Hamamatsu) and Openlab (Improvision) software.
The authors thank Kerry O’Neil and Margaret McFadyen for maintaining the B. pahangi lifecycle at the University of Glasgow, and Gillian Borland (University of Glasgow) for advice on cell transfections. This work was funded by a project grant awarded to ED and CB by the Wellcome Trust (WT 086823/Z/08/Z). AP, NEH and MB are funded by the Wellcome Trust (grant 098051), BR is funded by EC FP7 PARAVAC consortium, and RDE is funded by the University of Nottingham.
- Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP. The impact of microRNAs on protein output. Nature. 2008;455:64–71.View ArticlePubMed CentralPubMedGoogle Scholar
- Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993;75:843–54.View ArticlePubMedGoogle Scholar
- Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, et al. The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000;403:901–6.View ArticlePubMedGoogle Scholar
- Abbott AL, Alvarez-Saavedra E, Miska EA, Lau NC, Bartel DP, Horvitz HR, et al. The let-7 microRNA family members mir-48, mir-84, and mir-241 function together to regulate developmental timing in Caenorhabditis elegans. Dev Cell. 2005;9:403–14.View ArticlePubMed CentralPubMedGoogle Scholar
- de Lencastre A, Pincus Z, Zhou K, Kato M, Lee SS, Slack FJ. MicroRNAs both promote and antagonize longevity in C. elegans. Curr Biol. 2010;20:2159–68.View ArticlePubMed CentralPubMedGoogle Scholar
- Boulias K, Horvitz HR. The C. elegans microRNA mir-71 acts in neurons to promote germline-mediated longevity through regulation of DAF-16/FOXO. Cell Metab. 2012;15:439–50.View ArticlePubMed CentralPubMedGoogle Scholar
- Alvarez-Saavedra E, Horvitz HR. Many families of C. elegans microRNAs are not essential for development or viability. Curr Biol. 2010;20:367–73.View ArticlePubMed CentralPubMedGoogle Scholar
- Shaw WR, Armisen J, Lehrbach NJ, Miska EA. The conserved miR-51 microRNA family is redundantly required for embryonic development and pharynx attachment in Caenorhabditis elegans. Genetics. 2010;185:897–905.View ArticlePubMed CentralPubMedGoogle Scholar
- Wightman B, Ha I, Ruvkun G. Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans. Cell. 1993;75:855–62.View ArticlePubMedGoogle Scholar
- Moss EG, Lee RC, Ambros V. The cold shock domain protein LIN-28 controls developmental timing in C. elegans and is regulated by the lin-4 RNA. Cell. 1997;88:637–46.View ArticlePubMedGoogle Scholar
- Slack FJ, Basson M, Liu ZC, Ambros V, Horvitz HR, Ruvkun G. The lin-41 RBCC gene acts in the C. elegans heterochronic pathway between the let-7 regulatory RNA and the LIN-29 transcription factor. Mol Cell. 2000;5:659–69.View ArticlePubMedGoogle Scholar
- Roush S, Slack FJ. The let-7 family of microRNAs. Trends Cell Biol. 2008;18:505–16.View ArticlePubMedGoogle Scholar
- Grishok A. Biology and Mechanisms of Short RNAs in Caenorhabditis elegans. Adv Genet. 2013;83:1–69.View ArticlePubMedGoogle Scholar
- Winter AD, Weir W, Hunt M, Berriman M, Gilleard JS, Devaney E, et al. Diversity in parasitic nematode genomes: the microRNAs of Brugia pahangi and Haemonchus contortus are largely novel. BMC Genomics. 2012;13:4.View ArticlePubMed CentralPubMedGoogle Scholar
- Poole CB, Gu W, Kumar S, Jin J, Davis PJ, Bauche D, et al. Diversity and expression of microRNAs in the filarial parasite. Brugia malayi. PLoS One. 2014;9:e96498.View ArticlePubMed CentralPubMedGoogle Scholar
- Ottesen EA, Duke BOL, Karam M, Behbehani K. Strategies and tools for the control/elimination of lymphatic filariasis. Bull World Health Organ. 1997;75:491–503.PubMed CentralPubMedGoogle Scholar
- Hotez PJ, Brindley PJ, Bethony JM, King CH, Pearce EJ, Jacobson J. Helminth infections: the great neglected tropical diseases. J Clin Invest. 2008;118:1311–21.View ArticlePubMed CentralPubMedGoogle Scholar
- Subramanian S, Stolk WA, Ramaiah KD, Plaisier AP, Krishnamoorthy K, Van Oortmarssen GJ, et al. The dynamics of Wuchereria bancrofti infection: a model-based analysis of longitudinal data from Pondicherry, India. Parasitology. 2004;128:467–82.View ArticlePubMedGoogle Scholar
- Poole CB, Davis PJ, Jin JM, McReynolds LA. Cloning and bioinformatic identification of small RNAs in the filarial nematode, Brugia malayi. Mol Biochem Parasitol. 2010;169:87–94.View ArticlePubMedGoogle Scholar
- Zhu Y, Skogerbo G, Ning Q, Wang Z, Li B, Yang S, et al. Evolutionary relationships between miRNA genes and their activity. BMC Genomics. 2012;13:718.View ArticlePubMed CentralPubMedGoogle Scholar
- Zisoulis DG, Kai ZS, Chang RK, Pasquinelli AE. Autoregulation of microRNA biogenesis by let-7 and Argonaute. Nature. 2012;486:541–4.PubMed CentralPubMedGoogle Scholar
- Turner MJ, Jiao AL, Slack FJ. Autoregulation of lin-4 microRNA transcription by RNA activation (RNAa) in C. elegans. Cell Cycle. 2014;13:772–81.View ArticlePubMed CentralPubMedGoogle Scholar
- Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. The role of site accessibility in microRNA target recognition. Nat Genet. 2007;39:1278–84.View ArticlePubMedGoogle Scholar
- Wang J, Czech B, Crunk A, Wallace A, Mitreva M, Hannon GJ, et al. Deep small RNA sequencing from the nematode Ascaris reveals conservation, functional diversification, and novel developmental profiles. Genome Res. 2011;21:1462–77.View ArticlePubMed CentralPubMedGoogle Scholar
- Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, Vierstraete A, et al. A molecular evolutionary framework for the phylum Nematoda. Nature. 1998;392:71–5.View ArticlePubMedGoogle Scholar
- Ahmed R, Chang Z, Younis AE, Langnick C, Li N, Chen W, et al. Conserved miRNAs are candidate post-transcriptional regulators of developmental arrest in free-living and parasitic nematodes. Genome Biol Evol. 2013;5:1246–60.View ArticlePubMed CentralPubMedGoogle Scholar
- de Wit E, Linsen SEV, Cuppen E, Berezikov E. Repertoire and evolution of miRNA genes in four divergent nematode species. Genome Res. 2009;19:2064–74.View ArticlePubMed CentralPubMedGoogle Scholar
- Shi Z, Montgomery TA, Qi Y, Ruvkun G. High-throughput sequencing reveals extraordinary fluidity of miRNA, piRNA, and siRNA pathways in nematodes. Genome Res. 2013;23:497–508.View ArticlePubMed CentralPubMedGoogle Scholar
- Srinivasan J, Dillman AR, Macchietto MG, Heikkinen L, Lakso M, Fracchia KM, et al. The draft genome and transcriptome of Panagrellus redivivus are shaped by the harsh demands of a free-living lifestyle. Genetics. 2013;193:1279–95.View ArticlePubMed CentralPubMedGoogle Scholar
- Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31:3406–15.View ArticlePubMed CentralPubMedGoogle Scholar
- Jiang P, Wu H, Wang W, Ma W, Sun X, Lu Z. MiPred: classification of real and pseudo microRNA precursors using random forest prediction model with combined features. Nucleic Acids Res. 2007;35:W339–44.View ArticlePubMed CentralPubMedGoogle Scholar
- Buck AH, Coakley G, Simbari F, McSorley HJ, Quintana JF, Le Bihan T, et al. Exosomes secreted by nematode parasites transfer small RNAs to mammalian cells and modulate innate immunity. Nat Commun. 2014;5:5488.View ArticlePubMed CentralPubMedGoogle Scholar
- Rajewsky N. microRNA target predictions in animals. Nat Genet. 2006;38:S8–13.View ArticlePubMedGoogle Scholar
- Ritchie W, Rasko JE. Refining microRNA target predictions: sorting the wheat from the chaff. Biochem Biophys Res Commun. 2014;445:780–4.View ArticlePubMedGoogle Scholar
- Ghedin E, Wang SL, Spiro D, Caler E, Zhao Q, Crabtree J, et al. Draft genome of the filarial nematode parasite Brugia malayi. Science. 2007;317:1756–60.View ArticlePubMed CentralPubMedGoogle Scholar
- Godel C, Kumar S, Koutsovoulos G, Ludin P, Nilsson D, Comandatore F, et al. The genome of the heartworm, Dirofilaria immitis, reveals drug and vaccine targets. FASEB J. 2012;26:4650–61.View ArticlePubMed CentralPubMedGoogle Scholar
- Schnall-Levin M, Zhao Y, Perrimon N, Berger B. Conserved microRNA targeting in Drosophila is as widespread in coding regions as in 3’ UTRs. Proc Natl Acad Sci U S A. 2010;107:15751–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Ha I, Wightman B, Ruvkun G. A bulged lin-4/lin-14 RNA duplex is sufficient for Caenorhabditis elegans lin-14 temporal gradient formation. Genes Dev. 1996;10:3041–50.View ArticlePubMedGoogle Scholar
- Chi SW, Zang JB, Mele A, Darnell RB. Argonaute HITS-CLIP decodes microRNA-mRNA interaction maps. Nature. 2009;460:479–86.PubMed CentralPubMedGoogle Scholar
- Mangone M, Manoharan AP, Thierry-Mieg D, Thierry-Mieg J, Han T, Mackowiak SD, et al. The landscape of C. elegans 3’UTRs. Science. 2010;329:432–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5:R1.View ArticlePubMed CentralPubMedGoogle Scholar
- Remm M, Storm CE, Sonnhammer EL. Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001;314:1041–52.View ArticlePubMedGoogle Scholar
- Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007;27:91–105.View ArticlePubMed CentralPubMedGoogle Scholar
- Im SH, Lee J. Identification of HMG-5 as a double-stranded telomeric DNA-binding protein in the nematode Caenorhabditis elegans. FEBS Lett. 2003;554:455–61.View ArticlePubMedGoogle Scholar
- Ford L, Zhang J, Liu J, Hashmi S, Fuhrman JA, Oksov Y, et al. Functional analysis of the cathepsin-like cysteine protease genes in adult Brugia malayi using RNA interference. PLoS Negl Trop Dis. 2009;3:e377.View ArticlePubMed CentralPubMedGoogle Scholar
- Singh PK, Kushwaha S, Mohd S, Pathak M, Misra-Bhattacharya S. In vitro gene silencing of independent phosphoglycerate mutase (iPGM) in the filarial parasite Brugia malayi. Infect Dis Poverty. 2013;2:5.View ArticlePubMed CentralPubMedGoogle Scholar
- Kushwaha S, Singh PK, Shahab M, Pathak M, Bhattacharya SM. In vitro silencing of Brugia malayi trehalose-6-phosphate phosphatase impairs embryogenesis and in vivo development of infective larvae in jirds. PLoS Negl Trop Dis. 2012;6:e1770.View ArticlePubMed CentralPubMedGoogle Scholar
- Buck AH, Blaxter M. Functional diversification of Argonautes in nematodes: an expanding universe. Biochem Soc Trans. 2013;41:881–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Landmann F, Foster JM, Slatko BE, Sullivan W. Efficient in vitro RNA interference and immunofluorescence-based phenotype analysis in a human parasitic nematode. Brugia malayi. Parasit Vectors. 2012;5:16.View ArticlePubMed CentralPubMedGoogle Scholar
- Winter AD, McCormack G, Myllyharju J, Page AP. Prolyl 4-hydroxlase activity is essential for development and cuticle formation in the human infective parasitic nematode Brugia malayi. J Biol Chem. 2013;288:1750–61.View ArticlePubMed CentralPubMedGoogle Scholar
- Wu B, Novelli J, Jiang D, Dailey HA, Landmann F, Ford L, et al. Interdomain lateral gene transfer of an essential ferrochelatase gene in human parasitic nematodes. Proc Natl Acad Sci U S A. 2013;110:7748–53.View ArticlePubMed CentralPubMedGoogle Scholar
- Brenner JL, Jasiewicz KL, Fahley AF, Kemp BJ, Abbott AL. Loss of individual microRNAs causes mutant phenotypes in sensitized genetic backgrounds in C. elegans. Curr Biol. 2010;20:1321–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Jan CH, Friedman RC, Ruby JG, Bartel DP. Formation, regulation and evolution of Caenorhabditis elegans 3’UTRs. Nature. 2011;469:97–101.View ArticlePubMed CentralPubMedGoogle Scholar
- Karp X, Hammell M, Ow MC, Ambros V. Effect of life history on microRNA expression during C. elegans development. RNA. 2011;17:639–51.View ArticlePubMed CentralPubMedGoogle Scholar
- Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115:787–98.View ArticlePubMedGoogle Scholar
- Kar A, Gutierrez-Hartmann A. Molecular mechanisms of ETS transcription factor-mediated tumorigenesis. Crit Rev Biochem Mol Biol. 2013;48:522–43.View ArticlePubMed CentralPubMedGoogle Scholar
- Albino D, Longoni N, Curti L, Mello-Grand M, Pinton S, Civenni G, et al. ESE3/EHF controls epithelial cell differentiation and its loss leads to prostate tumors with mesenchymal and stem-like features. Cancer Res. 2012;72:2889–900.View ArticlePubMedGoogle Scholar
- Zisoulis DG, Lovci MT, Wilbert ML, Hutt KR, Liang TY, Pasquinelli AE, et al. Comprehensive discovery of endogenous Argonaute binding sites in Caenorhabditis elegans. Nat Struct Mol Biol. 2010;17:173–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang L, Ding L, Cheung TH, Dong MQ, Chen J, Sewell AK, et al. Systematic identification of C. elegans miRISC proteins, miRNAs, and mRNA targets by their interactions with GW182 proteins AIN-1 and AIN-2. Mol Cell. 2007;28:598–613.View ArticlePubMed CentralPubMedGoogle Scholar
- Reeves R. Nuclear functions of the HMG proteins. Biochim Biophys Acta. 2010;1799:3–14.View ArticlePubMed CentralPubMedGoogle Scholar
- Lee YS, Dutta A. The tumor suppressor microRNA let-7 represses the HMGA2 oncogene. Genes Dev. 2007;21:1025–30.View ArticlePubMed CentralPubMedGoogle Scholar
- Hunter SJ, Thompson FJ, Tetley L, Devaney E. Temperature is a cue for gene expression in the post-infective L3 of the parasitic nematode Brugia pahangi. Mol Biochem Parasitol. 2001;112:1–9.View ArticlePubMedGoogle Scholar
- May P, Liao W, Wu Y, Shuai B, McCombie WR, Zhang MQ, et al. The effects of carbon dioxide and temperature on microRNA expression in Arabidopsis development. Nat Commun. 2013;4:2145.View ArticlePubMedGoogle Scholar
- Fu Y, Lan J, Wu X, Yang D, Zhang Z, Nie H, et al. Identification of Dirofilaria immitis miRNA using illumina deep sequencing. Vet Res. 2013;44:3.View ArticlePubMed CentralPubMedGoogle Scholar
- Dieterich C, Clifton SW, Schuster LN, Chinwalla A, Delehaunty K, Dinkelacker I, et al. The Pristionchus pacificus genome provides a unique perspective on nematode lifestyle and parasitism. Nat Genet. 2008;40:1193–8.View ArticlePubMedGoogle Scholar
- Devaney E, Jecock RM. The expression of the Mr 30000 antigen in the third stage larvae of Brugia pahangi. Parasite Immunol. 1991;13:75–87.View ArticlePubMedGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.View ArticlePubMed CentralPubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics. 2014;31:166–9.View ArticlePubMed CentralPubMedGoogle Scholar
- McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40:4288–97.View ArticlePubMed CentralPubMedGoogle Scholar
- Horwich MD, Zamore PD. Design and delivery of antisense oligonucleotides to block microRNA function in cultured Drosophila and human cells. Nat Protoc. 2008;3:1537–49.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.