A locus conferring tolerance to Theileria infection in African cattle.

East Coast fever, a tick-borne cattle disease caused by the Theileria parva parasite, is among the biggest natural killers of cattle in East Africa, leading to over 1 million deaths annually. Here we report on the genetic analysis of a cohort of Bos indicus (Boran) cattle demonstrating heritable tolerance to infection with T. parva (h2 = 0.65, s.e. 0.57). Through a linkage analysis we identify a 6 Mb genomic region on bovine chromosome 15 that is significantly associated with survival outcome following T. parva exposure. Testing this locus in an independent cohort of animals replicates this association with survival following T. parva infection. A stop gained variant in a paralogue of the FAF1 gene in this region was found to be highly associated with survival across both related and unrelated animals, with only one of the 20 homozygote carriers (T/T) of this change succumbing to the disease in contrast to 44 out of 97 animals homozygote for the reference allele (C/C). Consequently, we present a genetic locus linked to tolerance of one of Africa's most important cattle diseases, raising the promise of marker-assisted selection for cattle that are less susceptible to infection by T. parva.

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.

Harnessing longitudinal information to identify genetic variation in tolerance of pigs to Porcine Reproductive and Respiratory Syndrome virus infection.

BACKGROUND: High resistance (the ability of the host to reduce pathogen load) and tolerance (the ability to maintain high performance at a given pathogen load) are two desirable host traits for producing animals that are resilient to infections. For Porcine Reproductive and Respiratory Syndrome (PRRS), one of the most devastating swine diseases worldwide, studies have identified substantial genetic variation in resistance of pigs, but evidence for genetic variation in tolerance has so far been inconclusive. Resistance and tolerance are usually considered as static traits. In this study, we used longitudinal viremia measurements of PRRS virus infected pigs to define discrete stages of infection based on viremia profile characteristics. These were used to investigate host genetic effects on viral load (VL) and growth at different stages of infection, to quantify genetic variation in tolerance at these stages and throughout the entire 42-day observation period, and to assess whether the single nucleotide polymorphism (SNP) WUR10000125 (WUR) with known large effects on resistance confers significant differences in tolerance. RESULTS: Genetic correlations between resistance and growth changed considerably over time. Individuals that expressed high genetic resistance early in infection tended to grow slower during that time-period, but were more likely to experience lower VL and recovery in growth by the later stage. The WUR genotype was most strongly associated with VL at early- to mid-stages of infection, and with growth at mid- to late-stages of infection. Both, single-stage and repeated measurements random regression models identified significant genetic variation in tolerance. The WUR SNP was significantly associated only with the overall tolerance slope fitted through all stages of infection, with the genetically more resistant AB pigs for the WUR SNP being also more tolerant to PRRS. CONCLUSIONS: The results suggest that genetic selection for improved tolerance of pigs to PRRS is possible in principle, but may be feasible only with genomic selection, requiring intense recording schemes that involve repeated measurements to reliably estimate genetic effects. In the absence of such records, consideration of the WUR genotype in current selection schemes appears to be a promising strategy to improve simultaneously resistance and tolerance of growing pigs to PRRS.

Genomic regions underlying susceptibility to bovine tuberculosis in Holstein-Friesian cattle.

BACKGROUND: The significant social and economic loss as a result of bovine tuberculosis (bTB) presents a continuous challenge to cattle industries in the UK and worldwide. However, host genetic variation in cattle susceptibility to bTB provides an opportunity to select for resistant animals and further understand the genetic mechanisms underlying disease dynamics. METHODS: The present study identified genomic regions associated with susceptibility to bTB using genome-wide association (GWA), regional heritability mapping (RHM) and chromosome association approaches. Phenotypes comprised de-regressed estimated breeding values of 804 Holstein-Friesian sires and pertained to three bTB indicator traits: i) positive reactors to the skin test with positive post-mortem examination results (phenotype 1); ii) positive reactors to the skin test regardless of post-mortem examination results (phenotype 2) and iii) as in (ii) plus non-reactors and inconclusive reactors to the skin tests with positive post-mortem examination results (phenotype 3). Genotypes based on the 50 K SNP DNA array were available and a total of 34,874 SNPs remained per animal after quality control. RESULTS: The estimated polygenic heritability for susceptibility to bTB was 0.26, 0.37 and 0.34 for phenotypes 1, 2 and 3, respectively. GWA analysis identified a putative SNP on Bos taurus autosomes (BTA) 2 associated with phenotype 1, and another on BTA 23 associated with phenotype 2. Genomic regions encompassing these SNPs were found to harbour potentially relevant annotated genes. RHM confirmed the effect of these genomic regions and identified new regions on BTA 18 for phenotype 1 and BTA 3 for phenotypes 2 and 3. Heritabilities of the genomic regions ranged between 0.05 and 0.08 across the three phenotypes. Chromosome association analysis indicated a major role of BTA 23 on susceptibility to bTB. CONCLUSION: Genomic regions and candidate genes identified in the present study provide an opportunity to further understand pathways critical to cattle susceptibility to bTB and enhance genetic improvement programmes aiming at controlling and eradicating the disease.

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.

