Posterior-based proposals for speeding up Markov chain Monte Carlo.

Markov chain Monte Carlo (MCMC) is widely used for Bayesian inference in models of complex systems. Performance, however, is often unsatisfactory in models with many latent variables due to so-called poor mixing, necessitating the development of application-specific implementations. This paper introduces 'posterior-based proposals' (PBPs), a new type of MCMC update applicable to a huge class of statistical models (whose conditional dependence structures are represented by directed acyclic graphs). PBPs generate large joint updates in parameter and latent variable space, while retaining good acceptance rates (typically 33%). Evaluation against other approaches (from standard Gibbs/random walk updates to state-of-the-art Hamiltonian and particle MCMC methods) was carried out for widely varying model types: an individual-based model for disease diagnostic test data, a financial stochastic volatility model, a mixed model used in statistical genetics and a population model used in ecology. While different methods worked better or worse in different scenarios, PBPs were found to be either near to the fastest or significantly faster than the next best approach (by up to a factor of 10). PBPs, therefore, represent an additional general purpose technique that can be usefully applied in a wide variety of contexts.

An analysis of effects of heterozygosity in dairy cattle for bovine tuberculosis resistance.

Genetic selection of cattle more resistant to bovine tuberculosis (bTB) may offer a complementary control strategy. Hypothesising underlying non-additive genetic variation, we present an approach using genome-wide high density markers to identify genomic loci with dominance effects on bTB resistance and to test previously published regions with heterozygote advantage in bTB. Our data comprised 1151 Holstein-Friesian cows from Northern Ireland, confirmed bTB cases and controls, genotyped with the 700K Illumina BeadChip. Genome-wide markers were tested for associations between heterozygosity and bTB status using marker-based relationships. Results were tested for robustness against genetic structure, and the genotypic frequencies of a significant locus were tested for departures from Hardy-Weinberg equilibrium. Genomic regions identified in our study and in previous publications were tested for dominance effects. Genotypic effects were estimated through ASReml mixed models. A SNP (rs43032684) on chromosome 6 was significant at the chromosome-wide level, explaining 1.7% of the phenotypic variance. In the controls, there were fewer heterozygotes for rs43032684 (P < 0.01) with the genotypic values suggesting that heterozygosity confers a heterozygote disadvantage. The region surrounding rs43032684 had a significant dominance effect (P < 0.01). SNP rs43032684 resides within a pseudogene with a parental gene involved in macrophage response to infection and within a copy-number-variation region previously associated with nematode resistance. No dominance effect was found for the region on chromosome 11, as indicated by a previous candidate region bTB study. These findings require further validation with large-scale data.

Fine-mapping host genetic variation underlying outcomes to Mycobacterium bovis infection in dairy cows.

BACKGROUND: Susceptibility to Mycobacterium bovis infection in cattle is governed in part by host genetics. However, cattle diagnosed as infected with M. bovis display varying signs of pathology. The variation in host response to infection could represent a continuum since time of exposure or distinct outcomes due to differing pathogen handling. The relationships between host genetics and variation in host response and pathological sequelae following M. bovis infection were explored by genotyping 1966 Holstein-Friesian dairy cows at 538,231 SNPs with three distinct phenotypes. These were: single intradermal cervical comparative tuberculin (SICCT) test positives with visible lesions (VLs), SICCT-positives with undetected visible lesions (NVLs) and matched controls SICCT-negative on multiple occasions. RESULTS: Regional heritability mapping identified three loci associated with the NVL phenotype on chromosomes 17, 22 and 23, distinct to the region on chromosome 13 associated with the VL phenotype. The region on chromosome 23 was at genome-wide significance and candidate genes overlapping the mapped window included members of the bovine leukocyte antigen class IIb region, a complex known for its role in immunity and disease resistance. Chromosome heritability analysis attributed variance to six and thirteen chromosomes for the VL and NVL phenotypes, respectively, and four of these chromosomes were found to explain a proportion of the phenotypic variation for both the VL and NVL phenotype. By grouping the M. bovis outcomes (VLs and NVLs) variance was attributed to nine chromosomes. When contrasting the two M. bovis infection outcomes (VLs vs NVLs) nine chromosomes were found to harbour heritable variation. Regardless of the case phenotype under investigation, chromosome heritability did not exceed 8% indicating that the genetic control of bTB resistance consists of variants of small to moderate effect situated across many chromosomes of the bovine genome. CONCLUSIONS: These findings suggest the host genetics of M. bovis infection outcomes is governed by distinct and overlapping genetic variants. Thus, variation in the pathology of M. bovis infected cattle may be partly genetically determined and indicative of different host responses or pathogen handling. There may be at least three distinct outcomes following M. bovis exposure in dairy cattle: resistance to infection, infection resulting in pathology or no detectable pathology.

