The distribution of SNP marker effects for faecal worm egg count in sheep, and the feasibility of using these markers to predict genetic merit for resistance to worm infections

Genomic Regions Associated with Resistance to Gastrointestinal Parasites in Australian Merino Sheep

The objective of this study was to identify genomic regions and genes associated with resistance to gastrointestinal nematodes in Australian Merino sheep in Uruguay, using the single-step GWAS methodology (ssGWAS), which is based on genomic estimated breeding values (GEBVs) obtained from a combination of pedigree, genomic, and phenotypic data. This methodology converts GEBVs into SNP effects. The analysis included 26,638 animals with fecal egg count (FEC) records obtained in two independent parasitic cycles (FEC1 and FEC2) and 1700 50K SNP genotypes. The comparison of genomic regions was based on genetic variances (gVar(%)) explained by non-overlapping regions of 20 SNPs. For FEC1 and FEC2, 18 and 22 genomic windows exceeded the significance threshold (gVar(%) ≥ 0.22%), respectively. The genomic regions with strong associations with FEC1 were located on chromosomes OAR 2, 6, 11, 21, and 25, and for FEC2 on OAR 5, 6, and 11. The proportion of genetic variance attributed to the top windows was 0.83% and 1.9% for FEC1 and FEC2, respectively. The 33 candidate genes shared between the two traits were subjected to enrichment analysis, revealing a marked enrichment in biological processes related to immune system functions. These results contribute to the understanding of the genetics underlying gastrointestinal parasite resistance and its implications for other productive and welfare traits in animal breeding programs.

Reducing fecal egg count through selective breeding alters dorper lamb response to Haemonchus contortus in an artificial challenge trial

Infection by gastrointestinal nematodes (GIN), particularly Haemonchus contortus, can be detrimental to sheep health and performance. Genetic susceptibility to GIN varies between breeds, with those lacking high levels of natural resistance often requiring frequent anthelmintic treatment when facing parasitic challenge. Genetic technology can serve as a tool to decrease GIN susceptibility via selection for sheep with reduced fecal egg count (FEC) estimated breeding values (EBVs). However, the physiological changes that result from implementation of this strategy are not well described. Additionally, there is a need for comparison of animals from recent selective breeding against breeds with inherent GIN resistance. In this study we administered a challenge of H. contortus to Dorper x White Dorper (DWD; n = 92) lambs that have been genetically selected for either low (DWD-) or high (DWD+) FEC EBVs and Barbados Blackbelly x Mouflon (BBM; n = 19) lambs from a genetically resistant breed backgrounds. Lamb FEC, packed-cell volume (PCV) and serum IgG were measured at intermittent levels over 5 weeks. At day 21 and day 35, the selectively bred DWD- had a lower mean FEC compared to DWD+, but were higher than BBM. Reductions in both PCV and serum IgG from initial day 0 levels were observed in DWD lambs, but not in BBM. Furthermore, from a subset of lambs (n = 24) harvested at day 21, DWD- only tended (p = 0.056) to have lower mean worm counts than DWD+, with BBM having the lowest mean worm count. Differentially expressed genes (DEGs) identified via RNA-sequencing of abomasal tissue at day 21 indicate a more pronounced Th2 immune response and more rapid worm expulsion occurred in iBBM than iDWD- and iDWD+ lambs. However, gene expression in DWD- suggests an association between reduced FEC EBV and gastric acid secretion and the ability to limit worm fecundity. Ultimately, selection of Dorper sheep for low FEC EBV can reduce susceptibility to GIN, but it will likely require multiple generations with this trait as a breeding priority before presenting a similar resistance level to Caribbean breeds.

Ovine PAPPA2 gene coding variants are linked to decreased fecal egg shedding in native Turkish sheep naturally infected with gastrointestinal nematodes

In this study, the association between PAPPA2 coding variants and gastrointestinal (GI) nematode fecal egg count (FEC) score in adult Turkish sheep was investigated. For this purpose, the FEC score was determined in adult sheep from six breeds: Karacabey Merino (n = 137), Kivircik (n = 116), Cine capari (n = 109), Karakacan (n = 102), Imroz (n = 73), and Chios (n = 50). Sheep were classified as shedders or non-shedders within breeds and flocks. The first group was the fecal egg shedders (> 50 per gram of feces), and the second group was the no fecal egg shedders (≤ 50 per gram of feces). The exon 1, exon 2, exon 5, exon 7, and a part of 5'UTR of the ovine PAPPA2 gene were genotyped by Sanger sequencing of these two groups. Fourteen synonymous and three non-synonymous single-nucleotide polymorphisms (SNPs) were found. The non-synonymous SNPs, D109N, D391H, and L409R variants, are reported for the first time. Two haplotype blocks were constructed on exon 2 and exon 7. The specific haplotype, C₃₉₁G₄₂₄G₄₄₉T₄₇₃C₅₁₅A₅₄₂ on the exon 2 that carries the 391H variant, was tested against four other common haplotypes. Our results indicate that C₃₉₁G₄₂₄G₄₄₉T₄₇₃C₅₁₅A₅₄₂ haplotype was significantly associated with fecal egg shedding status in adult Turkish sheep (p-value, 0.044).

Prevalence and resistance to gastrointestinal parasites in goats: A review

Gastrointestinal parasitism, particularly nematode infection, is a major health issue affecting goats worldwide, resulting in clinical diseases and productivity loss. Prevalent gastrointestinal parasites (GIPs) affecting goats in South Africa are the Strongyloides papillosus, Eimeria spp., and Strongyles, especially the Haemonchus contortus and Trichostrongylus spp. According to the issues discussed in this paper and by other authors, the prevalence and intensity of various GIPs vary with an animal's location, breed, age, sex, and season. Because GIPs easily develop resistance to chemical treatment, selecting and breeding genetically GIP-resistant animals would be a relatively simple and inexpensive strategy for reducing or eliminating the current reliance on chemotherapy. Potential phenotypic indicators for selecting GIP-resistant goats include parasitological, immunological, and pathological phenotypic markers. Synergistic use of these indicators should be encouraged for a more accurate simplified genotype selection of resistant animals. Genes with Mendelian inheritance, particularly those involved in immunoregulatory mechanisms, have been identified in goats. Exploring this knowledge base to develop cost-effective molecular tools that facilitate enhanced genetic improvement programs is a current challenge. Future statistical and biological models should investigate genetic variations within genomic regions and different candidate genes involved in immunoregulatory mechanisms, as well as the identification of single nucleotide polymorphisms known to affect GIP infection levels.

Variants Within Genes EDIL3 and ADGRB3 are Associated With Divergent Fecal Egg Counts in Katahdin Sheep at Weaning

Gastrointestinal nematodes (GIN) pose a severe threat to sheep production worldwide. Anthelmintic drug resistance coupled with growing concern regarding potential environmental effects of drug use have demonstrated the necessity of implementing other methods of GIN control. The aim of this study was to test for genetic variants associated with resistance or susceptibility to GIN in Katahdin sheep to improve the current understanding of the genetic mechanisms responsible for host response to GIN. Linear regression and case-control genome-wide association studies were conducted with high-density genotype data and cube-root transformed weaning fecal egg counts (tFEC) of 583 Katahdin sheep. The case-control GWAS identified two significant SNPs (P-values 1.49e-08 to 1.01e-08) within introns of the gene adhesion G protein-coupled receptor B3 (ADGRB3) associated with lower fecal egg counts. With linear regression, four significant SNPs (P-values 7.82e-08 to 3.34e-08) were identified within the first intron of the gene EGF-like repeats and discoidin domains 3 (EDIL3). These identified SNPs were in very high linkage disequilibrium (r² of 0.996–1), and animals with alternate homozygous genotypes had significantly higher median weaning tFEC phenotypes compared to all other genotypes. Significant SNPs were queried through public databases to identify putative transcription factor binding site (TFBS) and potential lncRNA differences between reference and alternate alleles. Changes in TFBS were predicted at two SNPs, and one significant SNP was found to be within a predicted lncRNA sequence with greater than 90% similarity to a known lncRNA in the bovine genome. The gene EDIL3 has been described in other species for its roles in the inhibition and resolution of inflammation. Potential changes of EDIL3 expression mediated through lncRNA expression and/or transcription factor binding may impact the overall immune response and reduce the ability of Katahdin sheep to control GIN infection. This study lays the foundation for further research of EDIL3 and ADGRB3 towards understanding genetic mechanisms of susceptibility to GIN, and suggests these SNPs may contribute to genetic strategies for improving parasite resistance traits in sheep.

