A 4 year observation of gastrointestinal nematode egg counts,nemabiomes and the benzimidazole resistance genotypes of Teladorsagia circumcincta on a Scottish sheep farm

Anthelmintic resistance threatens the sustainability of sheep production globally. Advice regarding strategies to reduce the development of anthelmintic resistance incorporates the outcomes of modelling exercises. Further understanding of gastrointestinal nematode species diversity, and population dynamics and genetics (which may vary between species) is required to refine these models; and field studies combining faecal egg outputs, species composition and resistance genetics are needed to calibrate them. In this study, faecal samples were taken from ewes and lambs on a commercial farm in south-eastern Scotland at approximately 3 t-4 week intervals between spring and autumn over a period of 4 years. Faecal egg counts were performed on these samples, and L₃ were collected from pooled coprocultures. Deep amplicon sequencing was used to determine both the species composition of these L₃ and the proportions of benzimidazole-resistant single nucleotide polymorphisms in the isotype-1 β-tubulin locus of the predominant species, Teladorsagia circumcincta L₃. Despite consistent management throughout the study, the results show variation in gastrointestinal nematode species composition with time and between age groups, that was potentially associated with weather conditions. The F200Y benzimidazole resistance mutation is close to genetic fixation in the T. circumcincta population on this farm. There was no evidence of variation in isotype-1 β-tubulin single nucleotide polymorphisms frequency between age groups, and no genetic evidence of reversion to benzimidazole susceptibility, despite targeted benzimidazole usage. This study highlights the need to include speciation when investigating gastrointestinal nematode epidemiology and anthelmintic resistance, and serves as an example of how genetic data may be analysed alongside species diversity and faecal egg counts, when markers for other anthelmintic classes are identified.     

Phylogenetic Characterization of β-Tubulins and Development of Pyrosequencing Assays for Benzimidazole Resistance in Cattle Nematodes

Control of helminth infections is a major task in livestock production to prevent health constraints and economic losses. However, resistance to established anthelmintic substances already impedes effective anthelmintic treatment in many regions worldwide. Thus, there is an obvious need for sensitive and reliable methods to assess the resistance status of at least the most important nematode populations. Several single nucleotide polymorphisms (SNPs) in the β-tubulin isotype 1 gene of various nematodes correlate with resistance to benzimidazoles (BZ), a major anthelmintic class. Here we describe the full-length β-tubulin isotype 1 and 2 and α-tubulin coding sequences of the cattle nematode Ostertagia ostertagi. Additionally, the Cooperia oncophora α-tubulin coding sequence was identified. Phylogenetic maximum-likelihood analysis revealed that both isotype 1 and 2 are orthologs to the Caenorhabditis elegans ben-1 gene which is also associated with BZ resistance upon mutation. In contrast, a Trichuris trichiura cDNA, postulated to be β-tubulin isotype 1 involved in BZ resistance in this human parasite, turned out to be closely related to C. elegans β-tubulins tbb-4 and mec-7 and would therefore represent the first non-ben-1-like β-tubulin to be under selection through treatment with BZs. A pyrosequencing assay was established to detect BZ resistance associated SNPs in β-tubulin isotype 1 codons 167, 198 and 200 of C. oncophora and O. ostertagi. PCR-fragments representing either of the two alleles were combined in defined ratios to evaluate the pyrosequencing assay. The correlation between the given and the measured allele frequencies of the respective SNPs was very high. Subsequently laboratory isolates and field populations with known resistance status were analyzed. With the exception of codon 167 in Cooperia, increases of resistance associated alleles were detected for all codons in at least one of the phenotypically resistant population. Pyrosequencing provides a fast, inexpensive and sensitive alternative to conventional resistance detection methods.     

The Emergence of Resistance to the Benzimidazole Anthlemintics in Parasitic Nematodes of Livestock Is Characterised by Multiple Independent Hard and Soft Selective Sweeps

