Comparative evaluation of different molecular methods for DNA extraction from individual Teladorsagia circumcincta nematodes.

BACKGROUND: The purpose of this study was to develop a reliable DNA extraction protocol to use on individual Teladorsagia circumcincta nematode specimens to produce high quality DNA for genome sequencing and phylogenetic analysis. Pooled samples have been critical in providing the groundwork for T. circumcincta genome construction, but there is currently no standard method for extracting high-quality DNA from individual nematodes. 11 extraction kits were compared based on DNA quality, yield, and processing time. RESULTS: 11 extraction protocols were compared, and the concentration and purity of the extracted DNA was quantified. Median DNA concentration among all methods measured on NanoDrop 2000™ ranged between 0.45-11.5 ng/µL, and on Qubit™ ranged between undetectable - 0.962 ng/µL. Median A260/280 ranged between 0.505-3.925, and median A260/230 ranged - 0.005 - 1.545. Larval exsheathment to remove the nematode cuticle negatively impacted DNA concentration and purity. CONCLUSIONS: A Schistosoma sp. DNA extraction method was determined as most suitable for individual T. circumcincta nematode specimens due to its resulting DNA concentration, purity, and relatively fast processing time.

Bioinformatic analysis of eosinophil activity and its implications for model and target species.

Eosinophils are important immune cells that have been implicated in resistance to gastrointestinal nematode (GIN) infections in both naturally and experimentally infected sheep. Proteins of particular importance appear to be IgA-Fc alpha receptor (FcαRI), C-C chemokine receptor type 3 (CCR3), proteoglycan 3 (PRG3, major basic protein 2) and EPX (eosinophil peroxidase). We used known human nucleotide sequences to search the ruminant genomes, followed by translation to protein and sequence alignments to visualize differences between sequences and species. Where a sequence was retrieved for cow, but not for sheep and goat, this was used additionally as a reference sequence. In this review, we show that eosinophil function varies among host species. Consequently, investigations into the mechanisms of ruminant immune responses to GIN should be conducted using the natural host. Specifically, we address differences in protein sequence and structure for eosinophil proteins.

Analysis of pooled genome sequences from Djallonke and Sahelian sheep of Ghana reveals co-localisation of regions of reduced heterozygosity with candidate genes for disease resistance and adaptation to a tropical environment.

BACKGROUND: The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually. RESULTS: Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis. CONCLUSIONS: Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases.

Implementation of an extended ZINB model in the study of low levels of natural gastrointestinal nematode infections in adult sheep.

BACKGROUND: In this study, two traits related with resistance to gastrointestinal nematodes (GIN) were measured in 529 adult sheep: faecal egg count (FEC) and activity of immunoglobulin A in plasma (IgA). In dry years, FEC can be very low in semi-extensive systems, such as the one studied here, which makes identifying animals that are resistant or susceptible to infection a difficult task. A zero inflated negative binomial model (ZINB) model was used to calculate the extent of zero inflation for FEC; the model was extended to include information from the IgA responses. RESULTS: In this dataset, 64% of animals had zero FEC while the ZINB model suggested that 38% of sheep had not been recently infected with GIN. Therefore 26% of sheep were predicted to be infected animals with egg counts that were zero or below the detection limit and likely to be relatively resistant to nematode infection. IgA activities of all animals were then used to decide which of the sheep with zero egg counts had been exposed and which sheep had not been recently exposed. Animals with zero FEC and high IgA activity were considered resistant while animals with zero FEC and low IgA activity were considered as not recently infected. For the animals considered as exposed to the infection, the correlations among the studied traits were estimated, and the influence of these traits on the discrimination between unexposed and infected animals was assessed. CONCLUSIONS: The model presented here improved the detection of infected animals with zero FEC. The correlations calculated here will be useful in the development of a reliable index of GIN resistance that could be of assistance for the study of host resistance in studies based on natural infection, especially in adult sheep, and also the design of breeding programs aimed at increasing resistance to parasites.

A comprehensive mapping of the structure and gene organisation in the sheep MHC class I region.

