Balancing selection at a premature stop mutation in the myostatin gene underlies a recessive leg weakness syndrome in pigs.

Balancing selection provides a plausible explanation for the maintenance of deleterious alleles at moderate frequency in livestock, including lethal recessives exhibiting heterozygous advantage in carriers. In the current study, a leg weakness syndrome causing mortality of piglets in a commercial line showed monogenic recessive inheritance, and a region on chromosome 15 associated with the syndrome was identified by homozygosity mapping. Whole genome resequencing of cases and controls identified a mutation causing a premature stop codon within exon 3 of the porcine Myostatin (MSTN) gene, similar to those causing a double-muscling phenotype observed in several mammalian species. The MSTN mutation was in Hardy-Weinberg equilibrium in the population at birth, but significantly distorted amongst animals still in the herd at 110 kg, due to an absence of homozygous mutant genotypes. In heterozygous form, the MSTN mutation was associated with a major increase in muscle depth and decrease in fat depth, suggesting that the deleterious allele was maintained at moderate frequency due to heterozygous advantage (allele frequency, q = 0.22). Knockout of the porcine MSTN by gene editing has previously been linked to problems of low piglet survival and lameness. This MSTN mutation is an example of putative balancing selection in livestock, providing a plausible explanation for the lack of disrupting MSTN mutations in pigs despite many generations of selection for lean growth.

The nedd-8 activating enzyme gene underlies genetic resistance to infectious pancreatic necrosis virus in Atlantic salmon.

Genetic resistance to infectious pancreatic necrosis virus (IPNV) in Atlantic salmon is a rare example of a trait where a single locus (QTL) explains almost all of the genetic variation. Genetic marker tests based on this QTL on salmon chromosome 26 have been widely applied in selective breeding to markedly reduce the incidence of the disease. In the current study, whole genome sequencing and functional annotation approaches were applied to characterise genes and variants in the QTL region. This was complemented by an analysis of differential expression between salmon fry of homozygous resistant and homozygous susceptible genotypes challenged with IPNV. These analyses pointed to the NEDD-8 activating enzyme 1 (nae1) gene as a putative functional candidate underlying the QTL effect. The role of nae1 in IPN resistance was further assessed via CRISPR-Cas9 knockout of the nae1 gene and chemical inhibition of the nae1 protein activity in Atlantic salmon cell lines, both of which resulted in highly significant reduction in productive IPNV replication. In contrast, CRISPR-Cas9 knockout of a candidate gene previously purported to be a cellular receptor for the virus (cdh1) did not have a major impact on productive IPNV replication. These results suggest that nae1 is the causative gene underlying the major QTL affecting resistance to IPNV in salmon, provide further evidence for the critical role of neddylation in host-pathogen interactions, and highlight the value in combining high-throughput genomics approaches with targeted genome editing to understand the genetic basis of disease resistance.

Estimating individuals' genetic and non-genetic effects underlying infectious disease transmission from temporal epidemic data.

Individuals differ widely in their contribution to the spread of infection within and across populations. Three key epidemiological host traits affect infectious disease spread: susceptibility (propensity to acquire infection), infectivity (propensity to transmit infection to others) and recoverability (propensity to recover quickly). Interventions aiming to reduce disease spread may target improvement in any one of these traits, but the necessary statistical methods for obtaining risk estimates are lacking. In this paper we introduce a novel software tool called SIRE (standing for "Susceptibility, Infectivity and Recoverability Estimation"), which allows for the first time simultaneous estimation of the genetic effect of a single nucleotide polymorphism (SNP), as well as non-genetic influences on these three unobservable host traits. SIRE implements a flexible Bayesian algorithm which accommodates a wide range of disease surveillance data comprising any combination of recorded individual infection and/or recovery times, or disease diagnostic test results. Different genetic and non-genetic regulations and data scenarios (representing realistic recording schemes) were simulated to validate SIRE and to assess their impact on the precision, accuracy and bias of parameter estimates. This analysis revealed that with few exceptions, SIRE provides unbiased, accurate parameter estimates associated with all three host traits. For most scenarios, SNP effects associated with recoverability can be estimated with highest precision, followed by susceptibility. For infectivity, many epidemics with few individuals give substantially more statistical power to identify SNP effects than the reverse. Importantly, precise estimates of SNP and other effects could be obtained even in the case of incomplete, censored and relatively infrequent measurements of individuals' infection or survival status, albeit requiring more individuals to yield equivalent precision. SIRE represents a new tool for analysing a wide range of experimental and field disease data with the aim of discovering and validating SNPs and other factors controlling infectious disease transmission.