Anthelmintic resistance is a major problem for the control of parasitic nematodes of livestock and of growing concern for human parasite control. However, there is little understanding of how resistance arises and spreads or of the “genetic signature” of selection for this group of important pathogens. We have investigated these questions in the system for which anthelmintic resistance is most advanced; benzimidazole resistance in the sheep parasites Haemonchus contortus and Teladorsagia circumcincta. Population genetic analysis with neutral microsatellite markers reveals that T. circumcincta has higher genetic diversity but lower genetic differentiation between farms than H. contortus in the UK. We propose that this is due to epidemiological differences between the two parasites resulting in greater seasonal bottlenecking of H. contortus. There is a remarkably high level of resistance haplotype diversity in both parasites compared with drug resistance studies in other eukaryotic systems. Our analysis suggests a minimum of four independent origins of resistance mutations on just seven farms for H. contortus, and even more for T. circumincta. Both hard and soft selective sweeps have occurred with striking differences between individual farms. The sweeps are generally softer for T. circumcincta than H. contortus, consistent with its higher level of genetic diversity and consequent greater availability of new mutations. We propose a model in which multiple independent resistance mutations recurrently arise and spread by migration to explain the widespread occurrence of resistance in these parasites. Finally, in spite of the complex haplotypic diversity, we show that selection can be detected at the target locus using simple measures of genetic diversity and departures from neutrality. This work has important implications for the application of genome-wide approaches to identify new anthelmintic resistance loci and the likelihood of anthelmintic resistance emerging as selection pressure is increased in human soil-transmitted nematodes by community wide treatment programs.     

Gene expression changes in a P-glycoprotein (Tci-pgp-9) putatively associated with ivermectin resistance in Teladorsagia circumcincta

Anthelmintic resistance in parasitic nematodes of small ruminants is widespread and, in some parts of the world, threatens the sustainability of sheep production. The genetic changes underlying resistance to anthelmintics, particularly ivermectin (IVM), remain to be determined. The majority of studies to date have investigated target site mutations; relatively little attention has been paid to the role of changes in gene expression. In this study, we investigated the expression of putative drug transporter molecules, P-glycoproteins (Pgps), in Teladorsagia circumcincta, the predominant parasitic nematode species of sheep in the UK and the major anthelmintic resistant species. Utilising a degenerate PCR approach, 11 partial Pgp sequences were identified. Constitutive differences in gene expression between an IVM-susceptible (MTci2) and a multidrug-resistant (MTci5) isolate were determined for 10 of the Pgps using the ΔΔCt TaqMan® real-time PCR method. Gene expression differences were particularly marked in one of these genes, namely Tci-pgp-9. In the MTci5 isolate, statistically significant increases in Tci-pgp-9 expression, at the mRNA level, were observed across all life-cycle stages and most notably in eggs (55-fold increase). Comparison of the partial Tci-pgp-9 nucleotide sequences from MTci2 and MTci5 also identified high levels of polymorphism. This work has shown that constitutively increased expression in Tci-pgp-9, coupled with increased sequence polymorphism, could play a role in allowing multidrug-resistant T. circumcincta to survive IVM exposure. The genetic changes underpinning these gene expression changes remain to be elucidated and need to be investigated in other isolates. These changes could form the basis of an IVM resistance marker to monitor the spread of resistance and to evaluate management practices aimed at delaying its spread.     

Identification of the amino acids in the Major Histocompatibility Complex class II region of Scottish Blackface sheep that are associated with resistance to nematode infection

Lambs with the Major Histocompatibility Complex DRB1*1101 allele have been shown to produce fewer nematode eggs following natural and deliberate infection. These sheep also possess fewer adult Teladorsagia circumcincta than sheep with alternative alleles at the DRB1 locus. However, it is unclear if this allele is responsible for the reduced egg counts or merely acts as a marker for a linked gene. This study defined the MHC haplotypes in a population of naturally infected Scottish Blackface sheep by PCR amplification and sequencing, and examined the associations between MHC haplotypes and faecal egg counts by generalised linear mixed modelling. The DRB1*1101 allele occurred predominately on one haplotype and a comparison of haplotypes indicated that the causal mutation or mutations occurred in or around this locus. Additional comparisons with another resistant haplotype indicated that mutations in or around the DQB2*GU191460 allele were also responsible for resistance to nematode infections. Further analyses identified six amino acid substitutions in the antigen binding site of DRB1*1101 that were significantly associated with reductions in the numbers of adult T. circumcincta.     