Genome-wide association studies of immune, disease and production traits in indigenous chicken ecotypes.

BACKGROUND: The majority of chickens in sub-Saharan Africa are indigenous ecotypes, well adapted to the local environment and raised in scavenging production systems. Although they are generally resilient to disease challenge, routine vaccination and biosecurity measures are rarely applied and infectious diseases remain a major cause of mortality and reduced productivity. Management and genetic improvement programmes are hampered by lack of routine data recording. Selective breeding based on genomic technologies may provide the means to enhance sustainability. In this study, we investigated the genetic architecture of antibody response to four major infectious diseases [infectious bursal disease (IBDV), Marek's disease (MDV), fowl typhoid (SG), fowl cholera (PM)] and resistance to Eimeria and cestode parasitism, along with two production traits [body weight and body condition score (BCS)] in two distinct indigenous Ethiopian chicken ecotypes. We conducted variance component analyses, genome-wide association studies, and pathway and selective sweep analyses. RESULTS: The large majority of birds was found to have antibody titres for all pathogens and were infected with both parasites, suggesting almost universal exposure. We derived significant moderate to high heritabilities for IBDV, MDV and PM antibody titres, cestodes infestation, body weight and BCS. We identified single nucleotide polymorphisms (SNPs) with genome-wide significance for each trait. Based on these associations, we identified for each trait, pathways, networks and functional gene clusters that include plausible candidate genes. Selective sweep analyses revealed a locus on chromosome 18 associated with viral antibody titres and resistance to Eimeria parasitism that is within a positive selection signal. We found no significant genetic correlations between production, immune and disease traits, implying that selection for altered antibody response and/or disease resistance will not affect production. CONCLUSIONS: We confirmed the presence of genetic variability and identified SNPs significantly associated with immune, disease and production traits in indigenous village chickens. Results underpin the feasibility of concomitant genetic improvement for enhanced antibody response, resistance to parasitism and productivity within and across indigenous chicken ecotypes.

Genomic prediction of host resistance to sea lice in farmed Atlantic salmon populations.

BACKGROUND: Sea lice have significant negative economic and welfare impacts on marine Atlantic salmon farming. Since host resistance to sea lice has a substantial genetic component, selective breeding can contribute to control of lice. Genomic selection uses genome-wide marker information to predict breeding values, and can achieve markedly higher accuracy than pedigree-based methods. Our aim was to assess the genetic architecture of host resistance to sea lice, and test the utility of genomic prediction of breeding values. Individual lice counts were measured in challenge experiments using two large Atlantic salmon post-smolt populations from a commercial breeding programme, which had genotypes for ~33 K single nucleotide polymorphisms (SNPs). The specific objectives were to: (i) estimate the heritability of host resistance; (ii) assess its genetic architecture by performing a genome-wide association study (GWAS); (iii) assess the accuracy of predicted breeding values using varying SNP densities (0.5 to 33 K) and compare it to that of pedigree-based prediction; and (iv) evaluate the accuracy of prediction in closely and distantly related animals. RESULTS: Heritability of host resistance was significant (0.22 to 0.33) in both populations using either pedigree or genomic relationship matrices. The GWAS suggested that lice resistance is a polygenic trait, and no genome-wide significant quantitative trait loci were identified. Based on cross-validation analysis, genomic predictions were more accurate than pedigree-based predictions for both populations. Although prediction accuracies were highest when closely-related animals were used in the training and validation sets, the benefit of having genomic-versus pedigree-based predictions within a population increased as the relationships between training and validation sets decreased. Prediction accuracy reached an asymptote with a SNP density of ~5 K within populations, although higher SNP density was advantageous for cross-population prediction. CONCLUSIONS: Host resistance to sea lice in farmed Atlantic salmon has a significant genetic component. Phenotypes relating to host resistance can be predicted with moderate to high accuracy within populations, with a major advantage of genomic over pedigree-based methods, even at relatively sparse SNP densities. Prediction accuracies across populations were low, but improved with higher marker densities. Genomic selection can contribute to lice control in salmon farming.

Genome wide association and genomic prediction for growth traits in juvenile farmed Atlantic salmon using a high density SNP array.

BACKGROUND: The genetic architecture of complex traits in farmed animal populations is of interest from a scientific and practical perspective. The use of genetic markers to predict the genetic merit (breeding values) of individuals is commonplace in modern farm animal breeding schemes. Recently, high density SNP arrays have become available for Atlantic salmon, which facilitates genomic prediction and association studies using genome-wide markers and economically important traits. The aims of this study were (i) to use a high density SNP array to investigate the genetic architecture of weight and length in juvenile Atlantic salmon; (ii) to assess the utility of genomic prediction for these traits, including testing different marker densities; (iii) to identify potential candidate genes underpinning variation in early growth. RESULTS: A pedigreed population of farmed Atlantic salmon (n = 622) were measured for weight and length traits at one year of age, and genotyped for 111,908 segregating SNP markers using a high density SNP array. The heritability of both traits was estimated using pedigree and genomic relationship matrices, and was comparable at around 0.5 and 0.6 respectively. The results of the GWA analysis pointed to a polygenic genetic architecture, with no SNPs surpassing the genome-wide significance threshold, and one SNP associated with length at the chromosome-wide level. SNPs surpassing an arbitrary threshold of significance (P < 0.005, ~ top 0.5 % of markers) were aligned to an Atlantic salmon reference transcriptome, identifying 109 SNPs in transcribed regions that were annotated by alignment to human, mouse and zebrafish protein databases. Prediction of breeding values was more accurate when applying genomic (GBLUP) than pedigree (PBLUP) relationship matrices (accuracy ~ 0.7 and 0.58 respectively) and 5,000 SNPs were sufficient for obtaining this accuracy increase over PBLUP in this specific population. CONCLUSIONS: The high density SNP array can effectively capture the additive genetic variation in complex traits. However, the traits of weight and length both appear to be very polygenic with only one SNP surpassing the chromosome-wide threshold. Genomic prediction using the array is effective, leading to an improvement in accuracy compared to pedigree methods, and this improvement can be achieved with only a small subset of the markers in this population. The results have practical relevance for genomic selection in salmon and may also provide insight into variation in the identified genes underpinning body growth and development in salmonid species.

The potential for use of haematological and anti-IgE humoral responses as phenotypic markers for tick resistance in cattle.

Approximately 80% of the global cattle population is at risk of infestation and infection by ticks and tick-borne diseases (TTBDs). The economic losses from animal mortality, reduced production, vector control costs and animal treatment are very substantial, hence there is an urgent need to develop and deploy alternative vector control strategies. Breeding for host tick resistance has the potential for sustainable large-scale TTBD control especially in cattle. The gold standard method for phenotyping tick resistance in cattle is by counting ticks on the body but is very laborious and subjective. Better methods for phenotyping tick resistance more objectively, faster and at scale, are essential for selecting host genetic resistance to ticks. This study investigated the correlation between haematological cellular profiles and immunological responses (immunoglobulin E, IgE) and full body tick counts in herds of Bos indicus and Bos taurus following artificial tick challenge with Rhipicephalus decoloratus larvae. Fifty-four Friesian and Ayrshire (Bos taurus) and 52 East African Zebu (Bos indicus) calves were each infested with ∼2500 larvae. Near-replete adult female ticks (≥ 4.5 mm) were counted daily from Day 20-25. Blood and serum samples were obtained from each animal on Days 0 and 23 for cellular blood and IgE titre analysis, respectively. The indicine cattle were refractory to R. decoloratus infestation in comparison with the taurine breed (P < 0.0001). Repeated measurements of blood components pre-infestation revealed a significant (P < 0.05) association with tick count in IgE and red blood cells, haematocrit, and haemoglobin post-infestation. There was also a strong positive correlation between the tick counts and red blood cell numbers, haemoglobin, haematocrit, and IgE concentration (P < 0.0001) following tick challenge. The application of this approach to phenotype host resistance needs to be assessed using higher cattle numbers and with different tick species or genera.

Investigative power of Genomic Informational Field Theory (GIFT) relative to GWAS for genotype-phenotype mapping.

Identifying associations between phenotype and genotype is the fundamental basis of genetic analyses. Inspired by frequentist probability and the work of R.A. Fisher, genome-wide association studies (GWAS) extract information using averages and variances from genotype-phenotype datasets. Averages and variances are legitimated upon creating distribution density functions obtained through the grouping of data into categories. However, as data from within a given category cannot be differentiated, the investigative power of such methodologies is limited. Genomic Informational Field Theory (GIFT) is a method specifically designed to circumvent this issue. The way GIFT proceeds is opposite to that of GWAS. Whilst GWAS determines the extent to which genes are involved in phenotype formation (bottom-up approach), GIFT determines the degree to which the phenotype can select microstates (genes) for its subsistence (top-down approach). Doing so requires dealing with new genetic concepts, a.k.a. genetic paths, upon which significance levels for genotype-phenotype associations can be determined. By using different datasets obtained in ovis aries related to bone growth (Dataset-1) and to a series of linked metabolic and epigenetic pathways (Dataset-2), we demonstrate that removing the informational barrier linked to categories enhances the investigative and discriminative powers of GIFT, namely that GIFT extracts more information than GWAS. We conclude by suggesting that GIFT is an adequate tool to study how phenotypic plasticity and genetic assimilation are linked.

Exploring the value of genomic predictions to simultaneously improve production potential and resilience of farmed animals.

Sustainable livestock production requires that animals have a high production potential but are also highly resilient to environmental challenges. The first step to simultaneously improve these traits through genetic selection is to accurately predict their genetic merit. In this paper, we used simulations of sheep populations to assess the effect of genomic data, different genetic evaluation models and phenotyping strategies on prediction accuracies and bias for production potential and resilience. In addition, we also assessed the effect of different selection strategies on the improvement of these traits. Results show that estimation of both traits greatly benefits from taking repeated measurements and from using genomic information. However, the prediction accuracy for production potential is compromised, and resilience estimates tends to be upwards biased, when families are clustered in groups even when genomic information is used. The prediction accuracy was also found to be lower for both traits, resilience and production potential, when the environment challenge levels are unknown. Nevertheless, we observe that genetic gain in both traits can be achieved even in the case of unknown environmental challenge, when families are distributed across a large range of environments. Simultaneous genetic improvement in both traits however greatly benefits from the use of genomic evaluation, reaction norm models and phenotyping in a wide range of environments. Using models without the reaction norm in scenarios where there is a trade-off between resilience and production potential, and phenotypes are collected from a narrow range of environments may result in a loss for one trait. The study demonstrates that genomic selection coupled with reaction-norm models offers great opportunities to simultaneously improve productivity and resilience of farmed animals even in the case of a trade-off.

Investigating volatile semiochemical production from Bos taurus and Bos indicus as a novel phenotype for breeding host resistance to ixodid ticks.

Ticks and tick-borne diseases cause significant loss in livestock production with about 80% world's cattle at risk. The cost of chemical control is high and there is an ever-increasing tick resistance to chemical acaricides. Genetic selection as alternative long-term control strategy is constrained by laborious phenotyping using tick counts or scores. This study explored the use of host volatile semiochemicals that may be attractants or repellents to ticks as a phenotype for new tick resistance, with potential to be used as a proxy in selection programmes. Approximately 100 young cattle composed of Bos indicus and Bos taurus were artificially infested with 2,500 African blue tick, Rhipicephalus decoloratus larvae, with daily female tick (4.5 mm) counts taken from day 20 post-infestation. Volatile organic compounds were sampled from cattle before and after tick infestation by dynamic headspace collection, analysed by high-resolution gas chromatography (GC) and subjected to multivariate statistical analysis. Using 6-day repeated measure analysis, three pre-infestation GC peaks (BI938 - unknown, BI966 - 6-methyl-5-hepten-2-one and BI995 - hexyl acetate) and one post-infestation GC peak (AI933 - benzaldehyde / (E)-2-heptenal) were associated with tick resistance (P < 0.01 and P < 0.05 respectively). The high correlation coefficients (r = 0.66) between repeated records with all volatile compounds support the potential predictive value for volatile compounds in selective breeding programmes for tick resistance in cattle.

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.

Multiple Country and Breed Genomic Prediction of Tick Resistance in Beef Cattle.

Ticks cause substantial production losses for beef and dairy cattle. Cattle resistance to ticks is one of the most important factors affecting tick control, but largely neglected due to the challenge of phenotyping. In this study, we evaluate the pooling of tick resistance phenotyped reference populations from multi-country beef cattle breeds to assess the possibility of improving host resistance through multi-trait genomic selection. Data consisted of tick counts or scores assessing the number of female ticks at least 4.5 mm length and derived from seven populations, with breed, country, number of records and genotyped/phenotyped animals being respectively: Angus (AN), Brazil, 2,263, 921/1,156, Hereford (HH), Brazil, 6,615, 1,910/2,802, Brangus (BN), Brazil, 2,441, 851/851, Braford (BO), Brazil, 9,523, 3,062/4,095, Tropical Composite (TC), Australia, 229, 229/229, Brahman (BR), Australia, 675, 675/675, and Nguni (NG), South Africa, 490, 490/490. All populations were genotyped using medium density Illumina SNP BeadChips and imputed to a common high-density panel of 332,468 markers. The mean linkage disequilibrium (LD) between adjacent SNPs varied from 0.24 to 0.37 across populations and so was sufficient to allow genomic breeding values (GEBV) prediction. Correlations of LD phase between breeds were higher between composites and their founder breeds (0.81 to 0.95) and lower between NG and the other breeds (0.27 and 0.35). There was wide range of estimated heritability (0.05 and 0.42) and genetic correlation (-0.01 and 0.87) for tick resistance across the studied populations, with the largest genetic correlation observed between BN and BO. Predictive ability was improved under the old-young validation for three of the seven populations using a multi-trait approach compared to a single trait within-population prediction, while whole and partial data GEBV correlations increased in all cases, with relative improvements ranging from 3% for BO to 64% for TC. Moreover, the multi-trait analysis was useful to correct typical over-dispersion of the GEBV. Results from this study indicate that a joint genomic evaluation of AN, HH, BN, BO and BR can be readily implemented to improve tick resistance of these populations using selection on GEBV. For NG and TC additional phenotyping will be required to obtain accurate GEBV.

Integrating Genetic and Genomic Analyses of Combined Health Data Across Ecotypes to Improve Disease Resistance in Indigenous African Chickens.

Poultry play an important role in the agriculture of many African countries. The majority of chickens in sub-Saharan Africa are indigenous, raised in villages under semi-scavenging conditions. Vaccinations and biosecurity measures rarely apply, and infectious diseases remain a major cause of mortality and reduced productivity. Genomic selection for disease resistance offers a potentially sustainable solution but this requires sufficient numbers of individual birds with genomic and phenotypic data, which is often a challenge to collect in the small populations of indigenous chicken ecotypes. The use of information across-ecotypes presents an attractive possibility to increase the relevant numbers and the accuracy of genomic selection. In this study, we performed a joint analysis of two distinct Ethiopian indigenous chicken ecotypes to investigate the genomic architecture of important health and productivity traits and explore the feasibility of conducting genomic selection across-ecotype. Phenotypic traits considered were antibody response to Infectious Bursal Disease (IBDV), Marek's Disease (MDV), Fowl Cholera (PM) and Fowl Typhoid (SG), resistance to Eimeria and cestode parasitism, and productivity [body weight and body condition score (BCS)]. Combined data from the two chicken ecotypes, Horro (n = 384) and Jarso (n = 376), were jointly analyzed for genetic parameter estimation, genome-wide association studies (GWAS), genomic breeding value (GEBVs) calculation, genomic predictions, whole-genome sequencing (WGS), and pathways analyses. Estimates of across-ecotype heritability were significant and moderate in magnitude (0.22-0.47) for all traits except for SG and BCS. GWAS identified several significant genomic associations with health and productivity traits. The WGS analysis revealed putative candidate genes and mutations for IBDV (TOLLIP, ANGPTL5, BCL9, THEMIS2), MDV (GRM7), SG (MAP3K21), Eimeria (TOM1L1) and cestodes (TNFAIP1, ATG9A, NOS2) parasitism, which warrant further investigation. Reliability of GEBVs increased compared to within-ecotype calculations but accuracy of genomic prediction did not, probably because the genetic distance between the two ecotypes offset the benefit from increased sample size. However, for some traits genomic prediction was only feasible in across-ecotype analysis. Our results generally underpin the potential of genomic selection to enhance health and productivity across-ecotypes. Future studies should establish the required minimum sample size and genetic similarity between ecotypes to ensure accurate joint genomic selection.

Traits associated with innate and adaptive immunity in pigs: heritability and associations with performance under different health status conditions.

There is a need for genetic markers or biomarkers that can predict resistance towards a wide range of infectious diseases, especially within a health environment typical of commercial farms. Such markers also need to be heritable under these conditions and ideally correlate with commercial performance traits. In this study, we estimated the heritabilities of a wide range of immune traits, as potential biomarkers, and measured their relationship with performance within both specific pathogen-free (SPF) and non-SPF environments. Immune traits were measured in 674 SPF pigs and 606 non-SPF pigs, which were subsets of the populations for which we had performance measurements (average daily gain), viz. 1549 SPF pigs and 1093 non-SPF pigs. Immune traits measured included total and differential white blood cell counts, peripheral blood mononuclear leucocyte (PBML) subsets (CD4+ cells, total CD8alpha+ cells, classical CD8alphabeta+ cells, CD11R1+ cells (CD8alpha+ and CD8alpha-), B cells, monocytes and CD16+ cells) and acute phase proteins (alpha-1 acid glycoprotein (AGP), haptoglobin, C-reactive protein (CRP) and transthyretin). Nearly all traits tested were heritable regardless of health status, although the heritability estimate for average daily gain was lower under non-SPF conditions. There were also negative genetic correlations between performance and the following immune traits: CD11R1+ cells, monocytes and the acute phase protein AGP. The strength of the association between performance and AGP was not affected by health status. However, negative genetic correlations were only apparent between performance and monocytes under SPF conditions and between performance and CD11R1+ cells under non-SPF conditions. Although we cannot infer causality in these relationships, these results suggest a role for using some immune traits, particularly CD11R1+ cells or AGP concentrations, as predictors of pig performance under the lower health status conditions associated with commercial farms.

Arginine to Glutamine Variant in Olfactomedin Like 3 (OLFML3) Is a Candidate for Severe Goniodysgenesis and Glaucoma in the Border Collie Dog Breed.

Goniodysgenesis is a developmental abnormality of the anterior chamber of the eye. It is generally considered to be congenital in dogs (Canis lupus familiaris), and has been associated with glaucoma and blindness. Goniodysgenesis and early-onset glaucoma initially emerged in Border Collies in Australia in the late 1990s and have subsequently been found in this breed in Europe and the USA. The objective of the present study was to determine the genetic basis of goniodysgenesis in Border Collies. Clinical diagnosis was based on results of examinations by veterinary ophthalmologists of affected and unaffected dogs from eleven different countries. Genotyping using the Illumina high density canine single nucleotide variant genotyping chip was used to identify a candidate genetic region. There was a highly significant peak of association over chromosome 17, with a p-value of 2 × 10-13 Expression profiles and evolutionary conservation of candidate genes were assessed using public databases. Whole genome sequences of three dogs with glaucoma, three severely affected by goniodysgenesis and three unaffected dogs identified a missense variant in the olfactomedin like 3 (OLFML3) gene in all six affected animals. This was homozygous for the risk allele in all nine cases with glaucoma and 12 of 14 other severely affected animals. Of 67 reportedly unaffected animals, only one was homozygous for this variant (offspring of parents both with goniodysgenesis who were also homozygous for the variant). Analysis of pedigree information was consistent with an autosomal recessive mode of inheritance for severe goniodysgenesis (potentially leading to glaucoma) in this breed. The identification of a candidate genetic region and putative causative variant will aid breeders to reduce the frequency of goniodysgenesis and the risk of glaucoma in the Border Collie population.

Genome-Wide Association and Genomic Selection for Resistance to Amoebic Gill Disease in Atlantic Salmon.

Amoebic gill disease (AGD) is one of the largest threats to salmon aquaculture, causing serious economic and animal welfare burden. Treatments can be expensive and environmentally damaging, hence the need for alternative strategies. Breeding for disease resistance can contribute to prevention and control of AGD, providing long-term cumulative benefits in selected stocks. The use of genomic selection can expedite selection for disease resistance due to improved accuracy compared to pedigree-based approaches. The aim of this work was to quantify and characterize genetic variation in AGD resistance in salmon, the genetic architecture of the trait, and the potential of genomic selection to contribute to disease control. An AGD challenge was performed in ∼1,500 Atlantic salmon, using gill damage and amoebic load as indicator traits for host resistance. Both traits are heritable (h2 ∼0.25-0.30) and show high positive correlation, indicating they may be good measurements of host resistance to AGD. While the genetic architecture of resistance appeared to be largely polygenic in nature, two regions on chromosome 18 showed suggestive association with both AGD resistance traits. Using a cross-validation approach, genomic prediction accuracy was up to 18% higher than that obtained using pedigree, and a reduction in marker density to ∼2,000 SNPs was sufficient to obtain accuracies similar to those obtained using the whole dataset. This study indicates that resistance to AGD is a suitable trait for genomic selection, and the addition of this trait to Atlantic salmon breeding programs can lead to more resistant stocks.

Genomic Selection for Growth Traits in Pacific Oyster (Crassostrea gigas): Potential of Low-Density Marker Panels for Breeding Value Prediction.

Pacific oysters are a key aquaculture species globally, and genetic improvement via selective breeding is a major target. Genomic selection has the potential to expedite genetic gain for key target traits of a breeding program, but has not yet been evaluated in oyster. The recent development of SNP arrays for Pacific oyster (Crassostrea gigas) raises the opportunity to test genomic selection strategies for polygenic traits. In this study, a population of 820 oysters (comprising 23 full-sibling families) were genotyped using a medium density SNP array (23 K informative SNPs), and the genetic architecture of growth-related traits [shell height (SH), shell length (SL), and wet weight (WW)] was evaluated. Heritability was estimated to be moderate for the three traits (0.26 ± 0.06 for SH, 0.23 ± 0.06 for SL and 0.35 ± 0.05 for WW), and results of a GWAS indicated that the underlying genetic architecture was polygenic. Genomic prediction approaches were used to estimate breeding values for growth, and compared to pedigree based approaches. The accuracy of the genomic prediction models (GBLUP) outperformed the traditional pedigree approach (PBLUP) by ∼25% for SL and WW, and ∼30% for SH. Further, reduction in SNP marker density had little impact on prediction accuracy, even when density was reduced to a few hundred SNPs. These results suggest that the use of genomic selection in oyster breeding could offer benefits for the selection of breeding candidates to improve complex economic traits at relatively modest cost.

The Genomic Architecture of Fowl Typhoid Resistance in Commercial Layers.

Salmonella enterica serovar Gallinarum causes devastating outbreaks of fowl typhoid across the globe, especially in developing countries. With the use of antimicrobial agents being reduced due to legislation and the absence of licensed vaccines in some parts of the world, an attractive complementary control strategy is to breed chickens for increased resistance to Salmonella. The potential for genetic control of salmonellosis has been demonstrated by experimental challenge of inbred populations. Quantitative trait loci (QTL) associated with resistance have been identified in many genomic regions. A major QTL associated with systemic salmonellosis has been identified in a region termed SAL1. In the present study, two outbreaks of fowl typhoid in 2007 and 2012 in the United Kingdom were used to investigate the genetic architecture of Salmonella resistance in commercial laying hens. In the first outbreak 100 resistant and 150 susceptible layers were genotyped using 11 single nucleotide polymorphism (SNP) and 3 microsatellite markers located in the previously identified SAL1 region on chromosome 5. From the second outbreak 100 resistant and 200 susceptible layers, belonging to a different line, were genotyped with a high-density (600 K) genome-wide SNP array. Substantial heritability estimates were obtained in both populations (h 2 = 0.22 and 0.26, for the layers in the first and second outbreak, respectively). Significant associations with three markers on chromosome 5 located close to AKT1 and SIVA1 genes, coding for RAC-alpha serine/threonine protein kinase, and the CD27-binding protein SIVA1, respectively, were identified in the first outbreak. From analysis of the second outbreak, eight genome-wide significant associations with Salmonella resistance were identified on chromosomes 1, 6, 7, 11, 23, 24, 26, 28 and several others with suggestive genome-wide significance were found. Pathway and network analysis revealed the presence of many innate immune pathways related to Salmonella resistance. Although, significant associations with SNPs located in the SAL1 locus were not identified by the genome-wide scan for layers from the second outbreak, pathway analysis revealed P13K/AKT signaling as the most significant pathway. In summary, resistance to fowl typhoid is a heritable polygenic trait that could possibly be enhanced through selective breeding.