Quantitative genetics of gastrointestinal strongyle burden and associated body condition in feral horses

This is the author accepted manuscript. The final version is available on open access from Elsevier via the DOI in this recordVariability in host resistance or tolerance to parasites is nearly ubiquitous, and is of key significance in understanding the evolutionary processes shaping host-parasite interactions. While ample research has been conducted on the genetics of parasite burden in livestock, relatively little has been done in free-living populations. Here, we investigate the sources of (co)variation in strongyle nematode faecal egg count (FEC) and body condition in Sable Island horses, a feral population in which parasite burden has previously been shown to negatively correlate with body condition. We used the quantitative genetic “animal model” to understand the sources of (co)variation in these traits, and tested for impacts of an important spatial gradient in habitat quality on the parameter estimates. Although FEC is significantly heritable (h 2 = 0.43 ± 0.11), there was no evidence for significant additive genetic variation in body condition (h 2 = 0.04 ± 0.07), and therefore there was also no significant genetic covariance between the two traits. The negative phenotypic covariance between these traits therefore does not derive principally from additive genetic effects. We also found that both FEC and body condition increase from east to west across the island, which indicates that the longitudinal environmental gradient is not responsible for the negative phenotypic association observed between these traits. There was also little evidence to suggest that quantitative genetic parameters were biased when an individual’s location along the island’s environmental gradient was not incorporated into the analysis. This research provides new and important insights into the genetic basis and adaptive potential of parasite resistance in free-living animals, and highlights the importance of environmental heterogeneity in modulating host-parasite interactions in wild vertebrate systems.Natural Sciences and Engineering Research Council of CanadaCanada Foundation for InnovationRoyal Society International ExchangeRyan Dubé and Heather Ryan Veterinary and Research FundNSERC-CREATE Host-Parasite Interactions Training ProgramLeverhulme Trus

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

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

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 significance of the Van Nuys prognostic index in the management of ductal carcinoma in situ

<p>Abstract</p> <p>Background</p> <p>Debate regarding the benefit of radiotherapy after local excision of ductal carcinoma <it>in situ </it>(DCIS) continues. The Van Nuys Prognostic Index (VNPI) is thought to be a useful aid in deciding which patients are at increased risk of local recurrence and who may benefit from adjuvant radiotherapy (RT). Recently published interim data from the Sloane project has showed that the VNPI score did significantly affect the chances of getting planned radiotherapy in the UK, suggesting that British clinicians may already be using this scoring system to assist in decision making. This paper independently assesses the prognostic validity of the VNPI in a British population.</p> <p>Patients and methods</p> <p>A retrospective review was conducted of all patients (n = 215) who underwent breast conserving surgery for DCIS at a single institution between 1997 – 2006. No patients included in the study received additional radiotherapy or hormonal treatment. Kaplan Meier survival curves were calculated, to determine disease free survival, for the total sample and a series of univariate analyses were performed to examine the value of various prognostic factors including the VNPI. The log-rank test was used to determine statistical significance of differential survival rates. Multivariate Cox regression analysis was performed to analyze the significance of the individual components of the VNPI. All analyses were conducted using SPSS software, version 14.5.</p> <p>Results</p> <p>The mean follow-up period was 53 months (range 12–97, SD19.9). Ninety five tumours were high grade (44%) and 84 tumours exhibited comedo necrosis (39%). The closest mean initial excision margin was 2.4 mm (range 0–22 mm, standard deviation 2.8) and a total of 72 tumours (33%) underwent further re-excision. The observed and the actuarial 8 year disease-free survival rates in this study were 91% and 83% respectively. The VNPI score and the presence of comedo necrosis were the only statistically significant prognostic indicators (P < 0.05).</p> <p>Conclusion</p> <p>This follow-up study of 215 patients with DCIS treated with local excision and observation alone is one of the largest series in which rates of recurrence are unaffected by radiation therapy, hormone manipulation or chemotherapy. It has afforded us the opportunity to assess the prognostic impact of patient and tumour characteristics free of any potentially confounding treatment related influences. The results suggest that the VNPI can be used to identify a subset of patients who are at risk of local recurrence and who may potentially benefit from RT.</p

Genomic-Bioinformatic Analysis of Transcripts Enriched in the Third-Stage Larva of the Parasitic Nematode Ascaris suum