Association analysis and functional annotation of imputed sequence data within genomic regions influencing resistance to gastro-intestinal parasites detected by an LDLA approach in a nucleus flock of Sarda dairy sheep

Background

Gastroinestinal nematodes (GIN) are one of the major health problem in grazing sheep. Although genetic variability of the resistance to GIN has been documented, traditional selection is hampered by the difficulty of recording phenotypes, usually fecal egg count (FEC). To identify causative mutations or markers in linkage disequilibrium (LD) to be used for selection, the detection of quantitative trait loci (QTL) for FEC based on linkage disequilibrium-linkage analysis (LDLA) was performed on 4097 ewes (from 181 sires) all genotyped with the OvineSNP50 Beadchip. Identified QTL regions (QTLR) were imputed from whole-genome sequences of 56 target animals of the population. An association analysis and a functional annotation of imputed polymorphisms in the identified QTLR were performed to pinpoint functional variants with potential impact on candidate genes identified from ontological classification or differentially expressed in previous studies.

Results

After clustering close significant locations, ten QTLR were defined on nine Ovis aries chromosomes (OAR) by LDLA. The ratio between the ANOVA estimators of the QTL variance and the total phenotypic variance ranged from 0.0087 to 0.0176. QTL on OAR4, 12, 19, and 20 were the most significant. The combination of association analysis and functional annotation of sequence data did not highlight any putative causative mutations. None of the most significant SNPs showed a functional effect on genes’ transcript. However, in the most significant QTLR, we identified genes that contained polymorphisms with a high or moderate impact, were differentially expressed in previous studies, contributed to enrich the most represented GO process (regulation of immune system process, defense response). Among these, the most likely candidate genes were: TNFRSF1B and SELE on OAR12, IL5RA on OAR19, IL17A, IL17F, TRIM26, TRIM38, TNFRSF21, LOC101118999, VEGFA, and TNF on OAR20.

Conclusions

This study performed on a large experimental population provides a list of candidate genes and polymorphisms which could be used in further validation studies. The expected advancements in the quality of the annotation of the ovine genome and the use of experimental designs based on sequence data and phenotypes from multiple breeds that show different LD extents and gametic phases may help to identify causative mutations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00690-7.

Genetic × environment variation in sheep lines bred for divergent resistance to strongyle infection

Drug-resistant parasites threaten livestock production. Breeding more resistant hosts could be a sustainable control strategy. Environmental variation linked to animal management practices or to parasite species turnover across farms may however alter the expression of genetic potential. We created sheep lines with high or low resistance to Haemonchus contortus and achieved significant divergence on both phenotypic and genetic scales. We exposed both lines to chronic stress or to the infection by another parasite Trichostrongylus colubriformis, to test for genotype-by-environment and genotype-by-parasite species interactions respectively. Between-line divergence remained significant following chronic stress exposure although between-family variation was found. Significant genotype-by-parasite interaction was found although H. contortus-resistant lambs remained more resistant against T. colubriformis. Growth curves were not altered by the selection process although resistant lambs were lighter after the second round of divergence, before any infection took place. Breeding for resistance is a sustainable strategy but allowance needs to be made for environmental perturbations and worm species.

Integrative biology defines novel biomarkers of resistance to strongylid infection in horses

The widespread failure of anthelmintic drugs against nematodes of veterinary interest requires novel control strategies. Selective treatment of the most susceptible individuals could reduce drug selection pressure but requires appropriate biomarkers of the intrinsic susceptibility potential. To date, this has been missing in livestock species. Here, we selected Welsh ponies with divergent intrinsic susceptibility (measured by their egg excretion levels) to cyathostomin infection and found that their divergence was sustained across a 10-year time window. Using this unique set of individuals, we monitored variations in their blood cell populations, plasma metabolites and faecal microbiota over a grazing season to isolate core differences between their respective responses under worm-free or natural infection conditions. Our analyses identified the concomitant rise in plasma phenylalanine level and faecal Prevotella abundance and the reduction in circulating monocyte counts as biomarkers of the need for drug treatment (egg excretion above 200 eggs/g). This biological signal was replicated in other independent populations. We also unravelled an immunometabolic network encompassing plasma beta-hydroxybutyrate level, short-chain fatty acid producing bacteria and circulating neutrophils that forms the discriminant baseline between susceptible and resistant individuals. Altogether our observations open new perspectives on the susceptibility of equids to strongylid infection and leave scope for both new biomarkers of infection and nutritional intervention.

Detection of genomic regions underlying milk production traits in Valle del Belice dairy sheep using regional heritability mapping

The aim of this study was to identify genomic regions underlying milk production traits in the Valle del Belice dairy sheep using regional heritability mapping (RHM). Repeated measurements for milk yield (MY), fat percentage and yield (F% and FY) and protein percentage and yield (P% and PY), collected over a period of 6 years (2006–2012) on 481 Valle del Belice ewes, were used for the analysis. Animals were genotyped with the Illumina 50k SNP chip. Variance components, heritabilities and repeatabilities within and across lactations were estimated, fitting parity, litter size, season of lambing and fortnights in milk, as fixed; and additive genetic, permanent environment within and across lactations, flock by test‐day interaction and residual as random effects. For the RHM analysis, the model included the same fixed and random effects as before, plus an additional regional genomic additive effect (specific for the region being tested) as random. While the whole genomic additive effect was estimated using the genomic relationship matrix (GRM) constructed from all SNPs, the regional genomic additive effect was estimated from a GRM matrix constructed from the SNPs within each region. Heritability estimates ranged between 0.06 and 0.15, with repeatabilities being between 0.14 and 0.24 across lactations and between 0.23 and 0.39 within lactation for all milk production traits. A substantial effect of flock‐test‐day on milk production traits was also estimated. Significant genomic regions at either genome‐wide (p < .05) or suggestive (i.e., one false positive per genome scan) level were identified on chromosome (OAR) 2, 3 and 20 for F% and on OAR3 for P%, with the regions on OAR3 in common between the two traits. Our results confirmed the role of LALBA and AQP genes, on OAR3, as candidate genes for milk production traits in sheep.

Genome-wide insights on gastrointestinal nematode resistance in autochthonous Tunisian sheep

Gastrointestinal nematode (GIN) infections have negative impacts on animal health, welfare and production. Information from molecular studies can highlight the underlying genetic mechanisms that enhance host resistance to GIN. However, such information often lacks for traditionally managed indigenous livestock. Here, we analysed 600 K single nucleotide polymorphism genotypes of GIN infected and non-infected traditionally managed autochthonous Tunisian sheep grazing communal natural pastures. Population structure analysis did not find genetic differentiation that is consistent with infection status. However, by contrasting the infected versus non-infected cohorts using ROH, LR-GWAS, F_ST and XP-EHH, we identified 35 candidate regions that overlapped between at least two methods. Nineteen regions harboured QTLs for parasite resistance, immune capacity and disease susceptibility and, ten regions harboured QTLs for production (growth) and meat and carcass (fatness and anatomy) traits. The analysis also revealed candidate regions spanning genes enhancing innate immune defence (SLC22A4, SLC22A5, IL-4, IL-13), intestinal wound healing/repair (IL-4, VIL1, CXCR1, CXCR2) and GIN expulsion (IL-4, IL-13). Our results suggest that traditionally managed indigenous sheep have evolved multiple strategies that evoke and enhance GIN resistance and developmental stability. They confirm the importance of obtaining information from indigenous sheep to investigate genomic regions of functional significance in understanding the architecture of GIN resistance.

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.