A novel ABCG-like transporter of Trypanosoma cruzi is involved in natural resistance to benznidazole

Benznidazole (BZ) is one of the two drugs used for Chagas disease treatment. Nevertheless therapeutic failures of BZ have been reported, which were mostly attributed to variable drug susceptibility among Trypanosoma cruzi strains. ATP-binding cassette (ABC) transporters are involved in a variety of translocation processes and some members have been implicated in drug resistance. Here we report the characterisation of the first T. cruzi ABCG transporter gene, named TcABCG1, which is over-expressed in parasite strains naturally resistant to BZ. Comparison of TcABCG1 gene sequence of two TcI BZ-resistant strains with CL Brener BZ-susceptible strain showed several single nucleotide polymorphisms, which determined 11 amino acid changes. CL Brener transfected with TcI transporter genes showed 40-47% increased resistance to BZ, whereas no statistical significant increment in drug resistance was observed when CL Brener was transfected with the homologous gene. Only in the parasites transfected with TcI genes there was 2-2.6-fold increased abundance of TcABCG1 transporter protein. The analysis in wild type strains also suggests that the level of TcABCG1 transporter is related to BZ natural resistance. The characteristics of untranslated regions of TcABCG1 genes of BZ-susceptible and resistant strains were investigated by computational tools.     

Beta-tubulin complementary DNA sequence variations observed between cyathostomins from benzimidazole-susceptible and -resistant populations

The molecular mechanism of benzimidazole (BZ) resistance in cyathostomins of horses is still unclear. Previous studies revealed that the TTC or TAC polymorphism in codon 200 of the beta-tubulin isotype 1 gene is not as strictly correlated with BZ resistance as in trichostrongyles in sheep. To identify further sites of polymorphism within the beta-tubulin gene related to BZ resistance, complete complementary DNAs (cDNAs) encoding beta-tubulin of adult worms of Cylicocyclus nassatus, Cyathostomum pateratum, Cyathostomum coronatum, Cyathostomum catinatum, Cylicostephanus longibursatus, and Cylicostephanus goldi of a BZ-resistant cyathostomin population were characterized using specific primers. The cDNA sequence of each species spans 1,429 bp, encoding a protein of 448 amino acids. The interspecific identities are 95.2-99.6% at the nucleotide and 98.7-100.0% at the peptide level. The comparison of the amino acid sequences of individuals isolated from the BZ-resistant cyathostomin population with those from individuals of Cc. nassatus, Cy. coronatum, Cy. pateratum, and Cy. catinatum of a BZ-susceptible one showed differing amino acids in 11 positions. The commonness of a phenylalanine to tyrosine mutation at position 167 in all the 6 cyathostomin species isolated from a BZ-resistant population suggests its involvement in the molecular mechanism in BZ resistance.     

A tandemly repetitive DNA sequence is present at diverse locations in the genome of Ostertagia circumcincta

A novel repetitive DNA sequence in the sheep parasitic nematode Ostertagia circumcincta was cloned and sequenced. This 1.2-kb sequence (Oc1B) was not found in the closely related cattle parasite Ostertagia ostertagi, nor in the more distantly related sheep parasites Haemonchus contortus or Trichostrongylus colubriformis. Sequences similar to Oc1B were found at various genomic locations and contained a pair of 33-bp direct repeats. Oc1B also contained a single copy of a 218-bp sequence (designated OcREP) which was present in 100 to 200 copies in the O. circumcincta genome and mostly organized in distinctive tandem arrays. The dual organizational pattern of OcREP as both a satellite-like sequence as well as interspersed as single copies amongst dissimilar sequences adds to the growing evidence for the fluidity of the parasitic nematode genome, and of eukaryotic genomes in general.     

Phenotypic analysis of host-parasite interactions in lambs infected with Teladorsagia circumcincta

