Parasitic nematodes are highly successful pathogens, inflicting disease on humans, animals and plants. Despite great differences in their life cycles, host preference and transmission modes, these parasites share a common capacity to manipulate their host's immune system. This is at least partly achieved through the release of excretory/secretory proteins, the most well-characterized component of nematode secretomes, that are comprised of functionally diverse molecules. In this work, we analyzed published protein secretomes of parasitic nematodes to identify common patterns as well as species-specific traits. The 20 selected organisms span 4 nematode clades, including plant pathogens, animal parasites, and the free-living species Caenorhabditis elegans. Transthyretin-like proteins were the only component common to all adult secretomes; many other protein classes overlapped across multiple datasets. The glycolytic enzymes aldolase and enolase were present in all parasitic species, but missing from C. elegans. Secretomes from larval stages showed less overlap between species. Although comparison of secretome composition across species and life-cycle stages is challenged by the use of different methods and depths of sequencing among studies, our workflow enabled the identification of conserved protein families and pinpointed elements that may have evolved as to enable parasitism. This strategy, extended to more secretomes, may be exploited to prioritize therapeutic targets in the future.

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.

Haemonchus contortus, one of the most economically important parasites of small ruminants, has become resistant to the anthelmintic ivermectin. Deciphering the role of P-glycoproteins in ivermectin resistance is desirable for understanding and overcoming this resistance. In the model nematode, Caenorhabditis elegans, P-glycoprotein-13 is expressed in the amphids, important neuronal structures for ivermectin activity. We have focused on its ortholog in the parasite, Hco-Pgp-13. A 3D model of Hco-Pgp-13, presenting an open inward-facing conformation, has been constructed by homology with the Cel-Pgp-1 crystal structure. In silico docking calculations predicted high affinity binding of ivermectin and actinomycin D to the inner chamber of the protein. Following in vitro expression, we showed that ivermectin and actinomycin D modulated Hco-Pgp-13 ATPase activity with high affinity. Finally, we found in vivo Hco-Pgp-13 localization in epithelial, pharyngeal and neuronal tissues. Taken together, these data suggest a role for Hco-Pgp-13 in ivermectin transport, which could contribute to anthelmintic resistance.

Background: Hemicelluloses are a diverse group of complex, non-cellulosic polysaccharides, which constitute approximately one-third of the plant cell wall and find use as dietary fibres, food additives and raw materials for biofuels. Genes involved in hemicellulose synthesis have not been extensively studied in small grain cereals. Results: In efforts to isolate the sequences for the cellulose synthase-like (Csl) gene family from wheat, we identified 108 genes (hereafter referred to as TaCsl). Each gene was represented by two to three homeoalleles, which are named as TaCslXY_ZA, TaCslXY_ZB, or TaCslXY_ZD, where X denotes the Csl subfamily, Y the gene number and Z the wheat chromosome where it is located. A quarter of these genes were predicted to have 2 to 3 splice variants, resulting in a total of 137 putative translated products. Approximately 45% of TaCsl genes were located on chromosomes 2 and 3. Sequences from the subfamilies C and D were interspersed between the dicots and grasses but those from subfamily A clustered within each group of plants. Proximity of the dicot-specific subfamilies B and G, to the grass-specific subfamilies H and J, respectively, points to their common origin. In silico expression analysis in different tissues revealed that most of the genes were expressed ubiquitously and some were tissue-specific. More than half of the genes had introns in phase 0, one-third in phase 2, and a few in phase 1. Conclusion: Detailed characterization of the wheat Csl genes has enhanced the understanding of their structural, functional, and evolutionary features. This information will be helpful in designing experiments for genetic manipulation of hemicellulose synthesis with the goal of developing improved cultivars for biofuel production and increased tolerance against various stresses.

Anthelmintic resistance is a serious problem for the control of equine gastrointestinal nematodes. In the present survey, 173 third stage larvae of cyathostomins were investigated from three different locations for the presence of the resistant genotype at codon 167 of the beta-tubulin gene, as this is the most prevalent mutation. The larvae from the state of Parana (n=67), Sao Paulo (n=54) and Santa Catarina (n=52), showed 61.2; 31.5 and 38.5% of the heterozygous resistant genotype - TTC/TAC, respectively. An unpublished mutation at codon 172 that results in a serine (S) to threonine (T) substitution was found in 17.9% (12/67) of samples from Parana; and 13.0% (7/54) of samples from Sao Paulo. We have compared the molecular diagnostic with the fecal egg count data (R2=−0.79) from the same farms, and consider that the use of routine molecular diagnostic in individual larva may help to determine the population genetic distribution that is associated with drug failure.

