Roz Laing

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

BackgroundThe 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.ResultsHere 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.ConclusionsThe H. contortus genome and transcriptome provide an essential platform for postgenomic research in this and other important strongylid parasites.

Characterization of the xenobiotic response of Caenorhabditis elegans to the anthelmintic drug albendazole and the identification of novel drug glucoside metabolites.

Knowledge of how anthelmintics are metabolized and excreted in nematodes is an integral part of understanding the factors that determine their potency, spectrum of activity and for investigating mechanisms of resistance. Although there is remarkably little information on these processes in nematodes, it is often suggested that they are of minimal importance for the major anthelmintic drugs. Consequently, we have investigated how the model nematode Caenorhabditis elegans responds to and metabolizes albendazole, one of the most important anthelmintic drugs for human and animal use. Using a mutant strain lacking the β-tubulin drug target to minimize generalized stress responses, we show that the transcriptional response is dominated by genes encoding XMEs (xenobiotic-metabolizing enzymes), particularly cytochrome P450s and UGTs (UDP-glucuronosyl transferases). The most highly induced genes are predominantly expressed in the worm intestine, supporting their role in drug metabolism. HPLC-MS/MS revealed the production of two novel glucoside metabolites in C. elegans identifying a major difference in the biotransformation of this drug between nematodes and mammals. This is the first demonstration of metabolism of a therapeutic anthelmintic in C. elegans and provides a framework for its use to functionally investigate nematode anthelmintic metabolism.

Ivermectin – Old Drug, New Tricks?

Ivermectin is one of the most important drugs in veterinary and human medicine for the control of parasitic infection and was the joint focus of the 2015 Nobel Prize in Physiology or Medicine, some 35 years after its remarkable discovery. Although best described for its activity on glutamate-gated chloride channels in parasitic nematodes, understanding of its mode of action remains incomplete. In the field of veterinary medicine, resistance to ivermectin is now widespread, but the mechanisms underlying resistance are unresolved. Here we discuss the history of this versatile drug and its use in global health. Based on recent studies in a variety of systems, we question whether ivermectin could have additional modes of action on parasitic nematodes.

The transcriptional response of Caenorhabditis elegans to Ivermectin exposure identifies novel genes involved in the response to reduced food intake.

We have examined the transcriptional response of Caenorhabditis elegans following exposure to the anthelmintic drug ivermectin (IVM) using whole genome microarrays and real-time QPCR. Our original aim was to identify candidate molecules involved in IVM metabolism and/or excretion. For this reason the IVM tolerant strain, DA1316, was used to minimise transcriptomic changes related to the phenotype of drug exposure. However, unlike equivalent work with benzimidazole drugs, very few of the induced genes were members of xenobiotic metabolising enzyme families. Instead, the transcriptional response was dominated by genes associated with fat mobilization and fatty acid metabolism including catalase, esterase, and fatty acid CoA synthetase genes. This is consistent with the reduction in pharyngeal pumping, and consequential reduction in food intake, upon exposure of DA1316 worms to IVM. Genes with the highest fold change in response to IVM exposure, cyp-37B1, mtl-1 and scl-2, were comparably up-regulated in response to short–term food withdrawal (4 hr) independent of IVM exposure, and GFP reporter constructs confirm their expression in tissues associated with fat storage (intestine and hypodermis). These experiments have serendipitously identified novel genes involved in an early response of C. elegans to reduced food intake and may provide insight into similar processes in higher organisms.

Annotation of Two Large Contiguous Regions from the Haemonchus contortus Genome Using RNA-seq and Comparative Analysis with Caenorhabditis elegans

The genomes of numerous parasitic nematodes are currently being sequenced, but their complexity and size, together with high levels of intra-specific sequence variation and a lack of reference genomes, makes their assembly and annotation a challenging task. Haemonchus contortus is an economically significant parasite of livestock that is widely used for basic research as well as for vaccine development and drug discovery. It is one of many medically and economically important parasites within the strongylid nematode group. This group of parasites has the closest phylogenetic relationship with the model organism Caenorhabditis elegans, making comparative analysis a potentially powerful tool for genome annotation and functional studies. To investigate this hypothesis, we sequenced two contiguous fragments from the H. contortus genome and undertook detailed annotation and comparative analysis with C. elegans. The adult H. contortus transcriptome was sequenced using an Illumina platform and RNA-seq was used to annotate a 409 kb overlapping BAC tiling path relating to the X chromosome and a 181 kb BAC insert relating to chromosome I. In total, 40 genes and 12 putative transposable elements were identified. 97.5% of the annotated genes had detectable homologues in C. elegans of which 60% had putative orthologues, significantly higher than previous analyses based on EST analysis. Gene density appears to be less in H. contortus than in C. elegans, with annotated H. contortus genes being an average of two-to-three times larger than their putative C. elegans orthologues due to a greater intron number and size. Synteny appears high but gene order is generally poorly conserved, although areas of conserved microsynteny are apparent. C. elegans operons appear to be partially conserved in H. contortus. Our findings suggest that a combination of RNA-seq and comparative analysis with C. elegans is a powerful approach for the annotation and analysis of strongylid nematode genomes.