Genetic evaluation for bovine tuberculosis resistance in dairy cattle.

Genetic evaluations for resistance to bovine tuberculosis (bTB) were calculated based on British national data including individual animal tuberculin skin test results, postmortem examination (presence of bTB lesions and bacteriological culture for Mycobacterium bovis), animal movement and location information, production history, and pedigree records. Holstein cows with identified sires in herds with bTB breakdowns (new herd incidents) occurring between the years 2000 and 2014 were considered. In the first instance, cows with a positive reaction to the skin test and a positive postmortem examination were defined as infected. Values of 0 and 1 were assigned to healthy and infected animal records, respectively. Data were analyzed with mixed models. Linear and logit function heritability estimates were 0.092 and 0.172, respectively. In subsequent analyses, breakdowns were split into 2-mo intervals to better model time of exposure and infection in the contemporary group. Intervals with at least one infected individual were retained and multiple intervals within the same breakdown were included. Healthy animal records were assigned values of 0, and infected records a value of 1 in the interval of infection and values reflecting a diminishing probability of infection in the preceding intervals. Heritability and repeatability estimates were 0.115 and 0.699, respectively. Reliabilities and across time stability of the genetic evaluation were improved with the interval model. Subsequently, 2 more definitions of "infected" were analyzed with the interval model: (1) all positive skin test reactors regardless of postmortem examination, and (2) all positive skin test reactors plus nonreactors with positive postmortem examination. Estimated heritability was 0.085 and 0.089, respectively; corresponding repeatability estimates were 0.701 and 0.697. Genetic evaluation reliabilities and across time stability did not change. Correlations of genetic evaluations for bTB with other traits in the current breeding goal were mostly not different from zero. Correlation with the UK Profitable Lifetime Index was moderate, significant, and favorable. Results demonstrated the feasibility of a national genetic evaluation for bTB resistance. Selection for enhanced resistance will have a positive effect on profitability and no antagonistic effects on current breeding goal traits. Official genetic evaluations are now based on the interval model and the last bTB trait definition.

Modelling gastrointestinal parasitism infection in a sheep flock over two reproductive seasons: in silico exploration and sensitivity analysis.

In reproducing ewes, a periparturient breakdown of immunity is often observed to result in increased fecal egg excretion, making them the main source of infection for their immunologically naive lambs. In this study, we expanded a simulation model previously developed for growing lambs to explore the impact of the genotype (performance and resistance traits) and host nutrition on the performance and parasitism of both growing lambs and reproducing ewes naturally infected with Teladorsagia circumcincta. Our model accounted for nutrient-demanding phases, such as gestation and lactation, and included a supplementary module to manage the age structure of the ewe flock. The model was validated by comparison with published data. Because model parameters were unknown or poorly estimated, detailed sensitivity analysis of the model was performed for the sheep mortality and the level of infection, following a preliminary screening step. The parameters with the greatest effect on parasite-related outputs were those driving animal growth and milk yield. Our model enables different parasite-control strategies (host nutrition, breeding for resistance and anthelmintic treatments) to be assessed on the long term in a sheep flock. To optimize in silico exploration, the parameters highlighted by the sensitivity analysis should be refined with real data.

Genome-wide association study of tick resistance in South African Nguni cattle.