Human onchocerciasis is a serious neglected tropical disease caused by the filarial nematode Onchocerca volvulus that can lead to blindness and chronic disability. Control of the disease relies largely on mass administration of a single drug, and the development of new drugs and vaccines depends on a better knowledge of parasite biology. Here, we describe the chromosomes of O. volvulus and its Wolbachia endosymbiont. We provide the highest-quality sequence assembly for any parasitic nematode to date, giving a glimpse into the evolution of filarial parasite chromosomes and proteomes. This resource was used to investigate gene families with key functions that could be potentially exploited as targets for future drugs. Using metabolic reconstruction of the nematode and its endosymbiont, we identified enzymes that are likely to be essential for O. volvulus viability. In addition, we have generated a list of proteins that could be targeted by Federal-Drug-Agency-approved but repurposed drugs, providing starting points for anti-onchocerciasis drug development.

Monepantel is a recently developed anthelmintic with a novel mode of action. Parasitic nematodes with reduced sensitivity to monepantel have led to the identification of MPTL-1, a ligand-gated ion-channel subunit of the parasitic nematode Haemonchus contortus, as a potential drug target. Homomeric MPTL-1 channels reconstituted in Xenopus oocytes are gated by µM concentrations of betaine and mM concentrations of choline. Measurement of reversal potentials indicated that the channel has a similar conductance for Na(+) and K(+) ions and does not permeate Ca(2+). Concentrations of monepantel (amino-acetonitrile derivative [AAD]-2225) >0.1 μM, but not its inactive enantiomer AAD-2224, induced channel opening in an irreversible manner. Currents elicited by monepantel alone were larger than the maximal current amplitudes achieved with betaine or choline, making monepantel a superagonist. Currents elicited by betaine or choline were allosterically potentiated by nM concentrations of monepantel and to a much smaller degree by AAD-2224. We have also reconstituted the Caenorhabditis elegans homomeric ACR-20 receptor in Xenopus oocytes. The acr-20 sequence has higher similarity to mptl-1 than acr-23, the primary target for monepantel mode of action in C. elegans. The ACR-20 channel is gated similarly as MPTL-1. Monepantel, but not AAD-2224, was able to induce channel opening in an irreversible manner at similar concentrations as for MPTL-1. Interestingly, the allosteric potentiation measured in the presence of betaine was much smaller than in MPTL-1 receptors. Together, these results establish the mode of action of monepantel in H. contortus and contribute to our understanding of the mode of action of this anthelmintic.

Haemonchus contortus is an abomasal nematode of ruminants that is widely present across the world. Its ability to cause death of infected animals and rapidly develop anthelmintic resistance makes it a dangerous pathogen. Ivermectin (IVM) and moxidectin (MOX) are macrocyclic lactones (MLs). They have been successfully used to treat parasitic nematodes over the last three decades. A genetic association between IVM selection and single nucleotide polymorphisms (SNPs) on the β-tubulin isotype 1 gene was reported in H. contortus. These SNPs result in replacing phenylalanine (F, TTC) with tyrosine (Y, TAC) at position 167 or 200 on the β-tubulin protein. Recently we reported a direct interaction of IVM with α- and β-tubulin. It had been hypothesized that the SNPs (F167Y and F200Y) may change tubulin dynamics and directly affect IVM binding. The goal of the current study was to observe the effects of SNPs (F167Y and F200Y) on tubulin polymerization and IVM binding. It was also of interest to evaluate the differences between IVM and MOX on tubulin polymerization. We conclude that the SNPs cause no difference in the polymerization of wild and mutant tubulins. Furthermore, neither of the SNPs reduced IVM binding. Varying results were obtained in the degree of polymerization of parasitic and mammalian tubulin for IVM and MOX, i.e., the extent of polymerization was greater for IVM compared with MOX, for H. contortus tubulin, and vice versa for mammalian tubulin. Molecular modeling showed that IVM and MOX docked into the taxane binding pocket of both mammalian and parasitic wild type and mutant tubulins. However the binding was stronger for mammalian tubulin as compared to parasitic tubulin.

Heavy reliance on macrocyclic lactones to treat parasitic nematodes has resulted in the evolution of widespread drug resistance that threatens human and animal health. Management strategies have been proposed that would slow the rise of resistance, however testing these strategies has been hampered by the lack of identified strong-effect resistance markers in parasites. We show that the Caenorhabditis elegans gene Cel_dyf-7, necessary for amphid sensory neuron development, also confers macrocyclic lactone sensitivity. In the sheep parasite Haemonchus contortus: (i) strains selected for macrocyclic lactone resistance were enriched in a Hco_dyf-7 haplotype that was rare in the drug-naïve population, (ii) the resistant haplotype correlated with the sensory neuron defects, and (iii) the resistant haplotype was associated with decreased Hco_dyf-7 expression. Resistant field isolates of H. contortus from five continents were enriched for the resistant haplotype, demonstrating the relevance of the Hco_dyf-7 haplotype to practise and indicating that it is a locus of strong effect. Hemizygosity resulting from sex linkage of dyf-7 likely contributes to the rise of resistance in treated populations.

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.