Detecting fecal egg count (FEC) for gastrointestinal nematodes of adult Turkish sheep with different scrapie related PRNP haplotypes

Scrapie is a transmissible spongiform encephalopathy caused by prions and leads to neurodegeneration in the Central Nervous System (CNS) of sheep and goats. Genetic resistance/susceptibility to scrapie is well studied and it is known that the variations of 136th, 154th and 171st codons at the ovine PRNP gene have a major effect on the development of the disease. Many studies demonstrated that selection for PRNP genotypes has not influenced other performance traits, nevertheless, there is a knowledge gap about the possible link between the PRNP gene and the status of the other important diseases that affect the sheep population worldwide. In the present study, we tested whether there is an association between scrapie-related PRNP genotypes and fecal egg count (FEC) of gastrointestinal nematodes in seven adult Turkish sheep breeds. For this purpose, FEC scores of studied sheep (n = 253) were determined and the same animals were genotyped for the PRNP gene. Finally, an association analysis was performed for scrapie resistant (ARR), susceptible (VRQ), and wild-type (ARQ) haplotypes. Based on our statistical analysis, it is concluded that PRNP genotypes have no positive or negative effect on the FEC scores of adult sheep.

Single nucleotide polymorphisms from candidate genes associated with nematode resistance and resilience in Corriedale and Pampinta sheep in Argentina

Selective breeding of genetically resistant animals is considered a promising strategy to face the problem of nematode resistance to anthelmintics and mitigate concerns about the presence of chemical residues in animal food products and the environment. Gastrointestinal nematode resistance is a complex, multifactorial trait related to host immunity. However, the mechanisms underlying host resistance and response to infection remain to be fully elucidated. In this context, the objective of this study was to provide insight into the chromosomal regions determining nematode resistance and resilience in Corriedale and resistance in Pampinta sheep breeds. A total of 170 single nucleotide polymorphisms (SNP) from 76 candidate genes for immune response were studied in 624 Corriedale and 304 Pampinta animals. Lambs underwent artificial or natural challenges with infective larvae mainly from Haemonchus contortus. Fecal egg counts, estimated breeding values for fecal egg counts, and rate of packed cell volume change and FAMACHA© score change over the challenge were used, when available, as indicators of host parasite resistance or resilience. Phenotype-genotype association studies were conducted and significance values obtained were adjusted for multiple testing errors. Eight SNPs, located on OARs 3, 6, 12, and 20, reached significance in Corriedale sheep under artificial challenge. Those SNP represent allelic variants from the MHC-Ovine Lymphocyte Antigen-DRA, two C-type lectin domain families, the Interleukin 2 receptor β, the Toll-like receptor 10, the Mannan binding lectin serine peptidase 2, and the NLR family, CARD domain containing 4 genes. On Pampinta lambs under natural challenge, we found three significant SNPs, located in the TIMP metallopeptidase inhibitor 3, the FBJ murine osteosarcoma viral oncogene homolog, and the Interleukin 20 receptor alpha genes, on OARs 3, 7, and 8, respectively. The results obtained herein confirm genomic regions previously reported as associated with nematode resistance in other sheep breeds, reinforcing their role in host response to parasites. These findings contribute to gain knowledge on parasite resistance and resilience in Corriedale sheep and report for the first time SNPs associated with resistance to gastrointestinal parasite infections in Pampinta breed.

Four single nucleotide polymorphisms (SNPs) are associated with resistance and resilience to Haemonchus contortus in Brazilian Morada Nova sheep

Gastrointestinal nematodes are a major constraint in sheep production. Breeding for resistance has proven to be an effective and feasible approach to address this problem. The use and investigation of genetic markers for resistance traits could accelerate genetic progress and lead to a better understanding of underlying molecular mechanisms. Thus, the aim of this study was to evaluate if five single nucleotide polymorphisms SNPs OAR2_14765360, OAR6_81718546, OAR11_62887032, OAR12_69606944 and OAR15_59871543 are associated with resistance and resilience traits in a flock of the Morada Nova sheep breed. Lambs were submitted to two consecutive parasite challenges by oral infection with 4000 infective larvae L₃) of Haemonchus contortus. Fecal egg counts (FEC), packed cell volume (PVC) and body weight were measured every one or two weeks for 42 days in each trial. DNA samples from 287 lambs, 131 ewes and 4 rams were amplified by ARMS-PCR or PCR-RFLP and genotypes were determined. Analysis of variance (ANOVA) was used for association analyses between genotypes and phenotypes. In case of significant association, the allele substitution effect was calculated based on a linear model. OAR2_14765360 and OAR12_69606944 were associated with FEC, and OAR12_69606944 also had significant effects on PCV and weight gain, showing favourable associations of the CC genotype with all evaluated traits. Both OAR6_81718546 and OAR11_62887032 were associated with weight gain, and OAR6_81718546 had an additional effect on PCV. OAR15_59871543 was not polymorphic in the population. OAR6_81718546 and OAR12_69606944 presented significant allele substitution effects of -1.06 ± 0.52 kg for the T allele on final body weight and 0.74 ± 0.32 for the C allele in PCV of the same sampling date, respectively. This is the first report of SNPs associated with gastrointestinal nematode resistance in this sheep breed. Our findings support the existence of quantitatice trait loci (QTL) for resistance and resilience in linkage disequilibrium with the polymorphic SNPs and suggest their future use for explorations of these traits in Morada Nova sheep.

Association study reveals Th17, Treg, and Th2 loci related to resistance to Haemonchus contortus in Florida Native sheep1

The aim of this study was to identify for the first time single nucleotide polymorphisms (SNPs) associated with Haemonchus contortus resistance in Florida Native sheep, using a targeted sequencing approach. One hundred and fifty-three lambs were evaluated in this study. At the start of the trial, phenotypic records for fecal egg count (FEC), FAMACHA score, body condition score (BCS), and weight were recorded and deworming of sheep with levamisole (18 mg/kg of body weight) was performed. Ten days post-deworming (baseline) and 28 d post-baseline, a full hematogram of each sheep was obtained and FEC, FAMACHA score, BCS, and weight were assessed. Average daily gain was calculated at the end of the trial. Out of 153 animals, 100 sheep were selected for genotyping using a targeted sequencing approach. Targeted sequencing panel included 100 candidate genes for immune response against H. contortus. SNPs were discarded if call rate <95% and minor allele frequency ≤0.05. A mixed model was used to analyze the response variables and included the identity by state matrix to control for population structure. A contemporary group (age, group, and sex) was included as fixed effect. Bonferroni correction was used to control for multiple testing. Eighteen SNPs on chromosomes 1, 2, 3, 4, 6, 7, 11, 15, 18, 20, 24, and 26 were significant for different traits. Our results suggest that loci related to Th17, Treg, and Th2 responses play an important role in the expression of resistant phenotypes. Several genes including ITGA4, MUC15, TLR3, PCDH7, CFI, CXCL10, TNF, CCL26, STAT3, GPX2, IL2RB, and STAT6 were identified as potential markers for resistance to natural H. contortus exposure. This is the first study that evaluates potential genetic markers for H. contortus resistance in Florida Native sheep.

Signatures of selection for resistance to Haemonchus contortus in sheep and goats

Background

Gastrointestinal nematode infection (GNI) is the most important disease affecting the small ruminant industry in U.S. The environmental conditions in the southern United States are ideal for the survival of the most pathogenic gastrointestinal nematode, Haemonchus contortus. Host genetic variation for resistance to H. contortus allows selective breeding for increased resistance of animals. This selection process increases the prevalence of particular alleles in sheep and goats and creates unique genetic patterns in the genome of these species. The aim of this study was to identify loci with divergent allelic frequencies in a candidate gene panel of 100 genes using two different approaches (frequentist and Bayesian) to estimate Fst outliers in three different breeds of sheep and goats exposed to H. contortus.

Results