Ticks and tick-borne diseases are among the main causes of economic loss in the South African cattle industry through high morbidity and mortality rates. Concerns of the general public regarding chemical residues may tarnish their perceptions of food safety and environmental health when the husbandry of cattle includes frequent use of acaricides to manage ticks. The primary objective of this study was to identify single nucleotide polymorphism (SNP) markers associated with host resistance to ticks in South African Nguni cattle. Tick count data were collected monthly from 586 Nguni cattle reared in four herds under natural grazing conditions over a period of two years. The counts were recorded for six species of ticks attached in eight anatomical locations on the animals and were summed by species and anatomical location. This gave rise to 63 measured phenotypes or traits, with results for 12 of these traits being reported here. Tick count (x) data were transformed using log10(x+1) and the resulting values were examined for normality. DNA was extracted from hair and blood samples and was genotyped using the Illumina BovineSNP50 assay. After quality control (call rate >90%, minor allele frequency >0.02), 40,436 SNPs were retained for analysis. Genetic parameters were estimated and association analysis for tick resistance was carried out using two approaches: a genome-wide association (GWA) analysis using the GenABEL package and a regional heritability mapping (RHM) analysis. The Bonferroni genome-wide (P<0.05) corrected significance threshold was 1.24×10(-6), with 2.47×10(-5) as the suggestive significance threshold (P<0.10) (i.e., one false positive per genome scan) in the GWA analysis. Likelihood ratio test (LRT) thresholds for genome-wide and suggestive significance were 13.5 and 9.15 for the RHM analysis. Six ixodid tick species were identified, with Amblyomma hebraeum (the vector for Heartwater disease) being the dominant species. Heritability estimates (h(2)) from the fitted animal and sire models ranged from 0.02±0.00 to 0.17±0.04 for the transformed tick count data. Several genomic regions harbouring quantitative trait loci (QTL) were identified for different tick count traits by both the GWA and RHM approaches. Three genome-wide significant regions on chromosomes 7, 10 and 19 were identified for total tick count on the head, total body A. hebraeum tick count and total A. hebraeum on the perineum region, respectively. Additional regions significant at the suggestive level were identified on chromosomes 1, 3, 6, 7, 8, 10, 11, 12, 14, 15, 17, 19 and 26 for several of the traits. The GWA approach identified more genomic regions than did the RHM approach. The chromosomal regions identified here as harbouring QTL underlying variation in tick burden form the basis for further analyses to identify specific candidate genes and polymorphisms related to cattle tick resistance and provide the potential for marker-assisted selection in Nguni cattle.

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.

Mapping and validation of a major QTL affecting resistance to pancreas disease (salmonid alphavirus) in Atlantic salmon (Salmo salar).

Pancreas disease (PD), caused by a salmonid alphavirus (SAV), has a large negative economic and animal welfare impact on Atlantic salmon aquaculture. Evidence for genetic variation in host resistance to this disease has been reported, suggesting that selective breeding may potentially form an important component of disease control. The aim of this study was to explore the genetic architecture of resistance to PD, using survival data collected from two unrelated populations of Atlantic salmon; one challenged with SAV as fry in freshwater (POP 1) and one challenged with SAV as post-smolts in sea water (POP 2). Analyses of the binary survival data revealed a moderate-to-high heritability for host resistance to PD in both populations (fry POP 1 h(2)~0.5; post-smolt POP 2 h(2)~0.4). Subsets of both populations were genotyped for single nucleotide polymorphism markers, and six putative resistance quantitative trait loci (QTL) were identified. One of these QTL was mapped to the same location on chromosome 3 in both populations, reaching chromosome-wide significance in both the sire- and dam-based analyses in POP 1, and genome-wide significance in a combined analysis in POP 2. This independently verified QTL explains a significant proportion of host genetic variation in resistance to PD in both populations, suggesting a common underlying mechanism for genetic resistance across lifecycle stages. Markers associated with this QTL are being incorporated into selective breeding programs to improve PD resistance.

Using model-based proposals for fast parameter inference on discrete state space, continuous-time Markov processes.