Female Blackface lambs expected to exhibit genetic variability for resistance to gastrointestinal nematodes, were either exposed to continuous experimental infections of Teladorsagia circumcincta or were sham-dosed to monitor phenotypic responses to infection. As a measure of parasitism and host response, worm-eggs in faeces (faecal egg count, FEC) were counted over a 3-month period and worm burdens were ascertained at post-mortem. The host response to the infection was also measured by differential counts of white blood cells, anti-T. circumcincta IgA antibody levels and body weight. Results suggest that nematode abundance (mean number of parasites per host) and prevalence (proportion of infected animals) were maximal shortly after the beginning of infection (21 days p.i.) when virtually all the infected animals were shedding worm eggs. Increasing anti-T. circumcincta IgA antibody and eosinophil concentrations were associated with a reduction in total numbers of adult worms and an increase in the frequency of early L4s. The data also suggest that genetic selection for an enhanced anti-T. circumcincta IgA response might complement selection based on a reduced FEC as a strategy to select for resistance to gastrointestinal nematodes.     

Haemonchus contortus and Trichostrongylus colubriformis did not adapt to long-term exposure to sheep that were genetically resistant or susceptible to nematode infections

We tested the hypothesis that Haemonchus contortus and Trichostrongylus colubriformis would adapt to long-term exposure to sheep that were either genetically resistant or susceptible to H. contortus. Sheep genotypes were from lines with 10 years prior selection for low (resistant, R) or high (susceptible, S) faecal worm egg count (WEC) following H. contortus infection. Long-term exposure of H. contortus and T.colubriformis to R or S genotypes was achieved using serial passage for up to 30 nematode generations. Thus, we generated four nematode strains; one strain of each species solely exposed to R sheep and one strain of each species solely exposed to S sheep. Considerable host genotype differences in mean WEC during serial passage confirmed adequate nematode selection pressure for both H. contortus (R 4900 eggs per gram (epg), S 19,900 epg) and T. colubriformis (R 5300 epg, S 13,500 epg). Adaptation of nematode strain to host genotype was tested using seven cross-classified tests for H. contortus, and two cross-classified and one outbred genotype test for T. colubriformis. In the cross-classified design, where each strain infects groups of R, S or randomly bred control sheep, parasite adaptation would be indicated by a significant host genotype by nematode strain interaction for traits indicating parasite reproductive success; specifically WEC and, for H. contortus strains, packed cell volume. We found no significant evidence of parasite adaptation to host genotype (P > 0.05) for either the H. contortus or T. colubriformis strains. Therefore, we argue that nematodes will not adapt quickly to sheep bred for nematode resistance, where selection is based on low WEC, although selecting sheep using a subset of immune functions may increase adaptation risk. Our results support the hypothesis that nematode resistance is determined by many genes each with relatively small effect. In conclusion, selection of sheep for nematode resistance using WEC should be sustainable in the medium to long-term.     

Contrasting genetic structures of two parasitic nematodes,determined on the basis of neutral microsatellite markers and selected anthelmintic resistance markers

For the first time, the neutral genetic relatedness of natural populations of Trichostrongylid nematodes was investigated in relation to polymorphism of the β‐tubulin gene, which is selected for anthelminthic treatments. The aim of the study was to assess the contribution of several evolutionary processes: migration and genetic drift by neutral genetic markers and selection by anthelminthic treatments on the presence of resistance alleles at β‐tubulin. We studied two nematode species (Teladorsagia circumcincta and Haemonchus contortus) common in temperate climatic zones; these species are characterized by contrasting life history traits. We studied 10 isolated populations of goat nematode parasites: no infected adult goat had been exchanged after the herds were established. Beta‐tubulin polymorphism was similar in these two species. One and two β‐tubulin alleles from T. circumcincta and H. contortus respectively were shared by several populations. Most of the β‐tubulin alleles were ‘private’ alleles. No recombination between alleles was detected in BZ‐resistant alleles from T. circumcincta and H. contortus. The T. circumcincta populations have not diverged much since their isolation (F_ST <0.08), whereas H. contortus displayed marked local genetic differentiation (F_ST ranging from 0.08 to 0.18). These findings suggest that there are severe bottlenecks in the H. contortus populations, possibly because of their reduced abundance during unfavourable periods and their high reproductive rate, which allows the species to persist even after severe population reduction. Overall, the data reported contradict the hypothesis of the origin of β‐tubulin resistance alleles in these populations from a single mutational event, but two other hypotheses (recurrent mutation generating new alleles in isolated populations and the introduction of existing alleles) emerge as equally likely.     