Monepantel is a member of the recently identified class of anthelmintics known as the amino-acetonitrile derivatives (AADs). Monepantel controls all major gastro-intestinal nematodes in sheep including those that are resistant to the classical anthelmintics. Previous studies have shown that the Caenorhabditis elegans acr-23 and the Haemonchus contortus Hco-mptl-1 genes may be prominent targets of monepantel. With this discovery it became possible to investigate the mode of action of monepantel in nematodes at the molecular level. In the present study, we show that a C. elegans mutant acr-23 strain is fully rescued by expressing the wild-type acr-23 gene. Moreover, we present a new mutant allele, and characterize acr-23 alleles genetically. We also show that acr-23 is expressed in body wall muscle cells, and provide therefore a possible explanation for the paralysis caused by monepantel. Furthermore, genetic evidence suggests that the chaperone RIC-3 is required for expression of full monepantel resistance. Finally, we present reconstitution of the C. elegans ACR-23 receptor in Xenopus laevis oocytes and provide direct evidence of its modulation by monepantel. Conversely, co-injection of the chaperone RIC-3 had no impact for channel reconstitution in X. laevis oocytes. These results reinforce the involvement of the ACR-23 family in the mode of action of monepantel and advance our understanding of this new class of anthelmintics.

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.

Resistance to macrocyclic lactones in the strongylid parasite of sheep, Haemonchus contortus, is associated with specific alleles of several genes, including those encoding ligand-gated chloride-channels. While previous functional studies of the ion-channels encoded by these resistant alleles have revealed alterations in ligand binding and response to the anthelmintics, we still do not know how these alleles are contributing to resistance in vivo. To understand the phenomenon of anthelmintic resistance in detail we need to link changes in the genes of individual parasites with their ability to resist the effects of anthelmintic exposure. We have determined the genotype of individual adult and larval H. contortus with respect to the glc-5 and lgc-37 genes linked with macrocyclic lactone resistance. In these same individuals, we measured feeding and movement, two characteristics targeted by the drug, which are believed to contribute to parasite killing. Both genes are linked with altered behavior in the absence of drug, providing evidence that genetic resistance may be associated with alterations in parasite biology. In the presence of macrocyclic lactones, both genes are associated with a degree of protection against drug action. Whether this protection may be effective under the application of anthelmintic as part of normal farm practice and whether this could explain the evolution of resistance remain unclear.

Mycoparasitism by antagonistic fungi involves changes in the biochemistry and physiology of both partners. Analysis of genes that are expressed during mycoparasite–host interaction represents a powerful strategy to obtain insight into the molecular events underlying these changes. The aim of this study is to identify genes whose expression is upregulated when the mycoparasite Stachybotrys elegans is in direct confrontation with its host Rhizoctonia solani. Suppression subtractive hybridization (SSH) was used to create a subtracted cDNA library, and differential screening was applied to identify the over-expressed transcripts. We report the analysis of 2,166 clones, among which 47% were upregulated during mycoparasitism. Two hundred and sixty-one clones were sequenced that corresponded to 94 unique genes. Forty-four of these were identified as novel genes, while the remainder showed similarity to a broad diversity of genes with putative functions related to toxin production, pathogenicity, and metabolism. As a result of mycoparasitism, 15 genes belonged to R. solani among which 9 genes were assigned putative functions. Quantitative RT-PCR was used to examine the upregulation of 12 genes during the course of mycoparasitism. Seven genes showed significant upregulation at least at one-time point during interaction of the mycoparasite with its host. This study describes a first step toward knowledge of S. elegans genome. The results present the useful application of EST analysis on S. elegans and provide preliminary indication of gene expression putatively involved in mycoparasitism.

Resistance to ivermectin and related drugs is an increasing problem for parasite control. The mechanism of ivermectin resistance in nematode parasites is currently unknown. Some P-glycoproteins and multidrug resistance proteins have been found to act as membrane transporters which pump drugs from the cell. A disruption of the mdr1a gene, which encodes a P-glycoprotein in mice, results in hypersensitivity to ivermectin. Genes encoding members of the P-glycoprotein family are known to exist in nematodes but the involvement of P-glycoprotein in nematode ivermectin-resistance has not been described. Our data suggest that a P-glycoprotein may play a role in ivermectin resistance in the sheep nematode parasite Haemonchus contortus. A full length P-glycoprotein cDNA from H. contortus has been cloned and sequenced. Analysis of the sequence showed 61–65% homology to other P-glycoprotein/multidrug resistant protein sequences, such as mice, human and Caenorhabditis elegans. Expression of P-glycoprotein mRNA was higher in ivermectin-selected than unselected strains of H. contortus. An alteration in the restriction pattern was also found for the genomic locus of P-glycoprotein derived from ivermectin-selected strains of H. contortus compared with unselected strains. P-glycoprotein gene structure and/or its transcription are altered in ivermectin-selected H. contortus. The multidrug resistance reversing agent, verapamil, increased the efficacy of ivermectin and moxidectin against a moxidectin-selected strain of this nematode in jirds (Meriones unguiculatus). These data indicate that a P-glycoprotein may be involved in resistance to ivermectin and other macrocyclic lactones in H. contortus.