Bayesian statistics provides a framework for the integration of dynamic models with incomplete data to enable inference of model parameters and unobserved aspects of the system under study. An important class of dynamic models is discrete state space, continuous-time Markov processes (DCTMPs). Simulated via the Doob-Gillespie algorithm, these have been used to model systems ranging from chemistry to ecology to epidemiology. A new type of proposal, termed 'model-based proposal' (MBP), is developed for the efficient implementation of Bayesian inference in DCTMPs using Markov chain Monte Carlo (MCMC). This new method, which in principle can be applied to any DCTMP, is compared (using simple epidemiological SIS and SIR models as easy to follow exemplars) to a standard MCMC approach and a recently proposed particle MCMC (PMCMC) technique. When measurements are made on a single-state variable (e.g. the number of infected individuals in a population during an epidemic), model-based proposal MCMC (MBP-MCMC) is marginally faster than PMCMC (by a factor of 2-8 for the tests performed), and significantly faster than the standard MCMC scheme (by a factor of 400 at least). However, when model complexity increases and measurements are made on more than one state variable (e.g. simultaneously on the number of infected individuals in spatially separated subpopulations), MBP-MCMC is significantly faster than PMCMC (more than 100-fold for just four subpopulations) and this difference becomes increasingly large.

Genetic analysis of reproductive traits and antibody response in a PRRS outbreak herd.

Porcine reproductive and respiratory syndrome (PRRS) is the most economically significant disease impacting pig production in North America, Europe, and Asia, causing reproductive losses such as increased rates of stillbirth and mummified piglets. The objective of this study was to explore the genetic basis of host response to the PRRS virus (PRRSV) in a commercial multiplier sow herd before and after a PRRS outbreak, using antibody response and reproductive traits. Reproductive data comprising number born alive (NBA), number alive at 24 h (NA24), number stillborn (NSB), number born mummified (NBM), proportion born dead (PBD), number born dead (NBD), number weaned (NW), and number of mortalities through weaning (MW) of 5,227 litters from 1,967 purebred Landrace sows were used along with a pedigree comprising 2,995 pigs. The PRRS outbreak date was estimated from rolling averages of farrowing traits and was used to split the data into a pre-PRRS phase and a PRRS phase. All 641 sows in the herd during the outbreak were blood sampled 46 d after the estimated outbreak date and were tested for anti-PRRSV IgG using ELISA (sample-to-positive [S/P] ratio). Genetic parameters of traits were estimated separately for the pre-PRRS and PRRS phase data sets. Sows were genotyped using the PorcineSNP60 BeadChip, and genome-wide association studies (GWAS) were performed using method Bayes B. Heritability estimates for reproductive traits ranged from 0.01 (NBM) to 0.12 (NSB) and from 0.01 (MW) to 0.12 (NBD) for the pre-PRRS and PRRS phases, respectively. S/P ratio had heritability (0.45) and strong genetic correlations with most traits, ranging from -0.72 (NBM) to 0.73 (NBA). In the pre-PRRS phase, regions associated with NSB and PBD explained 1.6% and 3% of the genetic variance, respectively. In the PRRS phase, regions associated with NBD, NSB, and S/P ratio explained 0.8%, 11%, and 50.6% of the genetic variance, respectively. For S/P ratio, 2 regions on SSC 7 (SSC7) separated by 100 Mb explained 40% of the genetic variation, including a region encompassing the major histocompatibility complex, which explained 25% of the genetic variance. These results indicate a significant genomic component associated with PRRSV antibody response and NSB in this data set. Also, the high heritability and genetic correlation estimates for S/P ratio during the PRRS phase suggest that S/P ratio could be used as an indicator of the impact of PRRS on reproductive traits.

Effect and mode of action of the Texel muscling QTL (TM-QTL) on carcass traits in purebred Texel lambs.