BACKGROUND: The major histocompatibility complex (MHC) is a chromosomal region that regulates immune responsiveness in vertebrates. This region is one of the most important for disease resistance because it has been associated with resistance or susceptibility to a wide variety of diseases and because the MHC often accounts for more of the variance than other loci. Selective breeding for disease resistance is becoming increasingly common in livestock industries, and it is important to determine how this will influence MHC polymorphism and resistance to diseases that are not targeted for selection. However, in sheep the order and sequence of the protein coding genes is controversial. Yet this information is needed to determine precisely how the MHC influences resistance and susceptibility to disease. METHODS: CHORI bacterial artificial chromosomes (BACs) known to contain sequences from the sheep MHC class I region were sub-cloned, and the clones partially sequenced. The resulting sequences were analysed and re-assembled to identify gene content and organisation within each BAC. The low resolution MHC class I physical map was then compared to the cattle reference genome, the Chinese Merino sheep MHC map published by Gao, et al. (2010) and the recently available sheep reference genome. RESULTS: Immune related class I genes are clustered into 3 blocks; beta, kappa and a novel block not previously identified in other organisms. The revised map is more similar to Bovidae maps than the previous sheep maps and also includes several genes previously not annotated in the Chinese Merino BAC assembly and others not currently annotated in the sheep reference chromosome 20. In particular, the organisation of nonclassical MHC class I genes is similar to that present in the cattle MHC. Sequence analysis and prediction of amino acid sequences of MHC class I classical and nonclassical genes was performed and it was observed that the map contained one classical and eight nonclassical genes together with three possible pseudogenes. CONCLUSIONS: The comprehensive physical map of the sheep MHC class I region enhances our understanding of the genetic architecture of the class I MHC region in sheep and will facilitate future studies of MHC function.

The host immune response to gastrointestinal nematode infection in sheep.

Gastrointestinal nematode infection represents a major threat to the health, welfare and productivity of sheep populations worldwide. Infected lambs have a reduced ability to absorb nutrients from the gastrointestinal tract, resulting in morbidity and occasional mortality. The current chemo-dominant approach to nematode control is considered unsustainable due to the increasing incidence of anthelmintic resistance. In addition, there is growing consumer demand for food products from animals not subjected to chemical treatment. Future mechanisms of nematode control must rely on alternative, sustainable strategies such as vaccination or selective breeding of resistant animals. Such strategies take advantage of the host's natural immune response to nematodes. The ability to resist gastrointestinal nematode infection is considered to be dependent on the development of a protective acquired immune response, although the precise immune mechanisms involved in initiating this process remain to be fully elucidated. In this study, current knowledge on the innate and acquired host immune response to gastrointestinal nematode infection in sheep and the development of immunity is reviewed.

Genome-wide association and regional heritability mapping to identify loci underlying variation in nematode resistance and body weight in Scottish Blackface lambs.

The genetic architecture underlying nematode resistance and body weight in Blackface lambs was evaluated comparing genome-wide association (GWA) and regional heritability mapping (RHM) approaches. The traits analysed were faecal egg count (FEC) and immunoglobulin A activity against third-stage larvae from Teladorsagia circumcincta, as indicators of nematode resistance, and body weight in a population of 752 Scottish Blackface lambs, genotyped with the 50k single-nucleotide polymorphism (SNP) chip. FEC for both Nematodirus and Strongyles nematodes (excluding Nematodirus), as well as body weight were collected at approximately 16, 20 and 24 weeks of age. In addition, a weighted average animal effect was estimated for both FEC and body weight traits. After quality control, 44 388 SNPs were available for the GWA analysis and 42 841 for the RHM, which utilises only mapped SNPs. The same fixed effects were used in both analyses: sex, year, management group, litter size and age of dam, with day of birth as covariate. Some genomic regions of interest for both nematode resistance and body weight traits were identified, using both GWA and RHM approaches. For both methods, strong evidence for association was found on chromosome 14 for Nematodirus average animal effect, chromosome 6 for Strongyles FEC at 16 weeks and chromosome 6 for body weight at 16 weeks. Across the entire data set, RHM identified more regions reaching the suggestive level than GWA, suggesting that RHM is capable of capturing some of the variation not detected by GWA analyses.

Selective forces shaping diversity in the class I region of the major histocompatibility complex in dairy cattle.

The major histocompatibility complex (MHC) is one of the most diverse regions of the mammalian genome. Diversity in MHC genes is integral to their function in the immune system, and while pathogens play a key role in shaping this diversity, the contribution of other selective forces remains unclear. The controlled breeding of cattle offers an excellent model for the identification and exploration of these forces. We characterized the MHC class I genes present in a sample of Canadian Holstein A.I. bulls and compared the results with those obtained in an earlier study. No evidence for a reduction in MHC diversity over 20 years was observed, but the relative frequency of some haplotypes had changed: the formerly rare A12 (w12B) haplotype had become the most common, together with A15, while A19, which dominated the earlier sample, had significantly reduced in frequency. Only 7% of bulls in the current study were MHC homozygous compared with the 14% expected under Hardy-Weinberg. To identify the selective forces at work, a gene substitution model was used to calculate the effects of MHC on selection traits using estimated breeding values for each bull. Significant associations between MHC and production, disease and fertility traits were identified, suggesting that MHC diversity is not merely shaped by disease in this controlled breeding system. The decrease in a common haplotype, the reduced number of homozygous bulls and the associations with disease and production traits together indicate that MHC diversity in dairy cattle is maintained by heterozygote advantage.