Phenotypic and genetic variation in the response of chickens to Eimeria tenella induced coccidiosis.

BACKGROUND: Coccidiosis is a major contributor to losses in poultry production. With emerging constraints on the use of in-feed prophylactic anticoccidial drugs and the relatively high costs of effective vaccines, there are commercial incentives to breed chickens with greater resistance to this important production disease. To identify phenotypic biomarkers that are associated with the production impacts of coccidiosis, and to assess their covariance and heritability, 942 Cobb500 commercial broilers were subjected to a defined challenge with Eimeria tenella (Houghton). Three traits were measured: weight gain (WG) during the period of infection, caecal lesion score (CLS) post mortem, and the level of a serum biomarker of intestinal inflammation, i.e. circulating interleukin 10 (IL-10), measured at the height of the infection. RESULTS: Phenotypic analysis of the challenged chicken cohort revealed a significant positive correlation between CLS and IL-10, with significant negative correlations of both these traits with WG. Eigenanalysis of phenotypic covariances between measured traits revealed three distinct eigenvectors. Trait weightings of the first eigenvector, (EV1, eigenvalue = 59%), were biologically interpreted as representing a response of birds that were susceptible to infection, with low WG, high CLS and high IL-10. Similarly, the second eigenvector represented infection resilience/resistance (EV2, 22%; high WG, low CLS and high IL-10), and the third eigenvector tolerance (EV3, 19%; high WG, high CLS and low IL-10), respectively. Genome-wide association studies (GWAS) identified two SNPs that were associated with WG at the suggestive level. CONCLUSIONS: Eigenanalysis separated the phenotypic impact of a defined challenge with E. tenella on WG, caecal inflammation/pathology, and production of IL-10 into three major eigenvectors, indicating that the susceptibility-resistance axis is not a single continuous quantitative trait. The SNPs identified by the GWAS for body weight were located in close proximity to two genes that are involved in innate immunity (FAM96B and RRAD).

Gene expression comparison of resistant and susceptible Atlantic salmon fry challenged with Infectious Pancreatic Necrosis virus reveals a marked contrast in immune response.

BACKGROUND: Infectious Pancreatic Necrosis (IPN) is a highly contagious birnavirus disease of farmed salmonid fish, which often causes high levels of morbidity and mortality. A large host genetic component to resistance has been previously described for Atlantic salmon (Salmo salar L.), which mediates high mortality rates in some families and zero mortality in others. However, the molecular and immunological basis for this resistance is not yet fully known. This manuscript describes a global comparison of the gene expression profiles of resistant and susceptible Atlantic salmon fry following challenge with the IPN virus. RESULTS: Salmon fry from two IPNV-resistant and two IPNV-susceptible full sibling families were challenged with the virus and sampled at 1 day, 7 days and 20 days post-challenge. Significant viral titre was observed in both resistant and susceptible fish at all timepoints, although generally at higher levels in susceptible fish. Gene expression profiles combined with gene ontology and pathway analyses demonstrated that while a clear immune response was observed in both resistant and susceptible fish, there were striking differences between the two phenotypes. The susceptible fish showed marked up-regulation of genes related to cytokine activity and inflammatory response that evidently failed to protect against the virus. In contrast, the resistant fish demonstrated a less pronounced immune response including up-regulation of genes relating to the M2 macrophage system. CONCLUSIONS: While only the susceptible phenotype shows appreciable mortality levels, both resistant and susceptible fish can become infected with IPNV. Susceptible fish are characterized by a much larger, yet ineffective, immune response, largely related to cytokine and inflammatory systems. Resistant fish demonstrate a more moderate, putative macrophage-mediated inflammatory response, which may contribute to their survival.