TM-QTL is a quantitative trait locus (QTL) on ovine chromosome 18 (OAR18) known to affect loin muscling in Texel sheep. Previous work suggested that its mode of inheritance is consistent with paternal polar overdominance, but this has yet to be formally demonstrated. This study used purebred Texel sheep segregating for TM-QTL to confirm its presence in the chromosomal region in which it was first reported and to determine its pattern of inheritance. To do so, this study used the first available data from a Texel flock, which included homozygote TM-QTL carriers (TM/TM; n=34) in addition to homozygote non-carriers (+/+; n=40 and, heterozygote TM-QTL-carriers inheriting TM-QTL from their sire (TM/+; n=53) or their dam (+/TM; n=17). Phenotypes included a wide range of loin muscling, carcass composition and tissue distribution traits. The presence of a QTL affecting ultrasound muscle depth on OAR18 was confirmed with a paternal QTL effect ranging from +0.54 to +2.82 mm UMD (s.e. 0.37 to 0.57 mm) across the sires segregating for TM-QTL. Loin muscle width, depth and area, loin muscle volume and dissected M. longissimus lumborum weight were significantly greater for TM/+ than +/+ lambs (+2.9% to +7.9%; P<0.05). There was significant evidence that the effect of TM-QTL on the various loin muscling traits measured was paternally polar overdominant (P<0.05). In contrast, there was an additive effect of TM-QTL on both live weight at 20 weeks and carcass weight; TM/TM animals were significantly (P<0.05) heavier than +/+ (+11.1% and +7.3%, respectively) and +/TM animals (+11.9% and +11.7%, respectively), with TM/+ intermediate. Weights of the leg, saddle and shoulder region (corrected for carcass weight) were similar in the genotypic groups. There was a tendency for lambs inheriting TM-QTL from their sire to be less fat with slightly more muscle than non-carriers. For example, carcass muscle weight measured by live animal CT-scanning was 2.8% higher in TM/TM than +/+ lambs (P<0.05), carcass muscle weight measured by carcass CT-scanning was 1.36% higher in TM/+ than +/+ lambs (P<0.05), and weight of fat trimmed from the carcass cuts was significantly lower for TM/+ than +/+ lambs (-11.2%; P<0.05). No negative effects of TM-QTL on carcass traits were found. Optimal commercial use of TM-QTL within the sheep industry would require some consideration, due to the apparently different mode of action of the two main effects of TM-QTL (on growth and muscling).

Accuracy of genomic prediction within and across populations for nematode resistance and body weight traits in sheep.

Genomic prediction utilizes single nucleotide polymorphism (SNP) chip data to predict animal genetic merit. It has the advantage of potentially capturing the effects of the majority of loci that contribute to genetic variation in a trait, even when the effects of the individual loci are very small. To implement genomic prediction, marker effects are estimated with a training set, including individuals with marker genotypes and trait phenotypes; subsequently, genomic estimated breeding values (GEBV) for any genotyped individual in the population can be calculated using the estimated marker effects. In this study, we aimed to: (i) evaluate the potential of genomic prediction to predict GEBV for nematode resistance traits and BW in sheep, within and across populations; (ii) evaluate the accuracy of these predictions through within-population cross-validation; and (iii) explore the impact of population structure on the accuracy of prediction. Four data sets comprising 752 lambs from a Scottish Blackface population, 2371 from a Sarda×Lacaune backcross population, 1000 from a Martinik Black-Belly×Romane backcross population and 64 from a British Texel population were used in this study. Traits available for the analysis were faecal egg count for Nematodirus and Strongyles and BW at different ages or as average effect, depending on the population. Moreover, immunoglobulin A was also available for the Scottish Blackface population. Results show that GEBV had moderate to good within-population predictive accuracy, whereas across-population predictions had accuracies close to zero. This can be explained by our finding that in most cases the accuracy estimates were mostly because of additive genetic relatedness between animals, rather than linkage disequilibrium between SNP and quantitative trait loci. Therefore, our results suggest that genomic prediction for nematode resistance and BW may be of value in closely related animals, but that with the current SNP chip genomic predictions are unlikely to work across breeds.

Genome-wide association study identifies novel loci associated with resistance to bovine tuberculosis.

Tuberculosis (TB) caused by Mycobacterium bovis is a re-emerging disease of livestock that is of major economic importance worldwide, as well as being a zoonotic risk. There is significant heritability for host resistance to bovine TB (bTB) in dairy cattle. To identify resistance loci for bTB, we undertook a genome-wide association study in female Holstein-Friesian cattle with 592 cases and 559 age-matched controls from case herds. Cases and controls were categorised into distinct phenotypes: skin test and lesion positive vs skin test negative on multiple occasions, respectively. These animals were genotyped with the Illumina BovineHD 700K BeadChip. Genome-wide rapid association using linear and logistic mixed models and regression (GRAMMAR), regional heritability mapping (RHM) and haplotype-sharing analysis identified two novel resistance loci that attained chromosome-wise significance, protein tyrosine phosphatase receptor T (PTPRT; P=4.8 × 10(-7)) and myosin IIIB (MYO3B; P=5.4 × 10(-6)). We estimated that 21% of the phenotypic variance in TB resistance could be explained by all of the informative single-nucleotide polymorphisms, of which the region encompassing the PTPRT gene accounted for 6.2% of the variance and a further 3.6% was associated with a putative copy number variant in MYO3B. The results from this study add to our understanding of variation in host control of infection and suggest that genetic marker-based selection for resistance to bTB has the potential to make a significant contribution to bTB control.

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.

Genetic factors controlling wool shedding in a composite Easycare sheep flock.

Historically, sheep have been selectively bred for desirable traits including wool characteristics. However, recent moves towards extensive farming and reduced farm labour have seen a renewed interest in Easycare breeds. The aim of this study was to quantify the underlying genetic architecture of wool shedding in an Easycare flock. Wool shedding scores were collected from 565 pedigreed commercial Easycare sheep from 2002 to 2010. The wool scoring system was based on a 10-point (0-9) scale, with score 0 for animals retaining full fleece and 9 for those completely shedding. DNA was sampled from 200 animals of which 48 with extreme phenotypes were genotyped using a 50-k SNP chip. Three genetic analyses were performed: heritability analysis, complex segregation analysis to test for a major gene hypothesis and a genome-wide association study to map regions in the genome affecting the trait. Phenotypes were treated as a continuous or binary variable and categories. High estimates of heritability (0.80 when treated as a continuous, 0.65-0.75 as binary and 0.75 as categories) for shedding were obtained from linear mixed model analyses. Complex segregation analysis gave similar estimates (0.80 ± 0.06) to those above with additional evidence for a major gene with dominance effects. Mixed model association analyses identified four significant (P < 0.05) SNPs. Further analyses of these four SNPs in all 200 animals revealed that one of the SNPs displayed dominance effects similar to those obtained from the complex segregation analyses. In summary, we found strong genetic control for wool shedding, demonstrated the possibility of a single putative dominant gene controlling this trait and identified four SNPs that may be in partial linkage disequilibrium with gene(s) controlling shedding.

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.

The interaction of host genetics and disease processes in chronic livestock disease: a simulation model of ovine footrot.

A stochastic, individual-based, simulation model of footrot in a flock of 200 ewes was developed that included flock demography, disease processes, host genetic variation for traits influencing infection and disease processes, and bacterial contamination of the environment. Sensitivity analyses were performed using ANOVA to examine the contribution of unknown parameters to outcome variation. The infection rate and bacterial death rate were the most significant factors determining the observed prevalence of footrot, as well as the heritability of resistance. The dominance of infection parameters in determining outcomes implies that observational data cannot be used to accurately estimate the strength of genetic control of underlying traits describing the infection process, i.e. resistance. Further work will allow us to address the potential for genetic selection to control ovine footrot.

Estimation of historical effective population size using linkage disequilibria with marker data.

Theory hypothesizes that the rate of decline in linkage disequilibrium (LD) as a function of distance between markers, measured by r(2), can be used to estimate effective population size (N(e)) and how it varies over time. The development of high-density genotyping makes feasible the application of this theory and has provided an impetus to improve predictions. This study considers the impact of several developments on the estimation of N(e) using both simulated and equine high-density single-nucleotide polymorphism data, when N(e) is assumed to be constant a priori and when it is not. In all models, estimates of N(e) were highly sensitive to thresholds imposed upon minor allele frequency (MAF) and to a priori assumptions on the expected r(2) for adjacent markers. Where constant N(e) was assumed a priori, then estimates with the lowest mean square error were obtained with MAF thresholds between 0.05 and 0.10, adjustment of r(2) for finite sample size, estimation of a [the limit for r(2) as recombination frequency (c) approaches 0] and relating N(e) to c (1 - c/2). The findings for predicting N(e) from models allowing variable N(e) were much less clear, apart from the desirability of correcting for finite sample size, and the lack of consistency in estimating recent N(e) (<7 generations) where estimates use data with large c. The theoretical conflicts over how estimation should proceed and uncertainty over where predictions might be expected to fit well suggest that the estimation of N(e) when it varies be carried out with extreme caution.

In silico exploration of the impact of pasture larvae contamination and anthelmintic treatment on genetic parameter estimates for parasite resistance in grazing sheep.

A mathematical model was developed to investigate the impact of level of Teladorsagia circumcincta larval pasture contamination and anthelmintic treatment on genetic parameter estimates for performance and resistance to parasites in sheep. Currently great variability is seen for published correlations between performance and resistance, with estimates appearing to vary with production environment. The model accounted for host genotype and parasitism in a population of lambs, incorporating heritable between-lamb variation in host-parasite interactions, with genetic independence of input growth and immunological variables. An epidemiological module was linked to the host-parasite interaction module via food intake (FI) to create a grazing scenario. The model was run for a population of lambs growing from 2 mo of age, grazing on pasture initially contaminated with 0, 1,000, 3,000, or 5,000 larvae/kg DM, and given either no anthelmintic treatment or drenched at 30-d intervals. The mean population values for FI and empty BW (EBW) decreased with increasing levels of initial larval contamination (IL(0)), with non-drenched lambs having a greater reduction than drenched ones. For non-drenched lambs the maximum mean population values for worm burden (WB) and fecal egg count (FEC) increased and occurred earlier for increasing IL(0), with values being similar for all IL(0) at the end of the simulation. Drenching was predicted to suppress WB and FEC, and cause reduced pasture contamination. The heritability of EBW for non-drenched lambs was predicted to be initially high (0.55) and decreased over time with increasing IL(0), whereas drenched lambs remained high throughout. The heritability of WB and FEC for all lambs was initially low (∼0.05) and increased with time to ∼0.25, with increasing IL(0) leading to this value being reached at faster rates. The genetic correlation between EBW and FEC was initially ∼-0.3. As time progressed the correlation tended towards 0, before becoming negative by the end of the simulation for non-drenched lambs, with increasing IL(0) leading to increasingly negative correlations. For drenched lambs, the correlation remained close to 0. This study highlights the impact of IL(0) and anthelmintic treatment on genetic parameters for resistance. Along with factors affecting performance penalties due to parasitism and time of reporting, the results give plausible causes for variation in genetic parameter estimates previously reported.

Possibilities to breed for resistance to nematode parasite infections in small ruminants in tropical production systems.

Gastrointestinal nematode parasitism is the most important disease affecting livestock production systems in developing countries, particularly small ruminant production systems. Of particular importance are infections with the strongyle Haemonchus contortus. Integrated disease control strategies are required, including improved management, nutrition and wise use of anthelmintic chemicals. Increasingly, selection of sheep or goats for improved nematode resistance is viewed as a valuable option to complement other control measures. Breeding for resistance is possible because of the existence of extensive genetic variation in resistance, both within and between breeds of sheep and goats. Such breeding schemes are most likely to be based on choice of appropriate breeds adapted to the local environmental conditions, followed by phenotypic selection for resistance. Goal and selection objective traits are likely to include performance (e.g. growth rate) under conditions of parasite challenge, faecal egg count (FEC) and measures of anaemia. With current technologies, genetic markers are likely to be too expensive and logistically difficult to incorporate into breeding schemes in tropical or developing countries. Genotype by environment interactions may be expected, particularly when comparing animals in environments that differ in the extent of parasite challenge or differ in the quality of available nutrition. However, there is no reason to expect antagonistic genetic relationships between performance and resistance, and selection indices should be readily constructed that improve both performance and resistance. If FEC is decreased, then pasture contamination should also decrease, leading to additional benefits for all sheep grazing the same pasture. Finally, breeding for nematode resistance should lead to lasting and sustained improvements in resistance or tolerance. There is no empirical evidence to suggest that nematodes will evolve rapidly in response to resistant hosts, and mathematical models based on genetic and biological principles also suggest that resistance should be sustainable.