A combined microscopic-molecular method for the diagnosis of strongylid infections in sheep

We evaluated a combined microscopic-molecular approach for the diagnosis of key strongylid infections in sheep using panels of well-defined control and test samples. The method established is based on the separation of nematode eggs from faecal samples using a salt flotation procedure, the extraction and column-purification of genomic DNA, followed by real-time PCR and melting-curve analysis. Specific and semi-quantitative amplification from (a minimum of 0.1-2.0 pg) genomic DNA of Haemonchus contortus, Teladorsagia circumcincta, Trichostrongylus spp., Cooperia oncophora, Oesophagostomum columbianum, Oesophagostomum venulosum or Chabertia ovina is achieved using a specific, forward oligonucleotide primer located in the second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA (rDNA) together with a conserved reverse primer in the large subunit of rDNA. Using a panel of well-defined genomic DNA samples from eggs from sheep monospecifically infected with H. contortus or Te. circumcincta, there was a correlation between cycle threshold (Ct) values in the PCR and numbers of egg per gram of faeces, thus allowing the semi-quantitation of parasite DNA in faeces. The findings of the present study indicate that a microscopic-molecular approach provides a useful tool for diagnosis, for epidemiological and ecological surveys as well as for integration into parasite monitoring, drug resistance (i.e. ‘egg count reduction’) testing or control programmes, particularly following semi- or full-automation.     

Population genetic structure and history of a generalist parasite infecting multiple sympatric host species

Host specificity is predicted to shape patterns of parasite gene flow between host species; specialist parasites should have low gene flow between host species, while generalists are predicted to have high gene flow between species. However, even for generalist parasites external forces, including ecological differences between host species may sometimes intervene to limit gene flow and create genetic structure. To investigate the potential for cryptic parasite genetic structure to arise under such circumstances, we examined the population genetic structure and history of the generalist nematode, Trichostrongylus axei, infecting six sympatric wild ungulate species in North America. Using genotypes for 186 T. axei larvae at two mitochondrial genes, cox1 and nad4, we found that T. axei was completely panmictic across host species, with 0% of genetic variation structured between host species and 97% within individual hosts. In addition, T. axei showed no evidence of recent genetic bottlenecks, had high nucleotide diversities (above 2%), and an effective population size estimated to be in the tens of millions. Our result that T. axei maintains high rates of gene flow between multiple sympatric host species adds to a growing body of information on trichostrongylid population genetic structure in different ecological contexts. Furthermore, the high rates of gene flow, coupled with high levels of genetic diversity and large effective population size which we observed in T. axei, point to a potentially broad capacity for rapid evolutionary change in this parasite.     

Genetic Variation of Taenia Pisiformis Collected from Sichuan,China, Based on the Mitochondrial Cytochrome b gene

Taenia pisiformis is one of the most important parasites of canines and rabbits. T. pisiformis cysticercus (the larval stage) causes severe damage to rabbit breeding, which results in huge economic losses. In this study, the genetic variation of T. pisiformis was determined in Sichuan Province, China. Fragments of the mitochondrial cytochrome b (cytb) (922 bp) gene were amplified in 53 isolates from 8 regions of T. pisiformis. Overall, 12 haplotypes were found in these 53 cytb sequences. Molecular genetic variations showed 98.4% genetic variation derived from intra-region. F_ST and Nm values suggested that 53 isolates were not genetically differentiated and had low levels of genetic diversity. Neutrality indices of the cytb sequences showed the evolution of T. pisiformis followed a neutral mode. Phylogenetic analysis revealed no correlation between phylogeny and geographic distribution. These findings indicate that 53 isolates of T. pisiformis keep a low genetic variation, which provide useful knowledge for monitoring changes in parasite populations for future control strategies.     