Our results for sheep populations showed SNPs under selection in C3AR1, CSF3, SOCS2, NOS2, STAT5B, TGFB2 and IL2RA genes using frequentist and Bayesian approaches. For goats, SNPs in CD1D, ITGA9, IL12A, IL13RA1, CD86 and TGFB2 genes were under selection. Common signatures of selection in both species were observed in NOS2, TGFB2 and TLR4 genes. Directional selection was present in all SNPs evaluated in the present study.

Conclusions

A total of 13 SNPs within 7 genes of our candidate gene panel related to H. contortus exposure were identified under selection in sheep populations. For goats, 11 SNPs within 7 genes were identified under selection. Results from this study support the hypothesis that resistance to H. contortus is likely to be controlled by many loci. Shared signatures of selection related to mechanisms of immune protection against H. contortus infection in sheep and goats could be useful targets in breeding programs aimed to produce resistant animals with low FEC.

Detection of genomic regions underlying resistance to gastrointestinal parasites in Australian sheep

BACKGROUND

This study aimed at identifying genomic regions that underlie genetic variation of worm egg count, as an indicator trait for parasite resistance in a large population of Australian sheep, which was genotyped with the high-density 600 K Ovine single nucleotide polymorphism array. This study included 7539 sheep from different locations across Australia that underwent a field challenge with mixed gastrointestinal parasite species. Faecal samples were collected and worm egg counts for three strongyle species, i.e. Teladorsagia circumcincta, Haemonchus contortus and Trichostrongylus colubriformis were determined. Data were analysed using genome-wide association studies (GWAS) and regional heritability mapping (RHM).

RESULTS

Both RHM and GWAS detected a region on Ovis aries (OAR) chromosome 2 that was highly significantly associated with parasite resistance at a genome-wise false discovery rate of 5%. RHM revealed additional significant regions on OAR6, 18, and 24. Pathway analysis revealed 13 genes within these significant regions (SH3RF1, HERC2, MAP3K, CYFIP1, PTPN1, BIN1, HERC3, HERC5, HERC6, IBSP, SPP1, ISG20, and DET1), which have various roles in innate and acquired immune response mechanisms, as well as cytokine signalling. Other genes involved in haemostasis regulation and mucosal defence were also detected, which are important for protection of sheep against invading parasites.

CONCLUSIONS

This study identified significant genomic regions on OAR2, 6, 18, and 24 that are associated with parasite resistance in sheep. RHM was more powerful in detecting regions that affect parasite resistance than GWAS. Our results support the hypothesis that parasite resistance is a complex trait and is determined by a large number of genes with small effects, rather than by a few major genes with large effects.

Enhancing genetic disease control by selecting for lower host infectivity and susceptibility

Infectious diseases have a huge impact on animal health, production and welfare, and human health. Understanding the role of host genetics in disease spread is important for developing disease control strategies that efficiently reduce infection incidence and risk of epidemics. While heritable variation in disease susceptibility has been targeted in livestock breeding, emerging evidence suggests that there is additional genetic variation in host infectivity, but the potential benefits of including infectivity into selection schemes are currently unknown. A Susceptible-Infected-Recovered epidemiological model incorporating polygenic genetic variation in both susceptibility and infectivity was combined with quantitative genetics selection theory to assess the non-linear impact of genetic selection on field measures of epidemic risk and severity. Response to 20 generations of selection was calculated in large simulated populations, exploring schemes differing in accuracy and intensity. Assuming moderate genetic variation in both traits, 50% selection on susceptibility required seven generations to reduce the basic reproductive number R₀ from 7.64 to the critical threshold of <1, below which epidemics die out. Adding infectivity in the selection objective accelerated the decline towards R₀ < 1, to 3 generations. Our results show that although genetic selection on susceptibility reduces disease risk and prevalence, the additional gain from selection on infectivity accelerates disease eradication and reduces more efficiently the risk of new outbreaks, while it alleviates delays generated by unfavourable correlations. In conclusion, host infectivity was found to be an important trait to target in future genetic studies and breeding schemes, to help reducing the occurrence and impact of epidemics.

Exploring the molecular basis of resistance/ susceptibility to mixed natural infection of Haemonchus contortus in tropical Indian goat breed

The present investigation was carried out with the objective to identify putative candidate genes / Quantitative trait loci for resistance / susceptibility towards Haemonchus infestation in tropical goat breed (Rohilkhandi goat) of India. The mean faecal egg count (FEC) and packed cell volume (PCV) of the population were 142.78 ± 22.54 epg (eggs per gram) and 31.73% ± 0.49, respectively. Grouping of animals as per dot ELISA test showed 41.33% (n = 124) positive and 58.66% (n = 176) negative for Haemonchus infestation. The microsatellite loci DYA and ODRB1.2 were significantly associated (P ≤ 0.05) to parasite resistance. The locus DYA showed significant association with log FEC and dot ELISA and the locus ODRB1.2 showed significant association with log FEC, PCV and dot ELISA at P ≤ 0.05. Real time expression profiling revealed that the susceptible group (high FEC group) had 11.1-fold more expression of IFNγ mRNA (Th₁ cytokine) and 0.11-fold lower expression of IL-10 mRNA (Th₂ cytokine), which was found to be statistically significant (P ≤ 0.05).

PCR-based genotyping of SNP markers in sheep

Single nucleotide polymorphisms (SNPs) are the main type of variation in genome, enabling them to be associated with traits of economic importance in livestock. Genome-wide association studies (GWAS) have led to the discovery of SNPs associated with desirable traits in sheep. However, in these studies, SNPs are genotyped by high-throughput methods in genome scale, which are expensive and require sophisticated equipment and analysis methods. Therefore, the goal of this study was to develop a reliable, rapid, and inexpensive polymerase chain reaction (PCR)-based method to genotype a medium number of animals for a few candidate SNPs previously associated with desirable phenotypes in sheep by GWAS, using markers associated with gastrointestinal nematode resistance as a model. DNA extracted from white-blood cells of 150 sheep was submitted to PCR amplification followed by agarose gel electrophoresis and determination of banding pattern. Tetra-primer ARMS-PCR was successfully optimized after changes in annealing temperature; annealing and extension times; concentration of MgCl₂ and DNA; ratios of inner, outer, forward and reverse primer; and addition of adjuvants, for genotyping the OAR2_14765360, OAR6_81718546, OAR11_62887032, and OAR12_69606944 SNPs in sheep. An extensive optimization of tetra-primer ARMS-PCR resulted in a suitable, simple, cost-effective PCR-based method of genotyping four SNP markers previously detected by chip arrays.

Immunogenomics of gastrointestinal nematode infection in ruminants - breeding for resistance to produce food sustainably and safely

Gastrointestinal nematode (GIN) infection of ruminants represents a major health and welfare challenge for livestock producers worldwide. The emergence of anthelmintic resistance in important GIN species and the associated animal welfare concerns have stimulated interest in the development of alternative and more sustainable strategies aimed at the effective management of the impact of GINs. These integrative strategies include selective breeding using genetic/genomic tools, grazing management, biological control, nutritional supplementation, vaccination and targeted selective treatment. In this review, the logic of selecting for "resistance" to GIN infection as opposed to "resilience" or "tolerance" is discussed. This is followed by a review of the potential application of immunogenomics to genetic selection for animals that have the capacity to withstand the impact of GIN infection. Advances in relevant genomic technologies are highlighted together with how these tools can be advanced to support the integration of immunogenomic information into ruminant breeding programmes.

Endurance Exercise Ability in the Horse: A Trait with Complex Polygenic Determinism