Analysis of putative resistance gene loci in UK field populations of Haemonchus contortus after 6years of macrocyclic lactone use.

Graphical abstract

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.

A Genome Resequencing-Based Genetic Map Reveals the Recombination Landscape of an Outbred Parasitic Nematode in the Presence of Polyploidy and Polyandry

Abstract The parasitic nematode Haemonchus contortus is an economically and clinically important pathogen of small ruminants, and a model system for understanding the mechanisms and evolution of traits such as anthelmintic resistance. Anthelmintic resistance is widespread and is a major threat to the sustainability of livestock agriculture globally; however, little is known about the genome architecture and parameters such as recombination that will ultimately influence the rate at which resistance may evolve and spread. Here, we performed a genetic cross between two divergent strains of H. contortus, and subsequently used whole-genome resequencing of a female worm and her brood to identify the distribution of genome-wide variation that characterizes these strains. Using a novel bioinformatic approach to identify variants that segregate as expected in a pseudotestcross, we characterized linkage groups and estimated genetic distances between markers to generate a chromosome-scale F1 genetic map. We exploited this map to reveal the recombination landscape, the first for any helminth species, demonstrating extensive variation in recombination rate within and between chromosomes. Analyses of these data also revealed the extent of polyandry, whereby at least eight males were found to have contributed to the genetic variation of the progeny analyzed. Triploid offspring were also identified, which we hypothesize are the result of nondisjunction during female meiosis or polyspermy. These results expand our knowledge of the genetics of parasitic helminths and the unusual life-history of H. contortus, and enhance ongoing efforts to understand the genetic basis of resistance to the drugs used to control these worms and for related species that infect livestock and humans throughout the world. This study also demonstrates the feasibility of using whole-genome resequencing data to directly construct a genetic map in a single generation cross from a noninbred nonmodel organism with a complex lifecycle.

A major locus for ivermectin resistance in a parasitic nematode

Abstract Background Infections with helminths cause an enormous disease burden in billions of animals and plants worldwide. Large scale use of anthelmintics has driven the evolution of resistance in a number of species that infect livestock and companion animals, and there are growing concerns regarding the reduced efficacy in some human-infective helminths. Understanding the mechanisms by which resistance evolves is the focus of increasing interest; robust genetic analysis of helminths is challenging, and although many candidate genes have been proposed, the genetic basis of resistance remains poorly resolved. Results Here, we present a genome-wide analysis of two genetic crosses between ivermectin resistant and sensitive isolates of the parasitic nematode Haemonchus contortus, an economically important gastrointestinal parasite of small ruminants and a model for anthelmintic research. Whole genome sequencing of parental populations, and key stages throughout the crosses, identified extensive genomic diversity that differentiates populations, but after backcrossing and selection, a single genomic quantitative trait locus (QTL) localised on chromosome V was revealed to be associated with ivermectin resistance. This QTL was common between the two geographically and genetically divergent resistant populations and did not include any leading candidate genes, suggestive of a previously uncharacterised mechanism and/or driver of resistance. Despite limited resolution due to low recombination in this region, population genetic analyses and novel evolutionary models supported strong selection at this Q.TL, driven by at least partial dominance of the resistant allele, and that large resistance-associated haplotype blocks were enriched in response to selection. Conclusions We have described the genetic architecture and mode of ivermectin selection, revealing a major genomic locus associated with ivermectin resistance, the most conclusive evidence to date in any parasitic nematode. This study highlights a novel genome-wide approach to the analysis of a genetic cross in non-model organisms with extreme genetic diversity, and the importance of a high quality reference genome in interpreting the signals of selection so identified.Helminth infections cause an enormous disease burden in humans, livestock and companion animals. Helminth control is heavily dependent on the regular administration of anthelmintic drugs. Resistance to these drugs is widespread in many livestock systems, and there is concern that mass drug administration programs to control human helminths may be selecting for resistance. Here, we build on two previous genetic crosses between ivermectin resistant and sensitive isolates of the parasitic nematode Haemonchus contortus, an economically important gastrointestinal parasite of small ruminants that is a model for anthelmintic research. By combining genetic mapping and whole genome sequencing with novel analyses, we identify the first single genomic locus conclusively associated with ivermectin resistance in a parasitic nematode. This locus is common between two geographically and genetically divergent strains and does not include leading candidate genes. This work represents the first example of a major anthelmintic resistance locus being identified by genetic mapping in a parasitic nematode and represents the most comprehensive genetic analysis of this trait in a parasitic nematode to date.

Transcriptomic profiling of nematode parasites surviving vaccine exposure

Graphical abstract