Expression of genes in gastrointestinal and lymphatic tissues during parasite infection in sheep genetically resistant or susceptible to Trichostrongylus colubriformis and Haemonchus contortus

Resistance to an acute gastrointestinal nematode (GIN) infection is dependent on the ability of the host to recognise the parasite and mount a protective Th2 response. It is hypothesised that lambs which are genetically susceptible to GIN will differentially up-regulate Th1-type genes and therefore remain susceptible to chronic parasitism compared with genetically resistant lambs which will differentially up-regulate Th2-type genes and clear the parasite infection. Two selection flocks, in which lines of Merino sheep produced lambs genetically resistant or susceptible to GIN, were acutely challenged once or thrice with either Haemonchus contortus or Trichostrongylus colubriformis. Faecal-egg counts (FECs), and plasma and tissue anti-parasite (H. contortus or T. colubriformis) antibody isotype responses showed that resistant animals challenged three times with T. colubriformis established a protective Th2 response (negligible FEC, IgG1 and IgE) whereas susceptible animals required multiple challenges to establish a significant IgG1 response despite FECs remaining high. Trichostrongylus colubriformis elicited a more pronounced host response than H. contortus. RNA extracted from tissues at the site of each parasite infection and associated lymph nodes were interrogated by microarray and quantitative PCR analyses to correlate host gene expression to FECs and antibody responses. The IFN-γ inducible gene cxcl10 was up-regulated in the susceptible line of the Trichostrongylus selection flock sheep after a tertiary challenge with the parasites H. contortus and T. colubriformis. However, a uniform pattern of genes was not up-regulated in resistant animals from both selection flocks during both parasite infections, suggesting that the mode of host resistance to these parasites is different, although some similarities in host susceptibility were apparent.     

Global spread and genetic variants of the two CYP9M10 haplotype forms associated with insecticide resistance in Culex quinquefasciatus Say

Insecticide resistance develops as a genetic factor (allele) conferring lower susceptibility to insecticides proliferates within a target insect population under strong positive selection. Intriguingly, a resistance allele pre-existing in a population often bears a series of further adaptive allelic variants through new mutations. This phenomenon occasionally results in replacement of the predominating resistance allele by fitter new derivatives, and consequently, development of greater resistance at the population level. The overexpression of the cytochrome P450 gene CYP9M10 is associated with pyrethroid resistance in the southern house mosquito Culex quinquefasciatus. Previously, we have found two genealogically related overexpressing CYP9M10 haplotypes, which differ in gene copy number (duplicated and non-duplicated). The duplicated haplotype was derived from the non-duplicated overproducer probably recently. In the present study, we investigated allelic series of CYP9M10 involved in three C. quinquefasciatus laboratory colonies recently collected from three different localities. Duplicated and non-duplicated overproducing haplotypes coexisted in African and Asian colonies indicating a global distribution of both haplotype lineages. The duplicated haplotypes both in the Asian and African colonies were associated with higher expression levels and stronger resistance than non-duplicated overproducing haplotypes. There were slight variation in expression level among the non-duplicated overproducing haplotypes. The nucleotide sequences in coding and upstream regions among members of this group also showed a little diversity. Non-duplicated overproducing haplotypes with relatively higher expression were genealogically closer to the duplicated haplotypes than the other non-duplicated overproducing haplotypes, suggesting multiple cis-acting mutations before duplication.     

Genome-wide scan reveals important additive and non-additive genetic effects associated with resistance to Haemonchus contortus in Florida Native sheep