Endurance horses are able to run at more than 20 km/h for 160 km (in bouts of 30-40 km). This level of performance is based on intense aerobic metabolism, effective body heat dissipation and the ability to endure painful exercise. The known heritabilities of endurance performance and exercise-related physiological traits in Arabian horses suggest that adaptation to extreme endurance exercise is influenced by genetic factors. The objective of the present genome-wide association study (GWAS) was to identify single nucleotide polymorphisms (SNPs) related to endurance racing performance in 597 Arabian horses. The performance traits studied were the total race distance, average race speed and finishing status (qualified, eliminated or retired). We used three mixed models that included a fixed allele or genotype effect and a random, polygenic effect. Quantile-quantile plots were acceptable, and the regression coefficients for actual vs. expected log10p-values ranged from 0.865 to 1.055. The GWAS revealed five significant quantitative trait loci (QTL) corresponding to 6 SNPs on chromosomes 6, 1, 7, 16, and 29 (two SNPs) with corrected p-values from 1.7 × 10-6 to 1.8 × 10-5. Annotation of these 5 QTL revealed two genes: sortilin-related VPS10-domain-containing receptor 3 (SORCS3) on chromosome 1 is involved in protein trafficking, and solute carrier family 39 member 12 (SLC39A12) on chromosome 29 is active in zinc transport and cell homeostasis. These two coding genes could be involved in neuronal tissues (CNS). The other QTL on chromosomes 6, 7, and 16 may be involved in the regulation of the gene expression through non-coding RNAs, CpG islands and transcription factor binding sites. On chromosome 6, a new candidate equine long non-coding RNA (KCNQ1OT1 ortholog: opposite antisense transcript 1 of potassium voltage-gated channel subfamily Q member 1 gene) was predicted in silico and validated by RT-qPCR in primary cultures of equine myoblasts and fibroblasts. This lncRNA could be one element of the cardiac rhythm regulation. Our GWAS revealed that equine performance during endurance races is a complex polygenic trait, and is partially governed by at least 5 QTL: two coding genes involved in neuronal tissues and three other loci with many regulatory functions such as slowing down heart rate.

Exploring the Genetic Resistance to Gastrointestinal Nematodes Infection in Goat Using RNA-Sequencing

Gastrointestinal nematodes (GINs) are one of the most economically important parasites of small ruminants and a major animal health concern in many regions of the world. However, the molecular mechanisms of the host response to GIN infections in goat are still little known. In this study, two genetically distinct goat populations, one relatively resistant and the other susceptible to GIN infections, were identified in Yichang goat and then four individuals in each group were chosen to compare mRNA expression profiles using RNA-seq. Field experiment showed lower worm burden, delayed and reduced egg production in the relatively resistant group than the susceptible group. The analysis of RNA-seq showed that 2369 genes, 1407 of which were up-regulated and 962 down-regulated, were significantly (p < 0.001) differentially expressed between these two groups. Functional annotation of the 298 genes more highly expressed in the resistant group yielded a total of 46 significant (p < 0.05) functional annotation clusters including 31 genes (9 in innate immunity, 13 in immunity, and 9 in innate immune response) related to immune biosynthetic process as well as transforming growth factor (TGF)-β, mitogen-activated protein kinase (MAPK), and cell adhesion molecules (CAMs) pathways. Our findings provide insights that are immediately relevant for the improvement of host resistance to GIN infections and which will make it possible to know the mechanisms underlying the resistance of goats to GIN infections.

Breeding for resistance to gastrointestinal nematodes - the potential in low-input/output small ruminant production systems

The control of gastrointestinal nematodes (GIN) is mainly based on the use of drugs, grazing management, use of copper oxide wire particles and bioactive forages. Resistance to anthelmintic drugs in small ruminants is documented worldwide. Host genetic resistance to parasites, has been increasingly used as a complementary control strategy, along with the conventional intervention methods mentioned above. Genetic diversity in resistance to GIN has been well studied in experimental and commercial flocks in temperate climates and more developed economies. However, there are very few report outputs from the more extensive low-input/output smallholder systems in developing and emerging countries. Furthermore, results on quantitative trait loci (QTL) associated with nematode resistance from various studies have not always been consistent, mainly due to the different nematodes studied, different host breeds, ages, climates, natural infections versus artificial challenges, infection level at sampling periods, among others. The increasing use of genetic markers (Single Nucleotide Polymorphisms, SNPs) in GWAS or the use of whole genome sequence data and a plethora of analytic methods offer the potential to identify loci or regions associated nematode resistance. Genomic selection as a genome-wide level method overcomes the need to identify candidate genes. Benefits in genomic selection are now being realised in dairy cattle and sheep under commercial settings in the more advanced countries. However, despite the commercial benefits of using these tools, there are practical problems associated with incorporating the use of marker-assisted selection or genomic selection in low-input/output smallholder farming systems breeding schemes. Unlike anthelmintic resistance, there is no empirical evidence suggesting that nematodes will evolve rapidly in response to resistant hosts. The strategy of nematode control has evolved to a more practical manipulation of host-parasite equilibrium in grazing systems by implementation of various strategies, in which improvement of genetic resistance of small ruminant should be included. Therefore, selection for resistant hosts can be considered as one of the sustainable control strategy, although it will be most effective when used to complement other control strategies such as grazing management and improving efficiency of anthelmintics currently.

Transcriptional profiling of the ovine abomasal lymph node reveals a role for timing of the immune response in gastrointestinal nematode resistance

Gastrointestinal nematodes are a serious cause of morbidity and mortality in grazing ruminants. The major ovine defence mechanism is acquired immunity, with protective immunity developing over time in response to infection. Nematode resistance varies both within and between breeds and is moderately heritable. A detailed understanding of the genes and mechanisms involved in protective immunity, and the factors that regulate this response, is required to aid both future breeding strategies and the development of effective and sustainable nematode control methods. The aim of this study was to compare the abomasal lymph node transcriptome of resistant and susceptible lambs in order to determine biological processes differentially expressed between resistant and susceptible individuals. Scottish Blackface lambs, with divergent phenotypes for resistance, were challenged with 30,000 Teladorsagia circumcincta larvae (L3), and abomasal lymph nodes recovered at 7 and 14days post-infection (dpi). High-throughput sequencing of cDNA from the abomasal lymph node was used to quantitatively sample the transcriptome with an average of 32 million reads per sample. A total of 194 and 144 genes were differentially expressed between resistant and susceptible lambs at 7 and 14 dpi respectively. Differentially expressed networks and biological processes were identified using Ingenuity Pathway Analysis. Genes involved in the inflammatory response, attraction of T lymphocytes and binding of leukocytes were more highly expressed in resistant animals at 7 dpi and in susceptible animals at 14 dpi indicating that resistant animals respond to infection earlier than susceptible animals. Twenty-four Single Nucleotide Polymorphisms (SNP) within 11 differentially expressed genes, were tested for association with gastrointestinal nematode resistance in the Scottish Blackface lambs. Four SNP, in 2 genes (SLC30A2 and ALB), were suggestively associated with faecal egg count. In conclusion, a large number of genes were differentially expressed in the abomasal lymph node of resistant and susceptible lambs responding to gastrointestinal nematode challenge. Resistant Scottish Blackface lambs appear to generate an earlier immune response to T. circumcincta. In susceptible lambs this response appears to be delayed. SNP in 2 differentially expressed genes were suggestively associated with faecal egg count indicating that differentially expressed genes may be considered candidate loci for mediating nematode resistance.

Genome-wide association study for endocrine fertility traits using single nucleotide polymorphism arrays and sequence variants in dairy cattle

Endocrine fertility traits, which are defined from progesterone concentration levels in milk, are interesting indicators of dairy cow fertility because they more directly reflect the cows own reproductive physiology than classical fertility traits, which are more biased by farm management decisions. The aim of this study was to detect quantitative trait loci (QTL) for 7 endocrine fertility traits in dairy cows by performing a genome-wide association study with 85k single nucleotide polymorphisms (SNP), and then fine-map targeted QTL regions, using imputed sequence variants. Two classical fertility traits were also analyzed for QTL with 85k SNP. The association between a SNP and a phenotype was assessed by single-locus regression for each SNP, using a linear mixed model that included a random polygenic effect. A total of 2,447 Holstein Friesian cows with 5,339 lactations with both phenotypes and genotypes were used for association analysis. Heritability estimates ranged from 0.09 to 0.15 for endocrine fertility traits and 0.03 to 0.10 for classical fertility traits. The genome-wide association study identified 17 QTL regions for endocrine fertility traits on Bos taurus autosomes (BTA) 2, 3, 8, 12, 15, 17, 23, and 25. The highest number (5) of QTL regions from the genome-wide association study was identified for the endocrine trait "proportion of samples with luteal activity." Overlapping QTL regions were found between endocrine traits on BTA 2, 3, and 17. For the classical trait calving to first service, 3 QTL regions were identified on BTA 3, 15, and 23, and an overlapping region was identified on BTA 23 with endocrine traits. Fine-mapping target regions for the endocrine traits on BTA 2 and 3 using imputed sequence variants confirmed the QTL from the genome-wide association study, and identified several associated variants that can contribute to an index of markers for genetic improvement of fertility. Several potential candidate genes underlying endocrine fertility traits were also identified in the target regions and are discussed. However, due to high linkage disequilibrium, it was not possible to specify genes or polymorphisms as causal factors for any of the regions.