Differential transcription in Ascaris suum was investigated using a genomic-bioinformatic approach. A cDNA archive enriched for molecules in the infective third-stage larva (L3) of A. suum was constructed by suppressive-subtractive hybridization (SSH), and a subset of cDNAs from 3075 clones subjected to microarray analysis using cDNA probes derived from RNA from different developmental stages of A. suum. The cDNAs (n = 498) shown by microarray analysis to be enriched in the L3 were sequenced and subjected to bioinformatic analyses using a semi-automated pipeline (ESTExplorer). Using gene ontology (GO), 235 of these molecules were assigned to ‘biological process’ (n = 68), ‘cellular component’ (n = 50), or ‘molecular function’ (n = 117). Of the 91 clusters assembled, 56 molecules (61.5%) had homologues/orthologues in the free-living nematodes Caenorhabditis elegans and C. briggsae and/or other organisms, whereas 35 (38.5%) had no significant similarity to any sequences available in current gene databases. Transcripts encoding protein kinases, protein phosphatases (and their precursors), and enolases were abundantly represented in the L3 of A. suum, as were molecules involved in cellular processes, such as ubiquitination and proteasome function, gene transcription, protein–protein interactions, and function. In silico analyses inferred the C. elegans orthologues/homologues (n = 50) to be involved in apoptosis and insulin signaling (2%), ATP synthesis (2%), carbon metabolism (6%), fatty acid biosynthesis (2%), gap junction (2%), glucose metabolism (6%), or porphyrin metabolism (2%), although 34 (68%) of them could not be mapped to a specific metabolic pathway. Small numbers of these 50 molecules were predicted to be secreted (10%), anchored (2%), and/or transmembrane (12%) proteins. Functionally, 17 (34%) of them were predicted to be associated with (non-wild-type) RNAi phenotypes in C. elegans, the majority being embryonic lethality (Emb) (13 types; 58.8%), larval arrest (Lva) (23.5%) and larval lethality (Lvl) (47%). A genetic interaction network was predicted for these 17 C. elegans orthologues, revealing highly significant interactions for nine molecules associated with embryonic and larval development (66.9%), information storage and processing (5.1%), cellular processing and signaling (15.2%), metabolism (6.1%), and unknown function (6.7%). The potential roles of these molecules in development are discussed in relation to the known roles of their homologues/orthologues in C. elegans and some other nematodes. The results of the present study provide a basis for future functional genomic studies to elucidate molecular aspects governing larval developmental processes in A. suum and/or the transition to parasitism

Population genetics of benzimidazole-resistant Haemonchus contortus and Haemonchus placei from buffalo and cattle: implications for the emergence and spread of resistance mutations

The population genetics of nematode parasites are poorly understood with practical reference to the selection and spread of anthelmintic resistance mutations. Haemonchus species are important to study the nematode population genetics due to their clinical importance in ruminant livestock, and the availability of genomic resources. In the present study, it has been examined that Haemonchus contortus and Haemonchus placei populations from three buffalo and nine cattle hosts. Seventy-three individual adult worms of H. contortus and 148 of H. placei were analysed using a panel of seven microsatellite markers. The number of alleles per locus in H. contortus and H. placei indicated that all populations were polymorphic for the microsatellites used in the present study. Genetic diversity parameters included high levels of allelic richness and heterozygosity, indicating effective population sizes, high mutation rates and high transmission frequencies in the area. Genetic structure parameters revealed low genetic differentiation between and high levels of genetic variation within H. contortus and H. placei populations. Population dynamic analyses showed an absence of heterozygosity excess in both species, suggesting that there was no deviation from genetic drift equilibrium. Our results provide a proof of concept for better understanding of the consequences of specific control strategies, climatic change or management strategies on the population genetics of anthelmintic resistance alleles in Haemonchus spp. infecting co-managed buffalo and cattle

Adaptation of gastrointestinal nematode parasites to host genotype: single locus simulation models