Florida Native sheep is among the sheep breeds best adapted to humid and hot climatic conditions such as those of Florida, USA, and have shown a superior ability to regulate nematode burdens. This is one of the oldest sheep breeds in North America and is an endangered species. To ensure genetic diversity and long-term survival of the breed, protection of the current genetic stock is critical and conservation efforts are required to promote its breeding and production. The objective of the present study was to investigate the importance of additive and non-additive genetic effects on resistance to natural Haemonchus contortus infections in Florida Native sheep using a whole genome scan. A total of 200 sheep were evaluated in the present study. Phenotypic records included faecal egg count (FEC, eggs/gram), FAMACHA® score, packed cell volume (PCV, %), body condition score and average daily gain (ADG, kg). Sheep were genotyped using the GGP Ovine 50K SNP chip and 45.2 k single nucleotide polymorphism (SNP) markers spanning the entire genome were available for quality control procedures. Mixed models were used to analyse the response variables and included the identity by state matrix to control for population structure. Bonferroni correction was used to control for multiple testing and a second arbitrary threshold (0.1 × 10^?3) was used. Fifteen SNPs with additive and non-additive genetic effects and located in Ovis aries chromosome OAR1, 2, 3, 6, 8, 10, 11, 12, 13 and 21 were associated with FEC, FAMACHA® score, PCV and ADG. These SNPs could be potential genetic markers for resistance to natural H. contortus exposure in Florida Native sheep.     

Quantitative trait loci for resistance to trichostrongylid infection in Spanish Churra sheep

Background

For ruminants reared on grazing systems, gastrointestinal nematode (GIN) parasite infections represent the class of diseases with the greatest impact on animal health and productivity. Among the many possible strategies for controlling GIN infection, the enhancement of host resistance through the selection of resistant animals has been suggested by many authors. Because of the difficulty of routinely collecting phenotypic indicators of parasite resistance, information derived from molecular markers may be used to improve the efficiency of classical genetic breeding.

Methods

A total of 181 microsatellite markers evenly distributed along the 26 sheep autosomes were used in a genome scan analysis performed in a commercial population of Spanish Churra sheep to detect chromosomal regions associated with parasite resistance. Following a daughter design, we analysed 322 ewes distributed in eight half-sib families. The phenotypes studied included two faecal egg counts (LFEC0 and LFEC1), anti-Teladorsagia circumcincta LIV IgA levels (IgA) and serum pepsinogen levels (Peps).

Results

The regression analysis revealed one QTL at the 5% genome-wise significance level on chromosome 6 for LFEC1 within the marker interval BM4621-CSN3. This QTL was found to be segregating in three out of the eight families analysed. Four other QTL were identified at the 5% chromosome-wise level on chromosomes 1, 10 and 14. Three of these QTL influenced faecal egg count, and the other one had an effect on IgA levels.

Conclusion

This study has successfully identified segregating QTL for parasite resistance traits in a commercial population. For some of the QTL detected, we have identified interesting coincidences with QTL previously reported in sheep, although most of those studies have been focused on young animals. Some of these coincidences might indicate that some common underlying loci affect parasite resistance traits in different sheep breeds. The identification of new QTL may suggest the existence of complex host-parasite relationships that have unique features depending on the host-parasite combination, perhaps due to the different mechanisms underlying resistance in adult sheep (hypersensitivity reactions) and lambs (immunity). The most significant QTL identified on chromosome 6 for LFEC₁ may be the target for future fine-mapping research efforts.     

Genetic Diversity of Taenia asiatica from Thailand and Other Geographical Locations as Revealed by Cytochrome c Oxidase Subunit 1 Sequences

Twelve 924 bp cytochrome c oxidase subunit 1 (cox1) mitochondrial DNA sequences from Taenia asiatica isolates from Thailand were aligned and compared with multiple sequence isolates from Thailand and 6 other countries from the GenBank database. The genetic divergence of T. asiatica was also compared with Taenia saginata database sequences from 6 different countries in Asia, including Thailand, and 3 countries from other continents. The results showed that there were minor genetic variations within T. asiatica species, while high intraspecies variation was found in T. saginata. There were only 2 haplotypes and 1 polymorphic site found in T. asiatica, but 8 haplotypes and 9 polymorphic sites in T. saginata. Haplotype diversity was very low, 0.067, in T. asiatica and high, 0.700, in T. saginata. The very low genetic diversity suggested that T. asiatica may be at a risk due to the loss of potential adaptive alleles, resulting in reduced viability and decreased responses to environmental changes, which may endanger the species.