Immunity to Haemonchus contortus and Vaccine Development

Sheep are capable of developing protective immunity to Haemonchus contortus through repeated exposure to this parasite, although this immune protection is the result of a complex interaction among age, gender, physiological status, pregnancy, lactation, nutrition and innate and adaptive immunity in the host animal. There are multiple effectors of the protective immune response, which differ depending on the developmental stage of the parasite being targeted, and our understanding of the effector mechanisms has developed considerably in the 2000s. The rational design of vaccines based on 'natural' or 'exposed' antigens depends on an understanding of this exposure-induced immunity. However, the most effective current vaccines rely on protection via the induction of high circulating antibody levels to 'hidden' gut antigens of H. contortus. The success of this latter strategy has resulted in the launch of a vaccine, which is based on extracts of the parasite's gut, to aid in the control of Haemonchus in Australia. The development of recombinant subunit vaccines based on the components of the successful native vaccine has not yet been achieved and most of the recent successes with recombinant subunit vaccines have focussed on antigens unrelated to the gut antigens. The future integration of an understanding of the immunobiology of this parasite with advances in antigen identification, expression (or synthesis) and presentation is likely to be pivotal to the further development of these recombinant subunit vaccines. Recent progress in each of the components underpinning this integrated approach is summarized in this review.

Detection and replication of QTL underlying resistance to gastrointestinal nematodes in adult sheep using the ovine 50K SNP array

BACKGROUND

Persistence of gastrointestinal nematode (GIN) infection and the related control methods have major impacts on the sheep industry worldwide. Based on the information generated with the Illumina OvineSNP50 BeadChip (50 K chip), this study aims at confirming quantitative trait loci (QTL) that were previously identified by microsatellite-based genome scans and identifying new QTL and allelic variants that are associated with indicator traits of parasite resistance in adult sheep. We used a commercial half-sib population of 518 Spanish Churra ewes with available data for fecal egg counts (FEC) and serum levels of immunoglobulin A (IgA) to perform different genome scan QTL mapping analyses based on classical linkage analysis (LA), a combined linkage disequilibrium and linkage analysis (LDLA) and a genome-wide association study (GWAS).

RESULTS

For the FEC and IgA traits, we detected a total of three 5 % chromosome-wise significant QTL by LA and 63 significant regions by LDLA, of which 13 reached the 5 % genome-wise significance level. The GWAS also revealed 10 significant SNPs associated with IgAt, although no significant associations were found for LFEC. Some of the significant QTL for LFEC that were detected by LA and LDLA on OAR6 overlapped with a highly significant QTL that was previously detected in a different half-sib population of Churra sheep. In addition, several new QTL and SNP associations were identified, some of which show correspondence with effects that were reported for different populations of young sheep. Other significant associations that did not coincide with previously reported associations could be related to the specific immune response of adult animals.

DISCUSSION

Our results replicate a FEC-related QTL located on OAR6 that was previously reported in Churra sheep and provide support for future research on the identification of the allelic variant that underlies this QTL. The small proportion of genetic variance explained by the detected QTL and the large number of functional candidate genes identified here are consistent with the hypothesis that GIN resistance/susceptibility is a complex trait that is not determined by individual genes acting alone but rather by complex multi-gene interactions. Future studies that combine genomic variation analysis and functional genomic information may help elucidate the biology of GIN disease resistance in sheep.

Genomic prediction and genome-wide association study for dagginess and host internal parasite resistance in New Zealand sheep

Background

Dagginess (faecal soiling of the perineum region) and host nematode parasite resistance are important animal welfare traits in New Zealand sheep. Genomic prediction (GP) estimates the genetic merit, as a molecular breeding value (mBV), for each trait based on many SNPs. The additional information the mBV provides (as determined by its accuracy) has led to its incorporation into breeding schemes. Some GP methods give SNP effects, which provide additional information to identify genome-wide associations (GWAS) for a trait of interest. Here we report results from a GP and GWAS study for dagginess and host nematode parasite resistance in a New Zealand sheep industry resource.

Results

Genomic prediction analysis was performed using 50K SNP chip data and parent average-removed, de-regressed BVs for five traits, from a resource of 8705 pedigree recorded animals. The five traits were dag score at three and eight months (DAG3, DAG8) and nematode faecal egg count in summer (FEC1), autumn (FEC2) and as an adult (AFEC). The resource consisted of Romney, Coopworth, Perendale, Texel and various breed crosses (designated: CompRCP, CompRCPT and CompCRP). The pure breeds, apart from Texel, plus CompRCP were used to develop the GP. The resulting SNP effects were used to identify genetic regions associated with dagginess and parasite resistance. Accuracies of the weighted correlation between mBV and true BV ranged between −0.07 (Texel) and 0.56 (Coopworth) for DAG3 and DAG8. For FEC1, FEC2 and AFEC accuracies ranged between −0.22 (CompRCPT) and 0.69 (Coopworth). The weighted average individual accuracy (calculated from theory) ranges were 0.13 (Texel) to 0.52 (Coopworth) and 0.11 (Texel) to 0.55 (Coopworth) respectively, for dagginess and parasite traits. There was one SNP for DAG8 that reached Bonferroni significance threshold (P < 1 × 10⁻⁶) on OAR15, the same two SNPs for each of the parasite traits (OAR26) and none for DAG3. A notable peak was also observed on OAR7 for all the parasite traits, however, it did not reach the Bonferroni significance threshold.

Conclusions

This study presents the first results of a GWAS on dagginess and faecal egg count traits in New Zealand sheep. The results suggest that there are quantitative trait loci on OAR 15 for dagginess and on OAR26 and seven for faecal egg count.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2148-2) contains supplementary material, which is available to authorized users.

Genetic resistance to infections in sheep

This paper considers genetic resistance to infectious disease in sheep, with appropriate comparison with goats, and explores how such variation may be used to assist in disease control. Many studies have attempted to quantify the extent to which host animals differ genetically in their resistance to infection or in the disease side-effects of infection, using either recorded animal pedigrees or information from genetic markers to quantify the genetic variation. Across all livestock species, whenever studies are sufficiently well powered, then genetic variation in disease resistance is usually seen and such evidence is presented here for three infections or diseases of importance to sheep, namely mastitis, foot rot and scrapie. A further class of diseases of importance in most small ruminant production systems, gastrointestinal nematode infections, is outside the scope of this review. Existence of genetic variation implies the opportunity, at least in principle, to select animals for increased resistance, with such selection ideally used as part of an integrated control strategy. For each of the diseases under consideration, evidence for genetic variation is presented, the role of selection as an aid to disease control is outlined and possible side effects of selection in terms of effects in performance, effects on resistance to other diseases and potential parasite/pathogen coevolution risks are considered. In all cases, the conclusion is drawn that selection should work and it should be beneficial, with the main challenge being to define cost effective selection protocols that are attractive to sheep farmers.

Estimation of genetic parameters for resistance to gastro-intestinal nematodes in pure blood Arabian horses