Conserved haplotype blocks within the sheep MHC and low SNP heterozygosity in the Class IIa subregion.

This report describes single-nucleotide polymorphisms (SNPs) in the sheep major histocompatibility complex (MHC) class II and class III regions and provides insights into the internal structure of this important genomic complex. MHC haplotypes were deduced from sheep family trios based on genotypes from 20 novel SNPs representative of the class II region and 10 previously described SNPs spanning the class III region. All 30 SNPs exhibited Hardy-Weinberg proportions in the sheep population studied. Recombination within an extended sire haplotype was observed within the class II region for 4 of 20 sheep chromosomes, thereby supporting the presence of separated IIa and IIb subregions similar to those present in cattle. SNP heterozygosity varied across the class II and III regions. One segment of the class IIa subregion manifested very low heterozygosity for several SNPs spanning approximately 120 Kbp. This feature corresponds to a subregion within the human MHC class II region previously described as a 'SNP desert' because of its paucity of SNPs. Linkage disequilibrium (LD) was reduced at the junction separating the putative class IIb and IIa subregions and also between the class IIa and the class III subregions. The latter observation is consistent with either an unmapped physical separation at this location or more likely a boundary characterized by more frequent recombination between two conserved subregions, each manifesting high within-block LD. These results identify internal blocks of loci in the sheep MHC, within which recombination is relatively rare.

A single nomenclature and associated database for alleles at the major histocompatibility complex class II DRB1 locus of sheep.

The development of standardised nomenclatures with associated databases containing reference sequences for alleles at polymorphic loci within the major histocompatibility complex (MHC) has been facilitated by the development of the immuno polymorphism database (IPD). Recently, included within IPD-MHC is information on allelic diversity within sheep species (IPD-MHC-OLA). Here, we present the first report of progress in populating the sheep IPD-MHC database with alleles at the class II MHC DRB1 locus. The sequence of 63 Ovar-DRB1 alleles within 24 allelic families is now held within the database, each meeting the minimum requirement of a complete second exon. These sequences are derived from a combination of genomic and cDNA-based approaches and represent the most extensive collection of validated alleles at the sheep DRB1 locus yet described. Although these 63 alleles probably represent only a fraction of the DRB1 allelic diversity in sheep species worldwide, we encourage the research community to use the official allelic nomenclature and to contribute allelic sequences to the database via its web-based submission tool. In time, the IPD-MHC-OLA resource will underpin population-based MHC genotyping studies and help to simplify meta-analyses of multi-source data from wild and domestic sheep populations.

A mechanistic model of developing immunity to Teladorsagia circumcincta infection in lambs.

Acquired immunity influences the severity of parasitic disease, but modelling the effects of acquired immunity in helminth infections has proved challenging. This may be due to a lack of suitable immunological data, or to the perceived complexity of modelling the immune response. We have developed a model of T. circumcincta infection in domestic sheep that incorporates the effects of acquired immunity on parasite establishment and fecundity. A large data set from commercially managed populations of Scottish Blackface sheep was used, which included relationships between IgA activity and worm length, and between worm length and fecundity. Use was also made of a recently published meta-analysis of parasite establishment rates. This realistic but simple model of nematode infection emulates observed patterns of faecal egg counts. The end-of-season faecal egg counts are remarkably robust to perturbations in the majority of the parameters, possibly because of priming of the immune system early in the season, reducing parasite establishment and growth and, therefore, faecal egg counts. Lowering the amount of early infection leads to higher end-of-season egg counts. The periparturient rise in egg counts in ewes appears to have an important role in supplying infection for the priming of the immune response. This feedback in the immune priming suggests that nematode infections may be difficult to eliminate.

An evolutionary perspective on gastrointestinal nematodes of sheep.