Background: Breeding livestock for improved resistance to disease is an increasingly important selection goal. However, the risk of pathogens adapting to livestock bred for improved disease resistance is difficult to quantify. Here, we explore the possibility of gastrointestinal worms adapting to sheep bred for low faecal worm egg count using computer simulation. Our model assumes sheep and worm genotypes interact at a single locus, such that the effect of an A allele in sheep is dependent on worm genotype, and the B allele in worms is favourable for parasitizing the A allele sheep but may increase mortality on pasture. We describe the requirements for adaptation and test if worm adaptation (1) is slowed by non-genetic features of worm infections and (2) can occur with little observable change in faecal worm egg count. Results: Adaptation in worms was found to be primarily influenced by overall worm fitness, viz. the balance between the advantage of the B allele during the parasitic stage in sheep and its disadvantage on pasture. Genetic variation at the interacting locus in worms could be from de novo or segregating mutations, but de novo mutations are rare and segregating mutations are likely constrained to have (near) neutral effects on worm fitness. Most other aspects of the worm infection we modelled did not affect the outcomes. However, the host-controlled mechanism to reduce faecal worm egg count by lowering worm fecundity reduced the selection pressure on worms to adapt compared to other mechanisms, such as increasing worm mortality. Temporal changes in worm egg count were unreliable for detecting adaptation, despite the steady environment assumed in the simulations. Conclusions: Adaptation of worms to sheep selected for low faecal worm egg count requires an allele segregating in worms that is favourable in animals with improved resistance but less favourable in other animals. Obtaining alleles with this specific property seems unlikely. With support from experimental data, we conclude that selection for low faecal worm egg count should be stable over a short time frame (e.g. 20 years). We are further exploring model outcomes with multiple loci and comparing outcomes to other control strategies

Genome-wide analysis of ivermectin response by Onchocerca volvulus reveals that genetic drift and soft selective sweeps contribute to loss of drug sensitivity

Treatment of onchocerciasis using mass ivermectin administration has reduced morbidity and transmission throughout Africa and Central/South America. Mass drug administration is likely to exert selection pressure on parasites, and phenotypic and genetic changes in several Onchocerca volvulus populations from Cameroon and Ghana-exposed to more than a decade of regular ivermectin treatment-have raised concern that sub-optimal responses to ivermectin's anti-fecundity effect are becoming more frequent and may spread.Pooled next generation sequencing (Pool-seq) was used to characterise genetic diversity within and between 108 adult female worms differing in ivermectin treatment history and response. Genome-wide analyses revealed genetic variation that significantly differentiated good responder (GR) and sub-optimal responder (SOR) parasites. These variants were not randomly distributed but clustered in ~31 quantitative trait loci (QTLs), with little overlap in putative QTL position and gene content between the two countries. Published candidate ivermectin SOR genes were largely absent in these regions; QTLs differentiating GR and SOR worms were enriched for genes in molecular pathways associated with neurotransmission, development, and stress responses. Finally, single worm genotyping demonstrated that geographic isolation and genetic change over time (in the presence of drug exposure) had a significantly greater role in shaping genetic diversity than the evolution of SOR.This study is one of the first genome-wide association analyses in a parasitic nematode, and provides insight into the genomics of ivermectin response and population structure of O. volvulus. We argue that ivermectin response is a polygenically-determined quantitative trait (QT) whereby identical or related molecular pathways but not necessarily individual genes are likely to determine the extent of ivermectin response in different parasite populations. Furthermore, we propose that genetic drift rather than genetic selection of SOR is the underlying driver of population differentiation, which has significant implications for the emergence and potential spread of SOR within and between these parasite populations

Peptidases compartmentalized to the Ascaris suum intestinal lumen and apical intestinal membrane

The nematode intestine is a tissue of interest for developing new methods of therapy and control of parasitic nematodes. However, biological details of intestinal cell functions remain obscure, as do the proteins and molecular functions located on the apical intestinal membrane (AIM), and within the intestinal lumen (IL) of nematodes. Accordingly, methods were developed to gain a comprehensive identification of peptidases that function in the intestinal tract of adult female Ascaris suum. Peptidase activity was detected in multiple fractions of the A. suum intestine under pH conditions ranging from 5.0 to 8.0. Peptidase class inhibitors were used to characterize these activities. The fractions included whole lysates, membrane enriched fractions, and physiological- and 4 molar urea-perfusates of the intestinal lumen. Concanavalin A (ConA) was confirmed to bind to the AIM, and intestinal proteins affinity isolated on ConA-beads were compared to proteins from membrane and perfusate fractions by mass spectrometry. Twenty-nine predicted peptidases were identified including aspartic, cysteine, and serine peptidases, and an unexpectedly high number (16) of metallopeptidases. Many of these proteins co-localized to multiple fractions, providing independent support for localization to specific intestinal compartments, including the IL and AIM. This unique perfusion model produced the most comprehensive view of likely digestive peptidases that function in these intestinal compartments of A. suum, or any nematode. This model offers a means to directly determine functions of these proteins in the A. suum intestine and, more generally, deduce the wide array functions that exist in these cellular compartments of the nematode intestine