Equine internal parasites, mostly cyathostomins, affect both horse welfare and performance. The appearance of anthelmintic-resistant parasites creates a pressing need for optimising drenching schemes. This optimization may be achieved by identifying genetic markers associated with host susceptibility to infection and then to drench carriers of these markers. The aim of our study was to characterise the genetics of horse resistance to strongyle infection by estimating heritability of this trait in an Arabian pure blood population. A population of 789 Arabian pure blood horses from the Michałów stud farm, Poland were measured for strongyle egg excretion twice a year, over 8 years. Low repeatability values were found for faecal egg counts. Our analyses showed that less than 10% of the observed variation for strongyle faecal egg counts in this population had a genetic origin. However, additional analyses highlighted an age-dependent increase in heritability which was 0.04 (±0.02) in young horses (up to 3 years of age) but 0.21 (±0.04) in older ones. These results suggest that a significant part of the inter-individual variation has a genetic origin. This paves the way to a genomic dissection of horse-nematode interactions which might provide predictive markers of susceptibility, allowing individualised drenching schemes.

Selection methods for resistance to and tolerance of helminths in livestock

Helminthiases are among the most important livestock diseases worldwide, in particular for small ruminants, which are the focus of this review. Resource Allocation Theory implies that high-productivity farm animals proportionate insufficient resources for adequate coping with stressful conditions. Significant differences between breeds and within breeds are seen, as well as genotype vs. environment interactions. With improvement of genetic host resistance to infection, transmission of infection will be impacted. On the other hand, genetic improvement of resilience can lead to a reduction in clinical signs of disease, but not necessarily reduce transmission of infection to other animals. Faecal egg count (FEC) is the main measurement used to evaluate helminthiasis load, despite the fact that the protocols and analytical methods can affect the results, and the FEC data frequently shows aggregative, negative skewed distribution, and a high coefficient of variation. Mass selection where heritability is generally medium to low generally produces slow results and low economic returns. Many studies have been published linking resistance to nematodes in livestock to Quantitative Trait Loci and most studies have concentrated on chromosomes where the major histocompatibility complex region is located. Nevertheless, these complex traits have been seen to be affected by thousands of variants that each has a small effect. More recent studies have shown that genome-wide selection strategies can be useful in selecting animals for improved production and resistance traits in this case.

Current status for gastrointestinal nematode diagnosis in small ruminants: where are we and where are we going?

Gastrointestinal nematode (GIN) parasites pose a significant economic burden particularly in small ruminant production systems. Anthelmintic resistance is a serious concern to the effective control of GIN parasites and has fuelled the focus to design and promote sustainable control of practices of parasite control. Many facets of sustainable GIN parasite control programs rely on the ability to diagnose infection both qualitatively and quantitatively. Diagnostics are required to determine anthelmintic efficacies, for targeted treatment programs and selection of animals for parasite resistant breeding. This review describes much of the research investigated to date to improve the current diagnostic for the above practices which is based on counting the number of parasite eggs in faeces.

Genomics and disease resistance studies in livestock

This paper considers the application of genetic and genomic techniques to disease resistance, the interpretation of data arising from such studies and the utilisation of the research outcomes to breed animals for enhanced resistance. Resistance and tolerance are defined and contrasted, factors affecting the analysis and interpretation of field data presented, and appropriate experimental designs discussed. These general principles are then applied to two detailed case studies, infectious pancreatic necrosis in Atlantic salmon and bovine tuberculosis in dairy cattle, and the lessons learnt are considered in detail. It is concluded that the rate limiting step in disease genetic studies will generally be provision of adequate phenotypic data, and its interpretation, rather than the genomic resources. Lastly, the importance of cross-disciplinary dialogue between the animal health and animal genetics communities is stressed.

Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes

Background

Gastrointestinal nematodes are one of the most serious causes of disease in domestic ruminants worldwide. There is considerable variation in resistance to gastrointestinal nematodes within and between sheep breeds, which appears to be due to underlying genetic diversity. Through selection of resistant animals, rapid genetic progress has been demonstrated in both research and commercial flocks. Recent advances in genome sequencing and genomic technologies provide new opportunities to understand the ovine host response to gastrointestinal nematodes at the molecular level, and to identify polymorphisms conferring nematode resistance.

Results

Divergent lines of Romney and Perendale sheep, selectively bred for high and low faecal nematode egg count, were genotyped using the Illumina® Ovine SNP50 BeadChip. The resulting genome-wide SNP data were analysed for selective sweeps on loci associated with resistance or susceptibility to gastrointestinal nematode infection. Population differentiation using F_ST and Peddrift revealed sixteen regions, which included candidate genes involved in chitinase activity and the cytokine response. Two of the sixteen regions identified were contained within previously identified QTLs associated with nematode resistance.

Conclusions

In this study we identified fourteen novel regions associated with resistance or susceptibility to gastrointestinal nematodes. Results from this study support the hypothesis that host resistance to internal nematode parasites is likely to be controlled by a number of loci of moderate to small effects.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-637) contains supplementary material, which is available to authorized users.

Genomic variation and population structure detected by single nucleotide polymorphism arrays in Corriedale, Merino and Creole sheep

The aim of this study was to investigate the genetic diversity within and among three breeds of sheep: Corriedale, Merino and Creole. Sheep from the three breeds (Merino n = 110, Corriedale n = 108 and Creole n = 10) were genotyped using the Illumina Ovine SNP50 beadchip^®. Genetic diversity was evaluated by comparing the minor allele frequency (MAF) among breeds. Population structure and genetic differentiation were assessed using STRUCTURE software, principal component analysis (PCA) and fixation index (F_ST). Fixed markers (MAF = 0) that were different among breeds were identified as specific breed markers. Using a subset of 18,181 single nucleotide polymorphisms (SNPs), PCA and STUCTURE analysis were able to explain population stratification within breeds. Merino and Corriedale divergent lines showed high levels of polymorphism (89.4% and 86% of polymorphic SNPs, respectively) and moderate genetic differentiation (F_ST = 0.08) between them. In contrast, Creole had only 69% polymorphic SNPs and showed greater genetic differentiation from the other two breeds (F_ST = 0.17 for both breeds). Hence, a subset of molecular markers present in the OvineSNP50 is informative enough for breed assignment and population structure analysis of commercial and Creole breeds.

Natural selection on individual variation in tolerance of gastrointestinal nematode infection

A 25-year study of wild sheep shows that individuals vary in how quickly they lose weight as parasite infections increase, and that those who lose the least weight when heavily infected produce more offspring.

Hosts may mitigate the impact of parasites by two broad strategies: resistance, which limits parasite burden, and tolerance, which limits the fitness or health cost of increasing parasite burden. The degree and causes of variation in both resistance and tolerance are expected to influence host–parasite evolutionary and epidemiological dynamics and inform disease management, yet very little empirical work has addressed tolerance in wild vertebrates. Here, we applied random regression models to longitudinal data from an unmanaged population of Soay sheep to estimate individual tolerance, defined as the rate of decline in body weight with increasing burden of highly prevalent gastrointestinal nematode parasites. On average, individuals lost weight as parasite burden increased, but whereas some lost weight slowly as burden increased (exhibiting high tolerance), other individuals lost weight significantly more rapidly (exhibiting low tolerance). We then investigated associations between tolerance and fitness using selection gradients that accounted for selection on correlated traits, including body weight. We found evidence for positive phenotypic selection on tolerance: on average, individuals who lost weight more slowly with increasing parasite burden had higher lifetime breeding success. This variation did not have an additive genetic basis. These results reveal that selection on tolerance operates under natural conditions. They also support theoretical predictions for the erosion of additive genetic variance of traits under strong directional selection and fixation of genes conferring tolerance. Our findings provide the first evidence of selection on individual tolerance of infection in animals and suggest practical applications in animal and human disease management in the face of highly prevalent parasites.

