Annotation of two large contiguous regions from the Haemonchus contortus genome using RNA-seq and comparative analysis with Caenorhabditis elegans

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

Toward integrative 'omics of the barber's pole worm and related parasitic nematodes

Advances in nucleic acid sequencing, mass spectrometry and computational biology have facilitated the identification, annotation and analysis of genes, transcripts, proteins and metabolites in model nematodes (Caenorhabditis elegans and Pristionchus pacificus) and socioeconomically important parasitic nematodes (Clades I, III, IV and V). Significant progress has been made in genomics and transcriptomics as well as in the proteomics and lipidomics of Haemonchus contortus (the barber's pole worm) - one of the most pathogenic representatives of the order Strongylida. Here, we review salient aspects of genomics, transcriptomics, proteomics, lipidomics, glycomics and functional genomics, and discuss the rise of integrative 'omics of this economically important parasite. Although our knowledge of the molecular biology, genetics and biochemistry of H. contortus and related species has progressed significantly, much remains to be explored, particularly in areas such as drug resistance, unique/unknown genes, host-parasite interactions, parasitism and the pathogenesis of disease, by integrating the use of multiple 'omics methods. This approach should lead to a better understanding of H. contortus and its relatives at a 'systems biology' level, and should assist in developing new interventions against these parasites.

A perspective on the discovery of selected compounds with anthelmintic activity against the barber's pole worm-Where to from here?

Parasitic roundworms (nematodes) cause substantial morbidity and mortality in animals worldwide. Anthelmintic treatment is central to controlling these worms, but widespread resistance to most of the commercially available anthelmintics for veterinary and agricultural use is compromising control, such that there is an urgency to discover new and effective drugs. The purpose of this article is to review information on parasitic nematodes, the treatment and control of parasitic nematode infections and aspects of discovering new anthelmintics in the context of anthelmintic resistance problems, and then to discuss some progress that our group has made in identifying selected compounds with activity against nematodes. The focus of our recent work has been on discovering new chemical entities and known drugs with anthelmintic activities against Haemonchus contortus as well as other socioeconomically important parasitic nematodes for subsequent development. Using whole worm-based phenotypic assays, we have been screening compound collections obtained via product-development-partnerships and/or collaborators, and active compounds have been assessed for their potential as anthelmintic candidates. Following the screening of 15,333 chemicals from five distinct compound collections against H. contortus, we have discovered one new chemical entity (designated SN00797439), two human kinase inhibitors (SNS-032 and AG-1295), 14 tetrahydroquinoxaline analogues, one insecticide (tolfenpyrad) and two tolfenpyrad (pyrazole-5-carboxamide) derivatives (a-15 and a-17) with anthelmintic activity in vitro. Some of these 20 'hit' compounds have selectivity against H. contortus in vitro when compared to particular human cell lines. In our opinion, some of these compounds could represent starting points for 'lead' development. Accordingly, the next research steps to be pursued include: (i) chemical optimisation of representative chemicals via structure-activity relationship (SAR) evaluations; (ii) assessment of the breadth of spectrum of anthelmintic activity on a range of other parasitic nematodes, such as strongyloids, ascaridoids, enoplids and filarioids; (iii) detailed investigations of the absorption, distribution, metabolism, excretion and toxicity (ADMET) of optimised chemicals with broad nematocidal or nematostatic activity; and (iv) establishment of the modes of action of lead candidates.

The confounding effects of high genetic diversity on the determination and interpretation of differential gene expression analysis in the parasitic nematode Haemonchus contortus

Differential expression analysis between parasitic nematode strains is commonly used to implicate candidate genes in anthelmintic resistance or other biological functions. We have tested the hypothesis that the high genetic diversity of an organism such as Haemonchus contortus could complicate such analyses. First, we investigated the extent to which sequence polymorphism affects the reliability of differential expression analysis between the genetically divergent H. contortus strains MHco3(ISE), MHco4(WRS) and MHco10(CAVR). Using triplicates of 20 adult female worms from each population isolated under parallel experimental conditions, we found that high rates of sequence polymorphism in RNAseq reads were associated with lower efficiency read mapping to gene models under default TopHat2 parameters, leading to biased estimates of inter-strain differential expression. We then showed it is possible to largely compensate for this bias by optimising the read mapping single nucleotide polymorphism (SNP) allowance and filtering out genes with particularly high single nucleotide polymorphism rates. Once the sequence polymorphism biases were removed, we then assessed the genuine transcriptional diversity between the strains, finding ≥824 differentially expressed genes across all three pairwise strain comparisons. This high level of inter-strain transcriptional diversity not only suggests substantive inter-strain phenotypic variation but also highlights the difficulty in reliably associating differential expression of specific genes with phenotypic differences. To provide a practical example, we analysed two gene families of potential relevance to ivermectin drug resistance; the ABC transporters and the ligand-gated ion channels (LGICs). Over half of genes identified as differentially expressed using default TopHat2 parameters were shown to be an artifact of sequence polymorphism differences. This work illustrates the need to account for sequence polymorphism in differential expression analysis. It also demonstrates that a large number of genuine transcriptional differences can occur between H. contortus strains and these must be considered before associating the differential expression of specific genes with phenotypic differences between strains.

Probing the opportunities for designing anthelmintic leads by sub-structural topology-based QSAR modelling

A quantitative structure-activity (QSAR) model has been developed for enriched tubulin inhibitors, which were retrieved from sequence similarity searches and applicability domain analysis. Using partial least square (PLS) method and leave-one-out (LOO) validation approach, the model was generated with the correlation statistics of [Formula: see text] and [Formula: see text] of 0.68 and 0.69, respectively. The present study indicates that topological descriptors, viz. BIC, CH_3_C, IC, JX and Kappa_2 correlate well with biological activity. ADME and toxicity (or ADME/T) assessment showed that out of 260 molecules, 255 molecules successfully passed the ADME/T assessment test, wherein the drug-likeness attributes were exhibited. These results showed that topological indices and the colchicine binding domain directly influence the aetiology of helminthic infections. Further, we anticipate that our model can be applied for guiding and designing potential anthelmintic inhibitors.

Perusal of parasitic nematode 'omics in the post-genomic era

The advent of high-throughput, next-generation sequencing methods combined with advances in computational biology and bioinformatics have greatly accelerated discovery within biomedical research. This "post-genomics" era has ushered in powerful approaches allowing one to quantify RNA transcript and protein abundance for every gene in the genome - often for multiple conditions. Herein, we chronicle how the post-genomics era has advanced our overall understanding of parasitic nematodes through transcriptomics and proteomics and highlight some of the important advances made in each major nematode clade. We primarily focus on organisms relevant to human health, given that nematode infections significantly impact disability-adjusted life years (DALY) scores within the developing world, but we also discuss organisms of veterinary importance as well as those used as laboratory models. As such, we envision that this review will serve as a comprehensive resource for those seeking a better understanding of basic parasitic nematode biology as well as those interested in targets for vaccination and pharmacological intervention.

In silico analysis of the binding of anthelmintics to Caenorhabditis elegansP-glycoprotein 1

Macrocyclic lactones (ML) are important anthelmintics used in animals and humans against parasite nematodes, but their therapeutic success is compromised by the spread of ML resistance. Some ABC transporters, such as p-glycoproteins (Pgps), are selected and overexpressed in ML-resistant nematodes, supporting a role for some drug efflux proteins in ML resistance. However, the role of such proteins in ML transport remains to be clarified at the molecular level. Recently, Caenorhabditis elegans Pgp-1 (Cel-Pgp-1) has been crystallized, and its drug-modulated ATPase function characterized in vitro revealed Cel-Pgp-1 as a multidrug transporter. Using this crystal structure, we have developed an in silico drug docking model in order to study the binding of ML and other anthelmintic drugs to Cel-Pgp-1. All tested ML bound with high affinity in a unique site, within the inner chamber of the protein, supporting that ML may be transported by Cel-Pgp-1. Interestingly, interacting residues delineate a ML specific fingerprint involving H-bonds, including T1028. In particular, benzofurane and spiroketal moieties bound to specific sub-sites. When compared with the aglycone ML, such as moxidectin and ivermectin aglycone, avermectin anthelmintics have significant higher affinity for Cel-Pgp-1, likely due to the sugar substituent(s) that bind to a specific area involving H-bonds at Y771. Triclabendazole, closantel and emodepside bound with good affinities to different sub-sites in the inner chamber, partially overlapping with the ML binding site, suggesting that they could compete for Cel-Pgp-1-mediated ML transport. In conclusion, this work provides novel information on the role of nematode Pgps in transporting anthelmintics, and a valuable tool to predict drug-drug interactions and to rationally design new competitive inhibitors of clinically-relevant nematode Pgps, to improve anthelmintic therapeutics.

Evidence from two independent backcross experiments supports genetic linkage of microsatellite Hcms8a20, but not other candidate loci, to a major ivermectin resistance locus in Haemonchus contortus

Haemonchus contortus is the leading parasitic nematode species used to study anthelmintic drug resistance. A variety of candidate loci have been implicated as being associated with ivermectin resistance in this parasite but definitive evidence of their importance is still lacking. We have previously performed two independent serial backcross experiments to introgress ivermectin resistance loci from two H. contortus ivermectin-resistant strains - MHco4(WRS) and MHco10(CAVR) - into the genetic background of the ivermectin-susceptible genome reference strain MHco3(ISE). We have interrogated a number of candidate ivermectin resistance loci in the resulting backcross populations and assessed the evidence for their genetic linkage to an ivermectin resistance locus. These include the microsatellite marker Hcms8a20 and six candidate genes Hco-glc-5, Hco-avr-14, Hco-lgc-37 (previously designated Hco-hg-1), Hco-pgp-9 (previously designated Hco-pgp-1), Hco-pgp-2 and Hco-dyf-7. We have sampled the haplotype diversity of amplicon markers within, or adjacent to, each of these loci in the parental strains and fourth generation backcross populations to assess the evidence for haplotype introgression from the resistant parental strain into the genomic background of the susceptible parental strain in each backcross. The microsatellite Hcms8a20 locus showed strong evidence of such introgression in both independent backcrosses, suggesting it is linked to an important ivermectin resistance mutation in both the MHco4(WRS) and MHco10(CAVR) strains. In contrast, Hco-glc-5, Hco-avr-14, Hco-pgp-9 and Hco-dyf-7 showed no evidence of introgression in either backcross. Hco-lgc-37 and Hco-pgp-2 showed only weak evidence of introgression in the MHco3/4 backcross but not in the MHco3/10 backcross. Overall, these results suggest that microsatellite marker Hcms8a20, but not the other candidate genes tested, is linked to a major ivermectin resistance locus in the MHco4(WRS) and MHco10(CAVR) strains. This work also emphasises the need for genome-wide approaches to identify mutations responsible for the ivermectin resistance in this parasite.

Haemonchus contortus: Genome Structure, Organization and Comparative Genomics

One of the first genome sequencing projects for a parasitic nematode was that for Haemonchus contortus. The open access data from the Wellcome Trust Sanger Institute provided a valuable early resource for the research community, particularly for the identification of specific genes and genetic markers. Later, a second sequencing project was initiated by the University of Melbourne, and the two draft genome sequences for H. contortus were published back-to-back in 2013. There is a pressing need for long-range genomic information for genetic mapping, population genetics and functional genomic studies, so we are continuing to improve the Wellcome Trust Sanger Institute assembly to provide a finished reference genome for H. contortus. This review describes this process, compares the H. contortus genome assemblies with draft genomes from other members of the strongylid group and discusses future directions for parasite genomics using the H. contortus model.

The evolution of pentameric ligand-gated ion-channels and the changing family of anthelmintic drug targets

SUMMARY Pentameric ligand-gated ion-channels rapidly transduce synaptic neurotransmitter signals to an electrical response in post-synaptic neuronal or muscle cells and control the neuromusculature of a majority of multicellular animals. A wide range of pharmaceuticals target these receptors including ethanol, nicotine, anti-depressants and other mood regulating drugs, compounds that control pain and mobility and are targeted by a majority of anthelmintic drugs used to control parasitic infection of humans and livestock. Major advances have been made in recent years to our understanding of the structure, function, activity and the profile of compounds that can activate specific receptors. It is becoming clear that these anthelmintic drug targets are not fixed, but differ in significant details from one nematode species to another. Here we review what is known about the evolution of the pentameric ligand-gated ion-channels, paying particular attention to the nematodes, how we can infer the origins of such receptors and understand the factors that determine how they change both over time and from one species to another. Using this knowledge provides a biological framework in which to understand these important drug targets and avenues to identify new receptors and aid the search for new anthelmintic drugs.

Improved structural annotation of protein-coding genes in the Meloidogyne hapla genome using RNA-Seq

As high-throughput cDNA sequencing (RNA-Seq) is increasingly applied to hypothesis-driven biological studies, the prediction of protein coding genes based on these data are usurping strictly in silico approaches. Compared with computationally derived gene predictions, structural annotation is more accurate when based on biological evidence, particularly RNA-Seq data. Here, we refine the current genome annotation for the Meloidogyne hapla genome utilizing RNA-Seq data. Published structural annotation defines 14 420 protein-coding genes in the M. hapla genome. Of these, 25% (3751) were found to exhibit some incongruence with RNA-Seq data. Manual annotation enabled these discrepancies to be resolved. Our analysis revealed 544 new gene models that were missing from the prior annotation. Additionally, 1457 transcribed regions were newly identified on the ends of as-yet-unjoined contigs. We also searched for trans-spliced leaders, and based on RNA-Seq data, identified genes that appear to be trans-spliced. Four 22-bp trans-spliced leaders were identified using our pipeline, including the known trans-spliced leader, which is the M. hapla ortholog of SL1. In silico predictions of trans-splicing were validated by comparison with earlier results derived from an independent cDNA library constructed to capture trans-spliced transcripts. The new annotation, which we term HapPep5, is publically available at www.hapla.org.

Haemonchus contortus as a paradigm and model to study anthelmintic drug resistance

Anthelmintic resistance is a major problem for the control livestock parasites and a potential threat to the sustainability of community-wide treatment programmes being used to control human parasites in the developing world. Anthelmintic resistance is essentially a complex quantitative trait in which multiple mutations contribute to the resistance phenotype in an additive manner. Consequently, a combination of forward genetic and genomic approaches are needed to identify the causal mutations and quantify their contribution to the resistance phenotype. Therefore, there is a need to develop genetic and genomic approaches for key parasite species identified as relevant models. Haemonchus contortus, a gastro-intestinal parasite of sheep, has shown a remarkable propensity to develop resistance to all the drugs used in its control. Partly because of this, and partly because of its experimental amenability, research on this parasite has contributed more than any other to our understanding of anthelmintic resistance. H. contortus offers a variety of advantages as an experimental system including the ability to undertake genetic crosses; a prerequisite for genetic mapping. This review will discuss the current progress on developing H. contortus as a model system in which to study anthelmintic resistance.

The genome and transcriptome of Haemonchus contortus, a key model parasite for drug and vaccine discovery

BACKGROUND

The small ruminant parasite Haemonchus contortus is the most widely used parasitic nematode in drug discovery, vaccine development and anthelmintic resistance research. Its remarkable propensity to develop resistance threatens the viability of the sheep industry in many regions of the world and provides a cautionary example of the effect of mass drug administration to control parasitic nematodes. Its phylogenetic position makes it particularly well placed for comparison with the free-living nematode Caenorhabditis elegans and the most economically important parasites of livestock and humans.

RESULTS

Here we report the detailed analysis of a draft genome assembly and extensive transcriptomic dataset for H. contortus. This represents the first genome to be published for a strongylid nematode and the most extensive transcriptomic dataset for any parasitic nematode reported to date. We show a general pattern of conservation of genome structure and gene content between H. contortus and C. elegans, but also a dramatic expansion of important parasite gene families. We identify genes involved in parasite-specific pathways such as blood feeding, neurological function, and drug metabolism. In particular, we describe complete gene repertoires for known drug target families, providing the most comprehensive understanding yet of the action of several important anthelmintics. Also, we identify a set of genes enriched in the parasitic stages of the lifecycle and the parasite gut that provide a rich source of vaccine and drug target candidates.

CONCLUSIONS

The H. contortus genome and transcriptome provide an essential platform for postgenomic research in this and other important strongylid parasites.

Characterization and comparative analysis of the complete Haemonchus contortus β-tubulin gene family and implications for benzimidazole resistance in strongylid nematodes

Parasitic nematode β-tubulin genes are of particular interest because they are the targets of benzimidazole drugs. However, in spite of this, the full β-tubulin gene family has not been characterized for any parasitic nematode to date. Haemonchus contortus is the parasite species for which we understand benzimidazole resistance the best and its close phylogenetic relationship with Caenorhabditis elegans potentially allows inferences of gene function by comparative analysis. Consequently, we have characterized the full β-tubulin gene family in H. contortus. Further to the previously identified Hco-tbb-iso-1 and Hco-tbb-iso-2 genes, we have characterized two additional family members designated Hco-tbb-iso-3 and Hco-tbb-iso-4. We show that Hco-tbb-iso-1 is not a one-to-one orthologue with Cel-ben-1, the only β-tubulin gene in C. elegans that is a benzimidazole drug target. Instead, both Hco-tbb-iso-1 and Hco-tbb-iso-2 have a complex evolutionary relationship with three C. elegans β-tubulin genes: Cel-ben-1, Cel-tbb-1 and Cel-tbb-2. Furthermore, we show that both Hco-tbb-iso-1 and Hco-tbb-iso-2 are highly expressed in adult worms; in contrast, Hco-tbb-iso-3 and Hco-tbb-iso-4 are expressed only at very low levels and are orthologous to the Cel-mec-7 and Cel-tbb-4 genes, respectively, suggesting that they have specialized functional roles. Indeed, we have found that the expression pattern of Hco-tbb-iso-3 in H. contortus is identical to that of Cel-mec-7 in C. elegans, being expressed in just six "touch receptor" mechano-sensory neurons. These results suggest that further investigation is warranted into the potential involvement of strongylid isotype-2 β-tubulin genes in mechanisms of benzimidazole resistance.

RNAseq analysis of the parasitic nematode Strongyloides stercoralis reveals divergent regulation of canonical dauer pathways

The infectious form of many parasitic nematodes, which afflict over one billion people globally, is a developmentally arrested third-stage larva (L3i). The parasitic nematode Strongyloides stercoralis differs from other nematode species that infect humans, in that its life cycle includes both parasitic and free-living forms, which can be leveraged to investigate the mechanisms of L3i arrest and activation. The free-living nematode Caenorhabditis elegans has a similar developmentally arrested larval form, the dauer, whose formation is controlled by four pathways: cyclic GMP (cGMP) signaling, insulin/IGF-1-like signaling (IIS), transforming growth factor β (TGFβ) signaling, and biosynthesis of dafachronic acid (DA) ligands that regulate a nuclear hormone receptor. We hypothesized that homologous pathways are present in S. stercoralis, have similar developmental regulation, and are involved in L3i arrest and activation. To test this, we undertook a deep-sequencing study of the polyadenylated transcriptome, generating over 2.3 billion paired-end reads from seven developmental stages. We constructed developmental expression profiles for S. stercoralis homologs of C. elegans dauer genes identified by BLAST searches of the S. stercoralis genome as well as de novo assembled transcripts. Intriguingly, genes encoding cGMP pathway components were coordinately up-regulated in L3i. In comparison to C. elegans, S. stercoralis has a paucity of genes encoding IIS ligands, several of which have abundance profiles suggesting involvement in L3i development. We also identified seven S. stercoralis genes encoding homologs of the single C. elegans dauer regulatory TGFβ ligand, three of which are only expressed in L3i. Putative DA biosynthetic genes did not appear to be coordinately regulated in L3i development. Our data suggest that while dauer pathway genes are present in S. stercoralis and may play a role in L3i development, there are significant differences between the two species. Understanding the mechanisms governing L3i development may lead to novel treatment and control strategies.

Parasitic nematodes infect over one billion people worldwide and cause many diseases, including strongyloidiasis, filariasis, and hookworm disease. For many of these parasites, including Strongyloides stercoralis, the infectious form is a developmentally arrested and long-lived thirdstage larva (L3i). Upon encountering a host, L3i quickly resume development and mature into parasitic adults. In the free-living nematode Caenorhabditis elegans, a similar developmentally arrested third-stage larva, known as the dauer, is regulated by four key cellular mechanisms. We hypothesized that similar cellular mechanisms control L3i arrest and activation. Therefore, we used deep-sequencing technology to characterize the S. stercoralis transcriptome (RNAseq), which allowed us to identify S. stercoralis homologs of components of these four mechanisms and examine their temporal regulation. We found similar temporal regulation between S. stercoralis and C. elegans for components of two mechanisms, but dissimilar temporal regulation for two others, suggesting conserved as well as novel modes of developmental regulation for L3i. Understanding L3i development may lead to novel control strategies as well as new treatments for strongyloidiasis and other diseases caused by parasitic nematodes.

Introgression of ivermectin resistance genes into a susceptible Haemonchus contortus strain by multiple backcrossing

Anthelmintic drug resistance in livestock parasites is already widespread and in recent years there has been an increasing level of anthelmintic drug selection pressure applied to parasitic nematode populations in humans leading to concerns regarding the emergence of resistance. However, most parasitic nematodes, particularly those of humans, are difficult experimental subjects making mechanistic studies of drug resistance extremely difficult. The small ruminant parasitic nematode Haemonchus contortus is a more amenable model system to study many aspects of parasite biology and investigate the basic mechanisms and genetics of anthelmintic drug resistance. Here we report the successful introgression of ivermectin resistance genes from two independent ivermectin resistant strains, MHco4(WRS) and MHco10(CAVR), into the susceptible genome reference strain MHco3(ISE) using a backcrossing approach. A panel of microsatellite markers were used to monitor the procedure. We demonstrated that after four rounds of backcrossing, worms that were phenotypically resistant to ivermectin had a similar genetic background to the susceptible reference strain based on the bulk genotyping with 18 microsatellite loci and individual genotyping with a sub-panel of 9 microsatellite loci. In addition, a single marker, Hcms8a20, showed evidence of genetic linkage to an ivermectin resistance-conferring locus providing a starting point for more detailed studies of this genomic region to identify the causal mutation(s). This work presents a novel genetic approach to study anthelmintic resistance and provides a “proof-of-concept” of the use of forward genetics in an important model strongylid parasite of relevance to human hookworms. The resulting strains provide valuable resources for candidate gene studies, whole genome approaches and for further genetic analysis to identify ivermectin resistance loci.

The use of drugs (anthelmintics) to control nematode parasites (roundworms) is common in both humans and animals. This has led to the widespread development of drug resistance in livestock parasites and serious concerns regarding its emergence in human parasites. Haemonchus contortus is a parasitic nematode of sheep that has a high propensity to develop resistance and is the most widely used model system in which to study anthelmintic drug resistance. Ivermectin is an extremely important drug for parasite control in both humans and animals. Here, we report a novel approach using genetic crossing to transfer a region of the H. contortus genome containing ivermectin resistance genes from resistant strains into a susceptible strain. During our backcrossing approach, we have identified a genetic marker showing evidence of genetic linkage to ivermectin resistance. The susceptible strain we have used is currently having its complete genome sequenced making the information and strains generated here extremely valuable for the identification of ivermectin resistance genes. This work represents an important proof of concept for using genetic approaches to identify genomic regions containing drug resistant genes in parasitic nematodes.

Diversity in parasitic nematode genomes: the microRNAs of Brugia pahangi and Haemonchus contortus are largely novel

Background

MicroRNAs (miRNAs) play key roles in regulating post-transcriptional gene expression and are essential for development in the free-living nematode Caenorhabditis elegans and in higher organisms. Whether microRNAs are involved in regulating developmental programs of parasitic nematodes is currently unknown. Here we describe the the miRNA repertoire of two important parasitic nematodes as an essential first step in addressing this question.

Results

The small RNAs from larval and adult stages of two parasitic species, Brugia pahangi and Haemonchus contortus, were identified using deep-sequencing and bioinformatic approaches. Comparative analysis to known miRNA sequences reveals that the majority of these miRNAs are novel. Some novel miRNAs are abundantly expressed and display developmental regulation, suggesting important functional roles. Despite the lack of conservation in the miRNA repertoire, genomic positioning of certain miRNAs within or close to specific coding genes is remarkably conserved across diverse species, indicating selection for these associations. Endogenous small-interfering RNAs and Piwi-interacting (pi)RNAs, which regulate gene and transposon expression, were also identified. piRNAs are expressed in adult stage H. contortus, supporting a conserved role in germline maintenance in some parasitic nematodes.

Conclusions

This in-depth comparative analysis of nematode miRNAs reveals the high level of divergence across species and identifies novel sequences potentially involved in development. Expression of novel miRNAs may reflect adaptations to different environments and lifestyles. Our findings provide a detailed foundation for further study of the evolution and function of miRNAs within nematodes and for identifying potential targets for intervention.