The purpose of this paper was to discuss from an evolutionary perspective the interaction between domestic sheep (Ovis aries) and their gastrointestinal nematodes. Although evolution is the central theme of biology, there has been little attempt to consider how evolutionary forces have shaped and continue to shape the relationships between domestic animals and their parasite community. Mathematical modelling of the host-parasite relationship indicated that the system is remarkably robust to perturbations in its parameters. This robustness may be a consequence of the long coevolution of host and parasites. Although nematodes can potentially evolve faster than the host, coevolution is not dominated by the parasite and there are several examples where breeds of cattle or sheep have evolved high levels of resistance to disease. Coevolution is a more equal partnership between host and nematode than is commonly assumed. Coevolution between parasites and the host immune system is often described as an arms race where both host immune response genes and parasite proteins evolve rapidly in response to each other. However, initial results indicate that nematode antigens are not evolving rapidly; the arms race between the immune system and nematodes, if it exists, is happening very slowly. Fisher's fundamental theorem of natural selection states that genes with positive effects on fitness will be fixed by natural selection. Consequently, heritable variation in fitness traits is expected to be low. Contrary to this argument, there is considerable genetic variation in resistance to nematode infection. In particular, the heritabilities of nematode-specific IgA and IgE activity are moderate to high. The reasons for this apparent violation of the fundamental theorem of natural selection are not clear but several possible explanations are explored. Faecal nematode egg counts increase at the beginning of the grazing season - a phenomenon known as the periparturient rise. This increase benefits host and parasite and appears to be a consequence of coevolution. In conclusion, an evolutionary perspective can shed light on many aspects of the host-parasite relationship in domestic animals.

Genetic variation among lambs in peripheral IgE activity against the larval stages of Teladorsagia circumcincta.

IgA and IgE activity against Teladorsagia circumcincta was investigated in a flock of Texel lambs following natural, mixed nematode infection among lambs. The distribution of IgA activity was similar to a gamma distribution whereas IgE activity was different. Box-Cox analysis demonstrated that X0.25 was a suitable transformation to normalise IgE responses. The transformed IgE activity was under moderate to strong genetic control. Nine different allergens were identified by proteomic analysis. Tropomyosin was selected for further analysis. IgE activity against tropomyosin was moderately heritable and associated with decreased egg counts and with reduced body weight at the time of sampling.

Immunoglobulins as Biomarkers for Gastrointestinal Nematodes Resistance in Small Ruminants: A systematic review.

The rise of anthelmintic resistance worldwide has led to the development of alternative control strategies for gastrointestinal nematodes (GIN) infections, which are one of the main constraints on the health of grazing small ruminants. Presently, breeding schemes rely mainly on fecal egg count (FEC) measurements on infected animals which are time-consuming and requires expertise in parasitology. Identifying and understanding the role of immunoglobulins in the mechanisms of resistance could provide a more efficient and sustainable method of identifying nematode-resistant animals for selection. In this study we review the findings on immunoglobulin response to GIN in the literature published to date (june 2019) and discuss the potential to use immunoglobulins as biomarkers. The literature review revealed 41 studies which measured at least one immunoglobulin: 35 focused on lamb immune response (18 used non-naïve lambs) and 7 on yearlings. In this review we propose a conceptual model summarizing the role of immunoglobulins in resistance to GIN. We highlight the need for more carefully designed and documented studies to allow comparisons across different populations on the immunoglobulin response to GIN infection.

Genetic variation in resistance to mixed, predominantly Teladorsagia circumcincta nematode infections of sheep: from heritabilities to gene identification.

In cool temperate areas, such as Scotland, sheep are infected by a variety of nematodes but the dominant nematode is Teladorsagia circumcincta. Resistant animals have one or more of the following features: fewer adult nematodes, more inhibited larvae, shorter adult nematodes and decreased production of nematode eggs. In lambs at the end of the first grazing season, the heritability of adult worm length is very strong, whereas the heritability of egg production is moderate. The heritability of worm number is low while there is no detectable genetic variation in the number of inhibited larvae. The major mechanisms underlying resistance to T. circumcincta appear to be the IgA mediated suppression of worm growth and the mast cell mediated regulation of worm number. Mast cell responses are slow to develop, possibly because they are responsible for protein loss and reduced growth of the host. Two genes have been repeatedly associated with resistance to T. Circumcincta: the MHC class II DRB1 locus on chromosome 20 and the interferon-gamma locus on chromosome 3. Although the causative mutations are still unknown both genes are plausible candidates.

Effects of host characteristics and parasite intensity on growth and fecundity of Trichostrongylus retortaeformis infections in rabbits.