Animals can defend themselves against parasites through either resistance (reducing parasite numbers, for example, by killing them) or tolerance (maintaining health as infections levels increase, for example, by repairing damage). Resistance has been well-studied in wild animals, but tolerance has been less so. We analysed data on body weight collected over 25 years on a natural population of Soay sheep, infected with parasitic gut worms. As parasite burden increased, sheep lost weight. Crucially, there was variation among individuals: some lost weight rapidly with increasing infections (i.e., showed “low tolerance”), whereas others lost weight slowly (i.e., showed “high tolerance”). The least tolerant individuals lost 4.5 kg of body weight across the range of parasite burdens that we saw, whereas the most tolerant lost only around 0.36 kg. However, variation in tolerance did not have a heritable genetic basis, so that although tolerance varied between individuals, this was not due to genetic differences. Further analysis revealed that there was natural selection on tolerance. Individuals who were more tolerant of infection produced more offspring over the course of their lives. This study shows that natural selection can act upon resistance and tolerance simultaneously in nature, a result that has implications for both human health and livestock management.

Functional investigation of a QTL affecting resistance to Haemonchus contortus in sheep

This study reports a functional characterization of a limited segment (QTL) of sheep chromosome 12 associated with resistance to the abomasal nematode Haemonchus contortus. The first objective was to validate the identified QTL through the comparison of genetically susceptible (N) and resistant (R) sheep produced from Martinik × Romane back-cross sheep. The R and N genotype groups were then experimentally infected with 10 000 H. contortus larvae and measured for FEC (every three days from 18 to 30 days post-challenge), haematocrit, worm burden and fertility. Significant differences in FEC and haematocrit drop were found between R and N sheep. In addition, the female worms recovered from R sheep were less fecund. The second step of the characterization was to investigate functional mechanisms associated with the QTL, thanks to a gene expression analysis performed on the abomasal mucosa and the abomasal lymph node. The gene expression level of a candidate gene lying within the QTL region (PAPP-A2) was measured. In addition, putative interactions between the chromosome segment under study and the top ten differentially expressed genes between resistant MBB and susceptible RMN sheep highlighted in a previous microarray experiment were investigated. We found an induction of Th-2 related cytokine genes expression in the abomasal mucosa of R sheep. Down-regulation of the PAPP-A2 gene expression was observed between naïve and challenged sheep although no differential expression was recorded between challenged R and N sheep. The genotyping of this limited region should contribute to the ability to predict the intrinsic resistance level of sheep.

Review: Genetics of helminth resistance in sheep

Karrow, N. A., Goliboski, K., Stonos, N., Schenkel, F. and Peregrine, A. 2014. Review: Genetics of helminth resistance in sheep. Can. J. Anim. Sci. 94: 1–9. Gastrointestinal helminth parasites are an important source of economic loss to sheep producers. A rapid increase in anthelmintic resistance has occurred around the globe; therefore, the industry is exploring alternative strategies such as genetic selection to control losses attributed to helminth infection. Since helminths have co-evolved with sheep for millions of years, natural selection for enhanced helminth resistance has occurred within certain breeds from various parts of the world. These breeds of sheep are being used to better understand the genetic aspects of helminth resistance. If the genetic variants that contribute to this phenotype can be identified, it may be possible to use selection strategies to introduce resistance alleles into other breeds or to increase their frequency within breeds. This review will provide an up-to-date overview of the pathology of helminth disease, the immune response to helminth infection, and the search for genes that confer helminth resistance.

Functional study of a genetic marker allele associated with resistance to Ascaris suum in pigs

Two single nucleotide polymorphisms (SNP TXNIP and SNP ARNT), both on chromosome 4, have been reported to be associated with roundworm (Ascaris suum) burden in pigs. In the present study, we selected pigs with two SNP TXNIP genotypes (AA; n = 24 and AB; n = 24), trickle-infected them with A. suum from 8 weeks of age until necropsy 8 weeks later, and tested the hypothesis that pigs with the AA genotype would have higher levels of resistance than pigs of AB genotype. We used different indicators of resistance (worm burden, fecal egg counts (FEC), number of liver white spots and A. suum-specific serum IgG antibody levels). Pigs of the AA genotype had lower mean macroscopic worm burden (2·4 vs 19·3; P = 0·06), lower mean total worm burden (26·5 vs 70·1; P = 0·09) and excreted fewer A. suum eggs at week 8 PI (mean number of eggs/g feces: 238 vs 1259; P = 0·14) than pigs of the AB genotype, as expected based on prior associations. The pigs were also genotyped at another locus (SNP ARNT) which showed a similar trend. This study provides suggestive evidence that resistant pigs may be selected using a genetic marker, TXNIP, and provides further support to the quantitative trait locus on chromosome 4.

A joint analysis to identify loci underlying variation in nematode resistance in three European sheep populations

Gastrointestinal nematode infections are one of the main health/economic issues in sheep industries, worldwide. Indicator traits for resistance such as faecal egg count (FEC) are commonly used in genomic studies; however, published results are inconsistent among breeds. Meta (or joint)-analysis is a tool for aggregating information from multiple independent studies. The aim of this study was to identify loci underlying variation in FEC, as an indicator of nematode resistance, in a joint analysis using data from three populations (Scottish Blackface, Sarda × Lacaune and Martinik Black-Belly × Romane), genotyped with the ovine 50k SNP chip. The trait analysed was the average animal effect for Strongyles and Nematodirus FEC data. Analyses were performed with regional heritability mapping (RHM), fitting polygenic effects with either the whole genomic relationship matrix or matrices excluding the chromosome being interrogated. Across-population genomic covariances were set to zero. After quality control, 4123 animals and 38 991 SNPs were available for the analysis. RHM identified genome-wide significant regions on OAR4, 12, 14, 19 and 20, with the latter being the most significant. The OAR20 region is close to the major histocompatibility complex, which has often been proposed as a functional candidate for nematode resistance. This region was significant only in the Sarda × Lacaune population. Several other regions, on OAR1, 3, 4, 5, 7, 12, 19, 20 and 24, were significant at the suggestive level.

Can we use genetic and genomic approaches to identify candidate animals for targeted selective treatment

Estimated breeding values (EBV) for faecal egg count (FEC) and genetic markers for host resistance to nematodes may be used to identify resistant animals for selective breeding programmes. Similarly, targeted selective treatment (TST) requires the ability to identify the animals that will benefit most from anthelmintic treatment. A mathematical model was used to combine the concepts and evaluate the potential of using genetic-based methods to identify animals for a TST regime. EBVs obtained by genomic prediction were predicted to be the best determinant criterion for TST in terms of the impact on average empty body weight and average FEC, whereas pedigree-based EBVs for FEC were predicted to be marginally worse than using phenotypic FEC as a determinant criterion. Whilst each method has financial implications, if the identification of host resistance is incorporated into a wider genomic selection indices or selective breeding programmes, then genetic or genomic information may be plausibly included in TST regimes.

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.

Toward genomic prediction from whole-genome sequence data: impact of sequencing design on genotype imputation and accuracy of predictions

Genomic prediction from whole-genome sequence data is attractive, as the accuracy of genomic prediction is no longer bounded by extent of linkage disequilibrium between DNA markers and causal mutations affecting the trait, given the causal mutations are in the data set. A cost-effective strategy could be to sequence a small proportion of the population, and impute sequence data to the rest of the reference population. Here, we describe strategies for selecting individuals for sequencing, based on either pedigree relationships or haplotype diversity. Performance of these strategies (number of variants detected and accuracy of imputation) were evaluated in sequence data simulated through a real Belgian Blue cattle pedigree. A strategy (AHAP), which selected a subset of individuals for sequencing that maximized the number of unique haplotypes (from single-nucleotide polymorphism panel data) sequenced gave good performance across a range of variant minor allele frequencies. We then investigated the optimum number of individuals to sequence by fold coverage given a maximum total sequencing effort. At 600 total fold coverage (x 600), the optimum strategy was to sequence 75 individuals at eightfold coverage. Finally, we investigated the accuracy of genomic predictions that could be achieved. The advantage of using imputed sequence data compared with dense SNP array genotypes was highly dependent on the allele frequency spectrum of the causative mutations affecting the trait. When this followed a neutral distribution, the advantage of the imputed sequence data was small; however, when the causal mutations all had low minor allele frequencies, using the sequence data improved the accuracy of genomic prediction by up to 30%.