Quantitative genetic analysis of the bTB diagnostic single intradermal comparative cervical test (SICCT).

BACKGROUND: Bovine tuberculosis (bTB) is a disease of significant economic importance and is a persistent animal health problem with implications for public health worldwide. Control of bTB in the UK has relied on diagnosis through the single intradermal comparative cervical test (SICCT). However, limitations in the sensitivity of this test hinder successful eradication and the control of bTB remains a major challenge. Genetic selection for cattle that are more resistant to bTB infection can assist in bTB control. The aim of this study was to conduct a quantitative genetic analysis of SICCT measurements collected during bTB herd testing. Genetic selection for bTB resistance will be partially informed by SICCT-based diagnosis; therefore it is important to know whether, in addition to increasing bTB resistance, this might also alter genetically the epidemiological characteristics of SICCT. RESULTS: Our main findings are that: (1) the SICCT test is robust at the genetic level, since its hierarchy and comparative nature provide substantial protection against random genetic changes that arise from genetic drift and from correlated responses among its components due to either natural or artificial selection; (2) the comparative nature of SICCT provides effective control for initial skin thickness and age-dependent differences; and (3) continuous variation in SICCT is only lowly heritable and has a weak correlation with SICCT positivity among healthy animals which was not significantly different from zero (P > 0.05). These emerging results demonstrate that genetic selection for bTB resistance is unlikely to change the probability of correctly identifying non-infected animals, i.e. the test's specificity, while reducing the overall number of cases. CONCLUSIONS: This study cannot exclude all theoretical risks from selection on resistance to bTB infection but the role of SICCT in disease control is unlikely to be rapidly undermined, with any adverse correlated responses expected to be weak and slow, which allow them to be monitored and managed.

Genome-wide association reveals QTL for growth, bone and in vivo carcass traits as assessed by computed tomography in Scottish Blackface lambs.

BACKGROUND: Improving meat quality including taste and tenderness is critical to the protection and development of markets for sheep meat. Phenotypic selection for such measures of meat quality is constrained by the fact that these parameters can only be measured post-slaughter. Carcass composition has an impact on meat quality and can be measured on live animals using advanced imaging technologies such as X-ray computed tomography (CT). Since carcass composition traits are heritable, they are potentially amenable to improvement through marker-assisted and genomic selection. We conducted a genome-wide association study (GWAS) on about 600 Scottish Blackface lambs for which detailed carcass composition phenotypes, including bone, fat and muscle components, had been captured using CT and which were genotyped for ~40,000 single nucleotide polymorphisms (SNPs) using the Illumina OvineSNP50 chip. RESULTS: We confirmed that the carcass composition traits were heritable with moderate to high (0.19-0.78) heritabilities. The GWAS analyses revealed multiple SNPs and quantitative trait loci (QTL) that were associated with effects on carcass composition traits and were significant at the genome-wide level. In particular, we identified a region on ovine chromosome 6 (OAR6) associated with bone weight and bone area that harboured SNPs with p values of 5.55 × 10(-8) and 2.63 × 10(-9), respectively. The same region had effects on fat area, fat density, fat weight and muscle density. We identified plausible positional candidate genes for these OAR6 QTL. We also detected a SNP that reached the genome-wide significance threshold with a p value of 7.28 × 10(-7) and was associated with muscle density on OAR1. Using a regional heritability mapping approach, we also detected regions on OAR3 and 24 that reached genome-wide significance for bone density. CONCLUSIONS: We identified QTL on OAR1, 3, 24 and particularly on OAR6 that are associated with effects on muscle, fat and bone traits. Based on available evidence that indicates that these traits are genetically correlated with meat quality traits, these associated SNPs have potential applications in selective breeding for improved meat quality. Further research is required to determine whether the effects associated with the OAR6 QTL are caused by a single gene or several closely-linked genes.

Genetic and genomic basis of antibody response to porcine reproductive and respiratory syndrome (PRRS) in gilts and sows.