Host-mediated responses and parasite density-dependent processes can have a major influence on the growth and fecundity of parasitic nematodes. However, host characteristics and parasite intensity consistently change during the course of an infection and these could affect worm length and number of eggs in a non-constant way. We used a free-living population of rabbits naturally infected with Trichostronglyus retortateformis and examined how adult nematode length and the number of eggs in utero were affected by host characteristics (i.e. age, sex, breeding status) and parasite intensity, in a seasonal environment, between 2004 and 2007. Nematode length and number of eggs in utero decreased exponentially with host age; in contrast, parasite intensity increased, peaked in juveniles and decreased in older hosts. These patterns were consistent between rabbit cohorts. A negative relationship was observed between parasite intensity and nematode length, as well as number of eggs. Nematode length was strongly affected by nematode sex and host age, while the number of eggs was mainly influenced by nematode length. The direct influence of host-mediated effects appeared quantitatively more important than parasite density dependence in controlling length and egg production in naturally infected wild rabbits. However, their relative contribution changed during the course of the infection such that, while host immunity still influenced worm numbers, the direct effect of density-dependent interactions contributed the most at high parasite intensities.

The distribution of the pathogenic nematode Nematodirus battus in lambs is zero-inflated.

Understanding the frequency distribution of parasites and parasite stages among hosts is essential for efficient experimental design and statistical analysis, and is also required for the development of sustainable methods of controlling infection. Nematodirus battus is one of the most important organisms that infect sheep but the distribution of parasites among hosts is unknown. An initial analysis indicated a high frequency of animals without N. battus and with zero egg counts, suggesting the possibility of a zero-inflated distribution. We developed a Bayesian analysis using Markov chain Monte Carlo methods to estimate the parameters of the zero-inflated negative binomial distribution. The analysis of 3000 simulated data sets indicated that this method out-performed the maximum likelihood procedure. Application of this technique to faecal egg counts from lambs in a commercial upland flock indicated that N. battus counts were indeed zero-inflated. Estimating the extent of zero-inflation is important for effective statistical analysis and for the accurate identification of genetically resistant animals.

Associations between polymorphisms in the porcine haptoglobin gene and baseline levels of serum haptoglobin.

Haptoglobin (Hp) is an acute phase protein that is a marker in blood for clinical and subclinical disease in the pig. The aim of this study was to identify single nucleotide polymorphisms (SNPs) in the Hp gene and analyse their influence on baseline serum levels. DNA samples and serum were collected from 345 boars. Of 13 SNPs identified, 5 were genotyped using PCR-RFLP and Pyrosequencing. Serum Hp levels were measured using a biochemical assay. A general linear model was fitted with line and genotype as fixed effects. In addition, linkage disequilibrium (LD) was estimated between the 5 SNPs using r-square and D prime. Serum Hp concentrations in the population showed a skewed distribution with a mean of 0.34 g/L (range 0-2.65 g/L). Three SNPs were found to be associated with baseline Hp levels (p-value = 0.0093, 0.0051 and 0.0094). These 3 SNPs were also found to be in high LD with each other. This is the first study to find associations between polymorphisms in the porcine Hp gene and baseline Hp serum levels. The results have implications for breeding for resistance to infection.

Detection of genes with moderate effects on disease resistance using ovine mhc and resistance to nematodes as an example.

Detecting some of the genes that influence disease resistance would improve our understanding of the processes that cause disease and also simplify disease control. Genes within the major histocompatibility complex (mhc) are strong candidates for disease resistance and they have been intensely studied for the last 30 years. Recently, several groups working independently have reported the existence of alleles within the mhc that are associated with enhanced resistance to nematode infection. This article uses hindsight to describe some of the potential pitfalls that hinder the search for valid disease resistance genes. The search requires a good understanding of disease biology, molecular genetics, statistical genetics and especially, the design and analysis of experiments. The power to detect mhc effects is quite low and is quite sensitive to the frequency of the putative resistance alleles.

Influence of rearing conditions and respiratory disease on haptoglobin levels in the pig at slaughter.

Associations between serum concentrations of haptoglobin, pathological lung lesions indicative of Mycoplasma hyopneumoniae (EP) or Actinobacillus pleuropneumoniae (PL) infection at slaughter and previous rearing environment were investigated in 510 pigs (90-100 kg live weight) from 17 farms in England. Haptoglobin concentrations were significantly higher in pigs showing pathological signs of EP infection compared to those without signs of this disease (EP positive median 0.43 mg ml(-1) vs. EP negative median 0.26 mg ml(-1), p<0.01). However, there were no significant associations between serum haptoglobin concentrations and pathological signs of PL. The presence of solid partitions compared with barred or similar open partitions was associated with a decrease of 0.44 mg ml(-1) farm mean haptoglobin concentration, whilst an increase in pen size of 10 m(2) was associated with a decrease of 0.08 mg ml(-1) farm mean haptoglobin concentration. The findings indicate that pathological signs of EP were associated with increased serum haptoglobin at slaughter, which in turn was influenced by components of the farm environment.