BACKGROUND: Our recent research showed that antibody response to porcine reproductive and respiratory syndrome (PRRS), measured as sample-to-positive (S/P) ratio, is highly heritable and has a high genetic correlation with reproductive performance during a PRRS outbreak. Two major quantitative trait loci (QTL) on Sus scrofa chromosome 7 (SSC7; QTLMHC and QTL130) accounted for ~40 % of the genetic variance for S/P. Objectives of this study were to estimate genetic parameters for PRRS S/P in gilts during acclimation, identify regions associated with S/P, and evaluate the accuracy of genomic prediction of S/P across populations with different prevalences of PRRS and using different single nucleotide polymorphism (SNP) sets. METHODS: Phenotypes and high-density SNP genotypes of female pigs from two datasets were used. The outbreak dataset included 607 animals from one multiplier herd, whereas the gilt acclimation (GA) dataset included data on 2364 replacement gilts from seven breeding companies placed on health-challenged farms. Genomic prediction was evaluated using GA for training and validation, and using GA for training and outbreak for validation. Predictions were based on SNPs across the genome (SNPAll), SNPs in one (SNPMHC and SNP130) or both (SNPSSC7) QTL, or SNPs outside the QTL (SNPRest). RESULTS: Heritability of S/P in the GA dataset increased with the proportion of PRRS-positive animals in the herd (from 0.28 to 0.47). Genomic prediction accuracies ranged from low to moderate. Average accuracies were highest when using only the 269 SNPs in both QTL regions (SNPSSC7, with accuracies of 0.39 and 0.31 for outbreak and GA validation datasets, respectively. Average accuracies for SNPALL, SNPMHC, SNP130, and SNPRest were, respectively, 0.26, 0.39, 0.21, and 0.05 for the outbreak, and 0.28, 0.25, 0.22, and 0.12, for the GA validation datasets. CONCLUSIONS: Moderate genomic prediction accuracies can be obtained for PRRS antibody response using SNPs located within two major QTL on SSC7, while the rest of the genome showed limited predictive ability. Results were obtained using data from multiple genetic sources and farms, which further strengthens these findings. Further research is needed to validate the use of S/P ratio as an indicator trait for reproductive performance during PRRS outbreaks.

The genetic architecture of growth and fillet traits in farmed Atlantic salmon (Salmo salar).

BACKGROUND: Performance and quality traits such as harvest weight, fillet weight and flesh color are of economic importance to the Atlantic salmon aquaculture industry. The genetic factors underlying these traits are of scientific and commercial interest. However, such traits are typically polygenic in nature, with the number and size of QTL likely to vary between studies and populations. The aim of this study was to investigate the genetic basis of several growth and fillet traits measured at harvest in a large farmed salmon population by using SNP markers. Due to the marked heterochiasmy in salmonids, an efficient two-stage mapping approach was applied whereby QTL were detected using a sire-based linkage analysis, a sparse SNP marker map and exploiting low rates of recombination, while a subsequent dam-based analysis focused on the significant chromosomes with a denser map to confirm QTL and estimate their position. RESULTS: The harvest traits all showed significant heritability, ranging from 0.05 for fillet yield up to 0.53 for the weight traits. In the sire-based analysis, 1695 offspring with trait records and their 20 sires were successfully genotyped for the SNPs on the sparse map. Chromosomes 13, 18, 19 and 20 were shown to harbor genome-wide significant QTL affecting several growth-related traits. The QTL on chr. 13, 18 and 20 were detected in the dam-based analysis using 512 offspring from 10 dams and explained approximately 6-7 % of the within-family variation in these traits. CONCLUSIONS: We have detected several QTL affecting economically important complex traits in a commercial salmon population. Overall, the results suggest that the traits are relatively polygenic and that QTL tend to be pleiotropic (affecting the weight of several components of the harvested fish). Comparison of QTL regions across studies suggests that harvest trait QTL tend to be relatively population-specific. Therefore, the application of marker or genomic selection for improvement in these traits is likely to be most effective when the discovery population is closely related to the selection candidates (e.g. within-family genomic selection).

Linkage maps of the Atlantic salmon (Salmo salar) genome derived from RAD sequencing.

BACKGROUND: Genetic linkage maps are useful tools for mapping quantitative trait loci (QTL) influencing variation in traits of interest in a population. Genotyping-by-sequencing approaches such as Restriction-site Associated DNA sequencing (RAD-Seq) now enable the rapid discovery and genotyping of genome-wide SNP markers suitable for the development of dense SNP linkage maps, including in non-model organisms such as Atlantic salmon (Salmo salar). This paper describes the development and characterisation of a high density SNP linkage map based on SbfI RAD-Seq SNP markers from two Atlantic salmon reference families. RESULTS: Approximately 6,000 SNPs were assigned to 29 linkage groups, utilising markers from known genomic locations as anchors. Linkage maps were then constructed for the four mapping parents separately. Overall map lengths were comparable between male and female parents, but the distribution of the SNPs showed sex-specific patterns with a greater degree of clustering of sire-segregating SNPs to single chromosome regions. The maps were integrated with the Atlantic salmon draft reference genome contigs, allowing the unique assignment of ~4,000 contigs to a linkage group. 112 genome contigs mapped to two or more linkage groups, highlighting regions of putative homeology within the salmon genome. A comparative genomics analysis with the stickleback reference genome identified putative genes closely linked to approximately half of the ordered SNPs and demonstrated blocks of orthology between the Atlantic salmon and stickleback genomes. A subset of 47 RAD-Seq SNPs were successfully validated using a high-throughput genotyping assay, with a correspondence of 97% between the two assays. CONCLUSIONS: This Atlantic salmon RAD-Seq linkage map is a resource for salmonid genomics research as genotyping-by-sequencing becomes increasingly common. This is aided by the integration of the SbfI RAD-Seq SNPs with existing reference maps and the draft reference genome, as well as the identification of putative genes proximal to the SNPs. Differences in the distribution of recombination events between the sexes is evident, and regions of homeology have been identified which are reflective of the recent salmonid whole genome duplication.

Development and validation of a high density SNP genotyping array for Atlantic salmon (Salmo salar).

BACKGROUND: Dense single nucleotide polymorphism (SNP) genotyping arrays provide extensive information on polymorphic variation across the genome of species of interest. Such information can be used in studies of the genetic architecture of quantitative traits and to improve the accuracy of selection in breeding programs. In Atlantic salmon (Salmo salar), these goals are currently hampered by the lack of a high-density SNP genotyping platform. Therefore, the aim of the study was to develop and test a dense Atlantic salmon SNP array. RESULTS: SNP discovery was performed using extensive deep sequencing of Reduced Representation (RR-Seq), Restriction site-Associated DNA (RAD-Seq) and mRNA (RNA-Seq) libraries derived from farmed and wild Atlantic salmon samples (n = 283) resulting in the discovery of > 400 K putative SNPs. An Affymetrix Axiom® myDesign Custom Array was created and tested on samples of animals of wild and farmed origin (n = 96) revealing a total of 132,033 polymorphic SNPs with high call rate, good cluster separation on the array and stable Mendelian inheritance in our sample. At least 38% of these SNPs are from transcribed genomic regions and therefore more likely to include functional variants. Linkage analysis utilising the lack of male recombination in salmonids allowed the mapping of 40,214 SNPs distributed across all 29 pairs of chromosomes, highlighting the extensive genome-wide coverage of the SNPs. An identity-by-state clustering analysis revealed that the array can clearly distinguish between fish of different origins, within and between farmed and wild populations. Finally, Y-chromosome-specific probes included on the array provide an accurate molecular genetic test for sex. CONCLUSIONS: This manuscript describes the first high-density SNP genotyping array for Atlantic salmon. This array will be publicly available and is likely to be used as a platform for high-resolution genetics research into traits of evolutionary and economic importance in salmonids and in aquaculture breeding programs via genomic selection.

The utility of low-density genotyping for imputation in the Thoroughbred horse.

BACKGROUND: Despite the dramatic reduction in the cost of high-density genotyping that has occurred over the last decade, it remains one of the limiting factors for obtaining the large datasets required for genomic studies of disease in the horse. In this study, we investigated the potential for low-density genotyping and subsequent imputation to address this problem. RESULTS: Using the haplotype phasing and imputation program, BEAGLE, it is possible to impute genotypes from low- to high-density (50K) in the Thoroughbred horse with reasonable to high accuracy. Analysis of the sources of variation in imputation accuracy revealed dependence both on the minor allele frequency of the single nucleotide polymorphisms (SNPs) being imputed and on the underlying linkage disequilibrium structure. Whereas equidistant spacing of the SNPs on the low-density panel worked well, optimising SNP selection to increase their minor allele frequency was advantageous, even when the panel was subsequently used in a population of different geographical origin. Replacing base pair position with linkage disequilibrium map distance reduced the variation in imputation accuracy across SNPs. Whereas a 1K SNP panel was generally sufficient to ensure that more than 80% of genotypes were correctly imputed, other studies suggest that a 2K to 3K panel is more efficient to minimize the subsequent loss of accuracy in genomic prediction analyses. The relationship between accuracy and genotyping costs for the different low-density panels, suggests that a 2K SNP panel would represent good value for money. CONCLUSIONS: Low-density genotyping with a 2K SNP panel followed by imputation provides a compromise between cost and accuracy that could promote more widespread genotyping, and hence the use of genomic information in horses. In addition to offering a low cost alternative to high-density genotyping, imputation provides a means to combine datasets from different genotyping platforms, which is becoming necessary since researchers are starting to use the recently developed equine 70K SNP chip. However, more work is needed to evaluate the impact of between-breed differences on imputation accuracy.

A SNP in the 5' flanking region of the myostatin-1b gene is associated with harvest traits in Atlantic salmon (Salmo salar).

BACKGROUND: Myostatin (MSTN) belongs to the transforming growth factor-ß superfamily and is a potent negative regulator of skeletal muscle development and growth in mammals. Most teleost fish possess two MSTN paralogues. However, as a consequence of a recent whole genome-duplication event, salmonids have four: MSTN-1 (-1a and -1b) and MSTN-2 (-2a and -2b). Evidence suggests that teleost MSTN plays a role in the regulation of muscle growth. In the current study, the MSTN-1b gene was re-sequenced and screened for SNP markers in a commercial population of Atlantic salmon. After genotyping 4,800 progeny for the discovered SNPs, we investigated their association with eight harvest traits - four body-weight traits, two ratios of weight traits, flesh colour and fat percentage - using a mixed model association analysis. RESULTS: Three novel SNPs were discovered in the MSTN-1b gene of Atlantic salmon. One of the SNPs, located within the 5' flanking region (g.1086C > T), had a significant association with harvest traits (p < 0.05), specifically for: Harvest Weight (kg), Gutted Weight (kg), Deheaded Weight (kg) and Fillet Weight (kg). The haplotype-based association analysis was consistent with this result because the two haplotypes that showed a significant association with body-weight traits, hap4 and hap5 (p < 0.05 and p < 0.01, respectively), differ by a single substitution at the g.1086C > T locus. The alleles at g.1086C > T act in an additive manner and explain a small percentage of the genetic variation of these phenotypes. CONCLUSIONS: The association analysis revealed that g.1086C > T had a significant association with all body-weight traits under study. Although the SNP explains a small percentage of the variance, our results indicate that a variation in the 5' flanking region of the myostatin gene is associated with the genetic regulation of growth in Atlantic salmon.

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.

Modelling the short- and long-term impacts of drenching frequency and targeted selective treatment on the performance of grazing lambs and the emergence of anthelmintic resistance.

Refugia-based treatment strategies aim to prolong anthelmintic efficacy by maintaining a parasite population unexposed to anthelmintics. Targeted selective treatment (TST) achieves this by treating only animals that will benefit most from treatment, using a determinant criterion (DC). We developed a mathematical model to compare various traits proposed as DC, and investigate impacts of TST and drenching frequency on sheep performance and anthelmintic resistance. Short term, decreasing the proportion of animals drenched reduced benefits of anthelmintic treatment, assessed by empty body weight (EBW), but decreased the rate of anthelmintic resistance development; each consecutive drenching had a reduced impact on average EBW and an increased impact on the rate of anthelmintic resistance emergences. The optimal DC was fecal egg count, maintaining the highest average EBW when reducing the proportion of animals drenched. Long-term, reducing the proportion of animals drenched had little impact on total weight gain benefits, across animals and years, whilst reducing drenching frequency increased it. Decreasing the frequency and proportion of animals drenched were both predicted to increase the duration of anthelmintic efficacy but reduce the total number of drenches administered before resistance was observed. TST and frequency of drenching may lead to different benefits in the short versus long term.

Characterisation of QTL-linked and genome-wide restriction site-associated DNA (RAD) markers in farmed Atlantic salmon.

BACKGROUND: Restriction site-associated DNA sequencing (RAD-Seq) is a genome complexity reduction technique that facilitates large-scale marker discovery and genotyping by sequencing. Recent applications of RAD-Seq have included linkage and QTL mapping with a particular focus on non-model species. In the current study, we have applied RAD-Seq to two Atlantic salmon families from a commercial breeding program. The offspring from these families were classified into resistant or susceptible based on survival/mortality in an Infectious Pancreatic Necrosis (IPN) challenge experiment, and putative homozygous resistant or susceptible genotype at a major IPN-resistance QTL. From each family, the genomic DNA of the two heterozygous parents and seven offspring of each IPN phenotype and genotype was digested with the SbfI enzyme and sequenced in multiplexed pools. RESULTS: Sequence was obtained from approximately 70,000 RAD loci in both families and a filtered set of 6,712 segregating SNPs were identified. Analyses of genome-wide RAD marker segregation patterns in the two families suggested SNP discovery on all 29 Atlantic salmon chromosome pairs, and highlighted the dearth of male recombination. The use of pedigreed samples allowed us to distinguish segregating SNPs from putative paralogous sequence variants resulting from the relatively recent genome duplication of salmonid species. Of the segregating SNPs, 50 were linked to the QTL. A subset of these QTL-linked SNPs were converted to a high-throughput assay and genotyped across large commercial populations of IPNV-challenged salmon fry. Several SNPs showed highly significant linkage and association with resistance to IPN, and population linkage-disequilibrium-based SNP tests for resistance were identified. CONCLUSIONS: We used RAD-Seq to successfully identify and characterise high-density genetic markers in pedigreed aquaculture Atlantic salmon. These results underline the effectiveness of RAD-Seq as a tool for rapid and efficient generation of QTL-targeted and genome-wide marker data in a large complex genome, and its possible utility in farmed animal selection programs.

A genome-wide association study of osteochondritis dissecans in the Thoroughbred.

Osteochondrosis is a developmental orthopaedic disease that occurs in horses, other livestock species, companion animal species, and humans. The principal aim of this study was to identify quantitative trait loci (QTL) associated with osteochondritis dissecans (OCD) in the Thoroughbred using a genome-wide association study. A secondary objective was to test the effect of previously identified QTL in the current population. Over 300 horses, classified as cases or controls according to clinical findings, were genotyped for the Illumina Equine SNP50 BeadChip. An animal model was first implemented in order to adjust each horse's phenotypic status for average relatedness among horses and other potentially confounding factors which were present in the data. The genome-wide association test was then conducted on the residuals from the animal model. A single SNP on chromosome 3 was found to be associated with OCD at a genome-wide level of significance, as determined by permutation. According to the current sequence annotation, the SNP is located in an intergenic region of the genome. The effects of 24 SNPs, representing QTL previously identified in a sample of Hanoverian Warmblood horses, were tested directly in the animal model. When fitted alongside the significant SNP on ECA3, two of these SNPs were found to be associated with OCD. Confirmation of the putative QTL identified on ECA3 requires validation in an independent sample. The results of this study suggest that a significant challenge faced by equine researchers is the generation of sufficiently large data sets to effectively study complex diseases such as osteochondrosis.

Strengthening insights into host responses to mastitis infection in ruminants by combining heterogeneous microarray data sources.

BACKGROUND: Gene expression profiling studies of mastitis in ruminants have provided key but fragmented knowledge for the understanding of the disease. A systematic combination of different expression profiling studies via meta-analysis techniques has the potential to test the extensibility of conclusions based on single studies. Using the program Pointillist, we performed meta-analysis of transcription-profiling data from six independent studies of infections with mammary gland pathogens, including samples from cattle challenged in vivo with S. aureus, E. coli, and S. uberis, samples from goats challenged in vivo with S. aureus, as well as cattle macrophages and ovine dendritic cells infected in vitro with S. aureus. We combined different time points from those studies, testing different responses to mastitis infection: overall (common signature), early stage, late stage, and cattle-specific. RESULTS: Ingenuity Pathway Analysis of affected genes showed that the four meta-analysis combinations share biological functions and pathways (e.g. protein ubiquitination and polyamine regulation) which are intrinsic to the general disease response. In the overall response, pathways related to immune response and inflammation, as well as biological functions related to lipid metabolism were altered. This latter observation is consistent with the milk fat content depression commonly observed during mastitis infection. Complementarities between early and late stage responses were found, with a prominence of metabolic and stress signals in the early stage and of the immune response related to the lipid metabolism in the late stage; both mechanisms apparently modulated by few genes, including XBP1 and SREBF1.The cattle-specific response was characterized by alteration of the immune response and by modification of lipid metabolism. Comparison of E. coli and S. aureus infections in cattle in vivo revealed that affected genes showing opposite regulation had the same altered biological functions and provided evidence that E. coli caused a stronger host response. CONCLUSIONS: This meta-analysis approach reinforces previous findings but also reveals several novel themes, including the involvement of genes, biological functions, and pathways that were not identified in individual studies. As such, it provides an interesting proof of principle for future studies combining information from diverse heterogeneous sources.

Host inhibits replication of European porcine reproductive and respiratory syndrome virus in macrophages by altering differential regulation of type-I interferon transcriptional response.

Porcine reproductive and respiratory syndrome (PRRS) is an infectious disease caused by a positive RNA strand arterivirus. PRRS virus (PRRSV) interacts primarily with lung macrophages. Little is known how the virus subverts the innate immune response to initiate its replication in alveolar macrophages. Large-scale transcriptional responses of macrophages with different levels of susceptibility to PRRSV infection were compared over 30 h of infection. This study demonstrates a rapid and intense host transcriptional remodelling during the early phase of the replication of the virus which correlates with transient repression of type-I interferon transcript as early as 8 h post-infection. These results support the suggestion from previous studies that host innate immune response inhibits replication of European porcine reproductive and respiratory syndrome virus in macrophages by altering differential regulation of type-I interferon transcriptional response.

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.

SummaryGenetic resistance to gastrointestinal worms is a complex trait of great importance in both livestock and humans. In order to gain insights into the genetic architecture of this trait, a mixed breed population of sheep was artificially infected with Trichostrongylus colubriformis (n=3326) and then Haemonchus contortus (n=2669) to measure faecal worm egg count (WEC). The population was genotyped with the Illumina OvineSNP50 BeadChip and 48 640 single nucleotide polymorphism (SNP) markers passed the quality controls. An independent population of 316 sires of mixed breeds with accurate estimated breeding values for WEC were genotyped for the same SNP to assess the results obtained from the first population. We used principal components from the genomic relationship matrix among genotyped individuals to account for population stratification, and a novel approach to directly account for the sampling error associated with each SNP marker regression. The largest marker effects were estimated to explain an average of 0·48% (T. colubriformis) or 0·08% (H. contortus) of the phenotypic variance in WEC. These effects are small but consistent with results from other complex traits. We also demonstrated that methods which use all markers simultaneously can successfully predict genetic merit for resistance to worms, despite the small effects of individual markers. Correlations of genomic predictions with breeding values of the industry sires reached a maximum of 0·32. We estimate that effective across-breed predictions of genetic merit with multi-breed populations will require an average marker spacing of approximately 10 kbp.