Early life jumping traits: Are they good proxies for success in show jumping competitions in Belgian warmblood horses?

The main goal of the Belgian Warmblood horse studbook (BWP) is to breed successful competition horses, with emphasis on show jumping. However, competition results are only available later in life and competition traits are lowly heritable. Hence, the use of phenotypes that record performance-related traits at an early life stage could help increase genetic progress. In this study, we evaluated the potential of eleven linear scored early life jumping traits assessed during jumping in freedom (2-5 years old) or under the saddle (4-6 years old) as proxies for later success in show jumping competitions. To this end, we estimated their heritabilities and genetic correlations with the competition trait, adjusted fence height, by using 2170 free jumping records, 1588 jumping under saddle records, 674,527 show jumping competition records and almost 81,000 informative horses in the pedigree. As participation of young horses in these contests is on a voluntary basis, a pre-selection most probably exists. To verify this hypothesis, we investigated the association between participation to young horse contests and participation to show jumping competitions later on (called here start status phenotype). We also estimated heritabilities for "start status in free jumping contest", "start status in jumping under saddle contest" and "start status in free jumping or jumping under saddle contest" by fitting threshold models. Furthermore, we calculated genetic correlations between these traits and adjusted fence height and calculated the correlations between EBVs for start status in young horse contests and EBVs for success in competitions. Estimated heritabilities of early life jumping traits ranged between 0.05 and 0.30. Their genetic correlations with adjusted fence height were moderate to high (rg  = 0.37-0.63). Relatively more horses that participated in young horse contests competed later on compared to horses that did not participate in young horse contests (p-value < 0.001). They were also significantly more successful in show jumping competitions. Furthermore, start status in young horse contests was moderately heritable in BWP horses (h2  = 0.56-0.65) and moderately to highly correlated with later success in competitions (rg  = 0.30-0.77). Hence, we showed that ELJ traits are good proxies for later success in competitions and that a pre-selection of horses occurs in young horse contests. It is suggested to stimulate participation to young horse contests to achieve a more representative sample of the population. Early life jumping traits can therefore optimize the genetic progress for show jumping performance.

Effect of sire type and a by-product based diet on performance and meat quality in growing-finishing pigs

For many years, pig production has focused on maximizing performance by selecting for maximal muscle growth and feeding diets that allow the animals to express their genetic potential. However, it is unclear whether this selection for muscle deposition has affected the capacity of pigs to cope with by-product-based diets, which rely on fat as the primary energy source instead of starches and sugars. Therefore, an experiment was set up to investigate if different types of boars affect how their progeny cope with alternative ingredients in the diet, with a possible need for adapted breeding schemes. Two types of boars within the Piétrain sire line were used based on either a high or low estimated breeding value for daily feed intake (HFI: high feed intake, low feed intake). When their progeny reached 14 weeks of age, two dietary strategies were compared: a control (CON) vs a by-product-based diet high in fat and fiber (HFF). The CON diet was mainly based on cereals (corn, wheat, barley) and soybean meal. The HFF diet was formulated to contain the same net energy, CP and digestible amino acid levels without any cereals or soybean meal. In total 192 animals were included in the experiment (48 animals/type of boar/diet) and performance, digestibility, carcass and meat quality were compared. None of the parameters showed a significant interaction (P < 0.05) between the type of boar and diet, suggesting that shifting to diets that are less prone to feed-food competition is equally feasible in different types of pigs. Type of boar did affect performance, carcass quality and intramuscular fat content. HFI pigs showed higher daily feed intake (DFI) and daily gain (P < 0.001), with no significant difference in feed conversion ratio (P = 0.205), lower carcass quality (P < 0.001) and higher intramuscular fat content (P = 0.030). For both boar types, pigs fed the CON diet performed better, with a higher daily gain (P = 0.028), DFI (P = 0.011) and dressing yield (P = 0.009) and better digestibility (P < 0.001), but without differences in feed conversion ratio or meat quality. In conclusion, there was no indication that pigs differing in feed intake capacity cope differently with a high-fat, high-fiber diet based on by-products. Different types of pigs may cope well with diets that are less prone to feed-food competition.

The novel IMAGE001 genotyping array as a valuable alternative for genetic diversity screening in chicken: a demonstration in a local chicken breed in Belgium

Screening for genetic diversity in livestock species breeds is of utmost importance, especially for local, small populations that are at the risk of extinction. Luckily, recent developments in technology increase access to genotyping, also for numerically small breeds. One of these new technologies is the IMAGE001 single nucleotide polymorphism genotyping array that includes markers for 6 different species (cow, pig, sheep, chicken, horse and goat). For our current study, we studied the Turkey-headed Malines chicken, a local chicken breed in Belgium, for the first time. A total of 110 animals were genotyped, together with 29 samples from 4 supposedly related breeds. The genotypes were used to assess the genetic diversity of this local breed. Our analysis revealed an average inbreeding coefficient of 0.20 through runs of homozygosity analysis, and effective population size estimation based on linkage disequilibrium indicated a low genetic diversity (Ne = 34). Moreover, a principal component analysis and a genetic differentiation study (FST) were performed using these marker data to position the Turkey-headed Malines relative to the 4 other indigenous Belgian chicken breeds. Finally, we discussed the practical implications of the overlap between the IMAGE001 array and other existing chicken genotyping arrays. This study is the first use of the novel IMAGE001 array to evaluate a local chicken breed, and demonstrates it as a viable option for genomic characterization a breed. Moreover, with this research, we are able to provide a good basis for further evaluation of the Belgian chicken heritage.

The forgotten variable? Does the euthanasia method and sample storage condition influence an organisms transcriptome - a gene expression analysis on multiple tissues in pigs

BACKGROUND

Transcriptomic studies often require collection of fresh tissues post euthanasia. The chosen euthanasia method might have the potential to induce variations in gene expressions that are unlinked with the experimental design. The present study compared the suitability of 'nitrogen gas in foam' (ANOXIA) in comparison to a non-barbiturate anaesthetic, T-61® (T61), for euthanizing piglets used in transcriptome research. Further, the effect of common tissue storage conditions, RNAlater™ (RL) and snap freezing in liquid nitrogen (LN2), on gene expression profiles were also analysed.

RESULTS

On comparison of the 3'mRNA-Seq data generated from pituitary, hypothalamus, liver and lung tissues, no significant differential expression in the protein coding genes were detected between the euthanasia methods. This implies that the nitrogen anoxia method could be a suitable alternative for euthanasia of piglets used in transcriptomic research. However, small nuclear RNAs (snRNAs) that constitute the eukaryotic spliceosomal machinery were found to be significantly higher (log2fold change ≥ 2.0, and adjusted p value ≤ 0.1) in pituitary samples collected using ANOXIA. Non-protein coding genes like snRNAs that play an important role in pre-mRNA splicing can subsequently modify gene expression. Storage in RL was found to be superior in preserving RNA compared to LN2 storage, as evidenced by the significantly higher RIN values in representative samples. However, storage in RL as opposed to LN2, also influenced differential gene expression in multiple tissues, perhaps as a result of its inability to inhibit biological activity during storage. Hence such external sources of variations should be carefully considered before arriving at research conclusions.

CONCLUSIONS

Source of biological variations like euthanasia method and storage condition can confound research findings. Even if we are unable to prevent the effect of these external factors, it will be useful to identify the impact of these variables on the parameter under observation and thereby prevent misinterpretation of our results.

A promising resilience parameter for breeding: the use of weight and feed trajectories in growing pigs

BACKGROUND

Increasing resilience is a priority in modern pig breeding. Recent research shows that general resilience can be quantified via variability in longitudinal data. The collection of such longitudinal data on weight, feed intake and feeding behaviour in pigs has been facilitated by the development of technologies such as automated feeding stations. The goal of this study was to investigate resilience traits, which were estimated as deviations from longitudinal weight, feed intake and feeding behaviour data during the finishing phase. A dataset with 324,207 records between the age of 95 and 155 days on 5,939 Piétrain pigs with known pedigree and genomic information was used. We provided guidelines for a rigid quality control of longitudinal body weight data, as we found that outliers can significantly affect results. Gompertz growth curve analysis, linear modelling and trajectory analyses were used for quantifying resilience traits.

RESULTS

To our knowledge, this is the first study comparing resilience traits from longitudinal body weight, feed intake and feeding behaviour data in pigs. We demonstrated that the resilience traits are lowly to moderately heritable for deviations in body weight (h² = 2.9%-20.2%), in feed intake (9.4%-23.3%) and in feeding behaviour (16.2%-28.3%). Additionally, these traits have good predictive abilities in cross-validation analyses. Deviations in individual body weight and feed intake trajectories are highly correlated (r_g = 0.78) with low to moderate favourable genetic correlations with feed conversion ratio (r_g = 0.39-0.49). Lastly, we showed that some resilience traits, such as the natural logarithm of variances of observed versus predicted body weights (lnvar_weight), are more robust to lower observation frequencies and are repeatable over three different time periods of the finishing phase.

CONCLUSIONS

Our results will help future studies investigating resilience traits and resilience-related traits. Moreover, our study provides first results on standardization of quality control and efficient data sampling from automated feeding station data. Our findings will be valuable for breeding organizations as they offer evidence that pigs' general resilience can be selected on with good accuracy. Moreover, this methodology might be extended to other species to quantify resilience based on longitudinal data.

Adjusted fence height: an improved phenotype for the genetic evaluation of show jumping performance in Warmblood horses

BACKGROUND

Show jumping is one of the most popular disciplines in the horse sector, which makes success in show jumping competitions an important breeding goal for many studbooks. Therefore, the genetic evaluation of show jumping performance is of major interest and this is the case for two Belgian Warmblood studbooks: the Belgian Warmblood horse and Zangersheide. In this study, first an improved phenotype for show jumping performance was developed, i.e. adjusted fence height based on a new non-arbitrary method to scale ranking and competition level, which are two major components of success in competitions. Second, we assessed the importance of including a rider effect in genetic models for show jumping performance, this effect being under debate in sport horse breeding. Third, genetic models based on elementary performances and one model based on a summarized performance were compared in terms of model fit, heritabilities and the stability of estimated breeding values to define the most suitable one for the genetic evaluation of show jumping performance.

RESULTS

In this study, more than 600,000 Belgian competition records and almost 81,000 horses were used. Genetic evaluations were developed based on elementary performances (Blom-transformed ranking and adjusted fence height) and on a summarized performance (highest level achieved). Estimated heritabilities of Blom-transformed ranking, adjusted fence height and highest level achieved were 0.09, 0.12 and 0.39, respectively. Including a rider effect improved the models for genetic evaluations. Estimated genetic correlations between the studied models were moderate to high (r_g = 0.60-0.99). With the best fit model, the accuracy of the estimated breeding value (EBV) for adjusted fence height reached 0.70 for a larger number of stallions and for stallions that tended to be younger.

CONCLUSIONS

We recommend breeders to implement this new phenotype 'adjusted fence height' in breeding programs. It is moderately to highly correlated with Blom-transformed ranking and highest level achieved, a proxy for lifetime success, and is available for selection candidates from an early age onwards.

Estimating genetics of body dimensions and activity levels in pigs using automated pose estimation

Pig breeding is changing rapidly due to technological progress and socio-ecological factors. New precision livestock farming technologies such as computer vision systems are crucial for automated phenotyping on a large scale for novel traits, as pigs’ robustness and behavior are gaining importance in breeding goals. However, individual identification, data processing and the availability of adequate (open source) software currently pose the main hurdles. The overall goal of this study was to expand pig weighing with automated measurements of body dimensions and activity levels using an automated video-analytic system: DeepLabCut. Furthermore, these data were coupled with pedigree information to estimate genetic parameters for breeding programs. We analyzed 7428 recordings over the fattening period of 1556 finishing pigs (Piétrain sire x crossbred dam) with two-week intervals between recordings on the same pig. We were able to accurately estimate relevant body parts with an average tracking error of 3.3 cm. Body metrics extracted from video images were highly heritable (61–74%) and significantly genetically correlated with average daily gain (r_g = 0.81–0.92). Activity traits were low to moderately heritable (22–35%) and showed low genetic correlations with production traits and physical abnormalities. We demonstrated a simple and cost-efficient method to extract body dimension parameters and activity traits. These traits were estimated to be heritable, and hence, can be selected on. These findings are valuable for (pig) breeding organizations, as they offer a method to automatically phenotype new production and behavioral traits on an individual level.

Screening for antibodies against the sheep scab mite (Psoroptes ovis) Pso o 2 antigen in experimentally infested Swifter sheep may fail to identify affected animals

Sheep scab, caused by Psoroptes ovis mites, represents a significant threat to sheep health and welfare. Infestations are diagnosed by parasite identification in skin scrapings, and more recently with a commercial ELISA against serum antibodies to the Pso o 2 mite allergen. However, little is known about the performance of the ELISA in non-UK sheep populations. In this study, six Swifter sheep were experimentally infested with P. ovis. Lesion sizes were monitored and serum IgG against Pso o 2 and the novel Pso-EIP-1 antigens were measured by ELISA. Although all sheep showed signs of infestation, serum from two animals failed to react with Pso o 2. However, they did react to Pso-EIP-1. This indicates that cases of sheep scab in (Swifter) sheep may remain undetected using the Pso o 2 ELISA, which may have implications for routine screening of non- UK sheep breeds.

The influence of thyroid state on hypothalamic AMP-activated protein kinase pathways in broilers

To investigate whether there are important interactions in play in broilers between thyroid hormones and the central regulation of energy homeostasis through AMP-activated protein kinase (AMPK), we induced a functional hyperthyroid and hypothyroid state in broiler chicks, and quantified systemic and hypothalamic AMPK related gene expression and related protein. Thyroid state was manipulated through dietary supplementation of triiodothyronine (T₃) or methimazole (MMI) for 7 days. A hypothalamic AMPK suppressor, 0.1% α-lipoic acid (α-LA) was used to assess the effects of the T₃ and MMI feed formulations on the AMPK pathways. Feed intake and body weight were reduced in both hypothyroid and hyperthyroid conditions. In hyperthyroid conditions (T₃ supplementation) expression of the AMPKα1 subunit increased, while in hypothyroid conditions (MMI supplementation) active phosphorylated AMPK levels in the hypothalamus dropped, but gene expression of the AMPKα1 and α2 subunit increased. For FAS and ACC (involved in fatty acid metabolism), and CRH, TRH and CNR1 (anorexigenic neuropeptides stimulating energy expenditure) there were indications that their regulation in response to thyroid state might be modulated through AMPK pathways. Our results indicate that the expression of hypothalamic AMPK as well as that of several other genes from AMPK pathways are involved in thyroid-hormone-induced changes in appetite, albeit differently according to thyroid state.

On-farm prevalence of and potential risk factors for boar taint

Boar taint is an unpleasant taste and odor that can occur in entire male pigs and is caused by androstenone, skatole, and to a lesser extent indole accumulating in fat tissue. In the present observational study, we evaluated an extensive list of such potential risk factors which influence boar taint: social hierarchy and puberty attainment, housing, health, preslaughter conditions, season, feed, carcass composition, slaughter weight or age, and breed. Details on these factors were collected by interviews with the participating farmers, observations on each farm by trained observers and farmers, as well as slaughterhouse data. Twenty-two farms (in West- and East-Flanders, ranging from 160 to 600 sows, selected on suitability) raising entire male pigs were included in the study to evaluate the link between boar taint and potential risk factors related to the farm and slaughter batch (114 slaughter batches and 16 791 entire male pigs in total). Average olfactory boar taint prevalence was 1.8 ± 0.8%. Boar taint prevalence varied also within farms up to a maximum range between slaughter batches of 9.1% which suggests an effect of factors varying between slaughter batches such as season or other variables varying between slaughter batches. Less aggressive behavior at the end of fattening as well as lower skin lesion scores at fattening as well as at slaughter could be associated with less boar taint. The same might be said for sexual behavior, though less convincingly from this study. Measures that reduce aggression and stress have therefore have the potential to lower boar taint prevalence. The same might be said for sexual behavior, though less convincingly from this study. Furthermore, boar taint prevalence was generally higher in winter than in summer, which is relevant from a planning perspective for the slaughterhouses to seek alternative markets. Finally, increased CP gave significantly lower boar taint prevalences. This may to some extent be explained by the negative association between boar taint and lean meat percentage, as increased dietary CP levels promote the carcass lean meat percentages which can then be associated with lower boar taint levels.

A publicly available repository of ROH islands reveals signatures of selection in different livestock and pet species

BACKGROUND

Runs of homozygosity (ROH) have become the state-of-the-art method for analysis of inbreeding in animal populations. Moreover, ROH are suited to detect signatures of selection via ROH islands and are used in other applications, such as genomic prediction and genome-wide association studies (GWAS). Currently, a vast amount of single nucleotide polymorphism (SNP) data is available online, but most of these data have never been used for ROH analysis. Therefore, we performed a ROH analysis on large medium-density SNP datasets in eight animal species (cat, cattle, dog, goat, horse, pig, sheep and water buffalo; 442 different populations) and make these results publicly available.

RESULTS

The results include an overview of ROH islands per population and a comparison of the incidence of these ROH islands among populations from the same species, which can assist researchers when studying other (livestock) populations or when looking for similar signatures of selection. We were able to confirm many known ROH islands, for example signatures of selection for the myostatin (MSTN) gene in sheep and horses. However, our results also included multiple other ROH islands, which are common to many populations and not identified to date (e.g. on chromosomes D4 and E2 in cats and on chromosome 6 in sheep).

CONCLUSIONS

We are confident that our repository of ROH islands is a valuable reference for future studies. The discovered ROH island regions represent a unique starting point for new studies or can be used as a reference for future studies. Furthermore, we encourage authors to add their population-specific ROH findings to our repository.

How to study runs of homozygosity using PLINK? A guide for analyzing medium density SNP data in livestock and pet species

Background

PLINK is probably the most used program for analyzing SNP genotypes and runs of homozygosity (ROH), both in human and in animal populations. The last decade, ROH analyses have become the state-of-the-art method for inbreeding assessment. In PLINK, the --homozyg function is used to perform ROH analyses and relies on several input settings. These settings can have a large impact on the outcome and default values are not always appropriate for medium density SNP array data. Guidelines for a robust and uniform ROH analysis in PLINK using medium density data are lacking, albeit these guidelines are vital for comparing different ROH studies. In this study, 8 populations of different livestock and pet species are used to demonstrate the importance of PLINK input settings. Moreover, the effects of pruning SNPs for low minor allele frequencies and linkage disequilibrium on ROH detection are shown.

Results

We introduce the genome coverage parameter to appropriately estimate F_ROH and to check the validity of ROH analyses. The effect of pruning for linkage disequilibrium and low minor allele frequencies on ROH analyses is highly population dependent and such pruning may result in missed ROH. PLINK’s minimal density requirement is crucial for medium density genotypes and if set too low, genome coverage of the ROH analysis is limited. Finally, we provide recommendations for the maximal gap, scanning window length and threshold settings.

Conclusions

In this study, we present guidelines for an adequate and robust ROH analysis in PLINK on medium density SNP data. Furthermore, we advise to report parameter settings in publications, and to validate them prior to analysis. Moreover, we encourage authors to report genome coverage to reflect the ROH analysis’ validity. Implementing these guidelines will substantially improve the overall quality and uniformity of ROH analyses.

Unraveling the genetic diversity of Belgian Milk Sheep using medium-density SNP genotypes

The present study focuses on the Belgian Milk Sheep in Flanders (Belgium) and compares its genetic diversity and relationship with the Flemish Sheep, the Friesian Milk Sheep, the French Lacaune dairy sheep and other Northern European breeds. For this study, 94 Belgian Milk Sheep, 23 Flemish Sheep and 22 Friesian Milk Sheep were genotyped with the OvineSNP50 array. In addition, 29 unregistered animals phenotypically similar to Belgian Milk Sheep were genotyped using the 15K ISGC chip. Both Belgian and Friesian Milk Sheep as well as the East Friesian Sheep were found to be less diverse than the other seven breeds included in this study. Genomic inbreeding coefficients based on runs of homozygosity (ROH) were estimated at 14.5, 12.4 and 10.2% for Belgian Milk Sheep, Flemish Sheep and Friesian Milk Sheep respectively. Out of 29 unregistered Belgian Milk Sheep, 28 mapped in the registered Belgian Milk Sheep population. Ancestry analysis, PCA and F_ST calculations showed that Belgian Milk Sheep are more related to Friesian Milk Sheep than to Flemish Sheep, which was contrary to the breeders' expectations. Consequently, breeders may prefer to crossbreed Belgian Milk Sheep with Friesian sheep populations (Friesian Milk Sheep or East Friesian Sheep) in order to increase diversity. This research underlines the usefulness of SNP chip genotyping and ROH analyses for monitoring genetic diversity and studying genetic links in small livestock populations, profiting from internationally available genotypes. As assessment of genetic diversity is vital for long-term breed survival, these results will aid flockbooks to preserve genetic diversity.

SNP genotypes reveal breed substructure, selection signatures and highly inbred regions in Piétrain pigs

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed’s genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and F_ST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population‐specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34–126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.

Inter- and intra-breed genome-wide copy number diversity in a large cohort of European equine breeds

BACKGROUND

Copy Number Variation (CNV) is a common form of genetic variation underlying animal evolution and phenotypic diversity across a wide range of species. In the mammalian genome, high frequency of CNV differentiation between breeds may be candidates for population-specific selection. However, CNV differentiation, selection and its population genetics have been poorly explored in horses.

RESULTS

We investigated the patterns, population variation and gene annotation of CNV using the Axiom® Equine Genotyping Array (670,796 SNPs) from a large cohort of individuals (N = 1755) belonging to eight European horse breeds, varying from draught horses to several warmblood populations. After quality control, 152,640 SNP CNVs (individual markers), 18,800 segment CNVs (consecutive SNP CNVs of same gain/loss state or both) and 939 CNV regions (CNVRs; overlapping segment CNVs by at least 1 bp) compared to the average signal of the reference (Belgian draught horse) were identified. Our analyses showed that Equus caballus chromosome 12 (ECA12) was the most enriched in segment CNV gains and losses (~ 3% average proportion of the genome covered), but the highest number of segment CNVs were detected on ECA1 and ECA20 (regardless of size). The Friesian horses showed private SNP CNV gains (> 20% of the samples) on ECA1 and Exmoor ponies displayed private SNP CNV losses on ECA25 (> 20% of the samples). The Warmblood cluster showed private SNP CNV gains located in ECA9 and Draught cluster showed private SNP CNV losses located in ECA7. The length of the CNVRs ranged from 1 kb to 21.3 Mb. A total of 10,612 genes were annotated within the CNVRs. The PANTHER annotation of these genes showed significantly under- and overrepresented gene ontology biological terms related to cellular processes and immunity (Bonferroni P-value < 0.05). We identified 80 CNVRs overlapping with known QTL for fertility, coat colour, conformation and temperament. We also report 67 novel CNVRs.

CONCLUSIONS

This work revealed that CNV patterns, in the genome of some European horse breeds, occurred in specific genomic regions. The results provide support to the hypothesis that high frequency private CNVs residing in genes may potentially be responsible for the diverse phenotypes seen between horse breeds.

Genomic Regions Associated with IgE Levels against Culicoides spp. Antigens in Three Horse Breeds

Insect bite hypersensitivity (IBH), which is a cutaneous allergic reaction to antigens from Culicoides spp., is the most prevalent skin disorder in horses. Misdiagnosis is possible, as IBH is usually diagnosed based on clinical signs. Our study is the first to employ IgE levels against several recombinant Culicoides spp. allergens as an objective, independent, and quantitative phenotype to improve the power to detect genetic variants that underlie IBH. Genotypes of 200 Shetland ponies, 127 Icelandic horses, and 223 Belgian Warmblood horses were analyzed while using a mixed model approach. No single-nucleotide polymorphism (SNP) passed the Bonferroni corrected significance threshold, but several regions were identified within and across breeds, which confirmed previously identified regions of interest and, in addition, identifying new regions of interest. Allergen-specific IgE levels are a continuous and objective phenotype that allow for more powerful analyses when compared to a case-control set-up, as more significant associations were obtained. However, the use of a higher density array seems necessary to fully employ the use of IgE levels as a phenotype. While these results still require validation in a large independent dataset, the use of allergen-specific IgE levels showed value as an objective and continuous phenotype that can deepen our understanding of the biology underlying IBH.

Copy number variations in Friesian horses and genetic risk factors for insect bite hypersensitivity

Background

Many common and relevant diseases affecting equine welfare have yet to be tested regarding structural variants such as copy number variations (CNVs). CNVs make up a substantial proportion of total genetic variability in populations of many species, resulting in more sequence differences between individuals than SNPs. Associations between CNVs and disease phenotypes have been established in several species, but equine CNV studies have been limited. Aim of this study was to identify CNVs and to perform a genome-wide association (GWA) study in Friesian horses to identify genomic loci associated with insect bite hypersensitivity (IBH), a common seasonal allergic dermatitis observed in many horse breeds worldwide.

Results

Genotypes were obtained using the Axiom® Equine Genotyping Array containing 670,796 SNPs. After quality control of genotypes, 15,041 CNVs and 5350 CNV regions (CNVRs) were identified in 222 Friesian horses. Coverage of the total genome by CNVRs was 11.2% with 49.2% of CNVRs containing genes. 58.0% of CNVRs were novel (i.e. so far only identified in Friesian horses). A SNP- and CNV-based GWA analysis was performed, where about half of the horses were affected by IBH. The SNP-based analysis showed a highly significant association between the MHC region on ECA20 and IBH in Friesian horses. Associations between the MHC region on ECA20 and IBH were also detected based on the CNV-based analysis. However, CNVs associated with IBH in Friesian horses were not often in close proximity to SNPs identified to be associated with IBH.

Conclusions

CNVs were identified in a large sample of the Friesian horse population, thereby contributing to our knowledge on CNVs in horses and facilitating our understanding of the equine genome and its phenotypic expression. A clear association was identified between the MHC region on ECA20 and IBH in Friesian horses based on both SNP- and CNV-based GWA studies. These results imply that MHC contributes to IBH sensitivity in Friesian horses. Although subsequent analyses are needed for verification, nucleotide differences, as well as more complex structural variations like CNVs, seem to contribute to IBH sensitivity. IBH should be considered as a common disease with a complex genomic architecture.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0657-0) contains supplementary material, which is available to authorized users.

On farm intervention studies on reduction of boar taint prevalence: Feeding strategies, presence of gilts and time in lairage

One of the challenges in the production of entire male pigs is the occurrence of boar taint. We separately tested the effect of 3 management strategies to reduce boar taint on respectively 2, 3, and 6 Flemish pig farms: 1) adapted feeding strategies, 2) presence of gilts in the compartment, and 3) varying lairage duration at the slaughterhouse. A commercialized feed concept resulted in a significant reduction of olfactory boar taint prevalence when fed for 2 weeks (T2W) compared to control (T-CON) (P = 0.030). For T2W and when fed for 3 weeks (T3W), androstenone (AND) (P = 0.002 for T2W, P = 0.029 for T3W) and skatole (SKA) (P < 0.001 for T2W and T3W) were significantly reduced compared to T-CON. Olfactory boar taint prevalence was significantly reduced when feeding 5% dried chicory roots (FI5%) (P = 0.032), but not for 3% dried chicory roots (FI3%) (P = 0.958). SKA concentration was significantly lower when feeding FI5% (P < 0.001) and when feeding FI3% (P = 0.034). Rearing entire male pigs separately from gilts and increasing lairage duration from <1 h to >3 h did not significantly affect boar taint.

Assessing the relative importance of health and conformation traits in the cavalier king Charles spaniel

Background

The selection of a future breeding dog is a complicated task, in which disease characteristics and different traits have to be combined and weighed against one another. Truncation selection, that is the exclusion of affected animals, may be very inefficient when selecting on a large number of traits, and may result in a reduction of the genetic diversity in a population or breed. Selection could be facilitated by the use of a selection index that combines multiple traits or breeding values into one score. This however requires a consideration of their relative value according to their economic weight, which is difficult to express in monetary units for health traits. The use of a choice experiment to derive non-market values might be a solution to this problem. This is a pilot study to assess the potential use of choice experiments to ascertain the public preference and relative importance attached to health- and conformation traits in the selection of a Cavalier King Charles spaniel. The focus was on two prevalent disorders, mitral valve disease and syringomyelia, and on several important conformation traits such as muzzle length and eye shape. Based on available prior information, a Bayesian D-optimal design approach was used to develop a choice experiment and the resulting choice sets.

Results

Every participant (breeder or owner) in the choice experiment was presented with a total of 17 choice sets, in which at most four traits could vary to reduce the cognitive burden. A total of 114 respondents participated in the choice experiment and results showed that respondents (breeders/owners) current attitudes were directed towards health (syringomyelia and mitral valve disease), followed by eye shape and level of inbreeding.

Conclusions

This approach identifies the value breeders and owners attach to certain traits in the breeding objective. The resulting relative weights, represented as the logworths obtained from the choice experiment, could be an alternative to economic weights. They could be implemented as a weight when breeding values are available, but more study on this topic will be necessary. A challenge in this approach is to scale up the experiment with additional traits. Moreover, for other traits, the genetic parameters and correlations should be known first, in order to include them in the health selection index as well.

Twelve years of chiari-like malformation and syringomyelia scanning in Cavalier King Charles Spaniels in the Netherlands: Towards a more precise phenotype

Chiari-like malformation (CM), syringomyelia (SM) and middle ear effusion (also called PSOM) are three conditions that frequently occur in Cavalier King Charles Spaniels (CKCS). Both CM and SM are currently screened in the Netherlands prior to breeding and are graded according to the British Veterinary Association's Kennel Club (BVA/KC) scheme. This study evaluated the prevalence and estimated genetic parameter of CM, SM and middle ear effusion from 12 years of screening results. For SM, the classical method using the BVA/KC scheme, was compared with exact measuring of the central canal dilation. For CM, the BVA/KC scheme was compared with a more detailed scheme. Next to this the presence of microchip artifacts was assessed. 1249 screening of 1020 dogs were re-evaluated. Results indicated the presence of CM in all dogs, suggesting it has become a breed-specific characteristic. And although different grades of CM were observed, the condition did not deteriorate over time. SM was present in 39% of the dogs and a clear age effect was demonstrated, with SM increasing with age. This emphasizes the importance of screening at appropriate age, since SM can worsen with increasing age. One alternative is to promote repeated measures. The presence of middle ear effusion in this study was 19%-21% for dogs younger than 3 years, and 32%-38% for dogs older than 3 years. In as much as 60%, microchip artifacts were noticed, leading to the recommendation to place microchips in another location in breeds that are susceptible to developing SM. Finally, this study estimated the heritability of CM in this population, due to the lack of phenotypic variance, to be very low at 0.02-0.03. The heritability for SM central canal dilatation to be 0.30, compared to 0.13 for the classical BVA/KC method, using a model including the age effect and the combined effect of veterinary clinic and year of the evaluation. Genetic correlations were rather small, ranging from 0.16-0.33. As a conclusion, screening for SM and CM in the entire population should be maintained, and a selection scheme against SM should be based on estimated breeding values for the exact measurement of the central canal dilatation.

A two-level pen for fattening pigs: Effects on behavior, performance, and postslaughter measurements

Concurrent with a tendency toward higher slaughter weights of fattening pigs, minimum requirements for space allowance are increasing. Allowing pigs more space in existing standard pens, however, leads to a decrease in the number of pigs per pen, which jeopardizes the economic viability of the pig farm. A possible solution includes creating a two-level pen by constructing a second level in an existing pen, to enable an increase in space allowance per pig with the same number of pigs. We investigated the effect of such a pen on the behavior, performance, and postslaughter results of fattening pigs during the entire fattening period (30 to 110 kg). A total of 444 pigs were distributed over standard control pens (0.74 m/pig) and two-level pens (0.99 m/pig). Feed was provided only in a trough at the ground level of the pen. The results show that the increased space allowance and choice of levels in two-level pens influenced the pigs' behavior in a positive way, as they performed less manipulation of pen mates ( < 0.0001). Moreover, pigs lay down more ( = 0.0007) and showed less head knocking ( = 0.005) in two-level pens, suggesting that they were calmer. In line with the behavioral results, pigs in two-level pens had fewer lesions on tails ( = 0.006), ears ( = 0.008), and shoulders ( = 0.01). Growth performance was not affected ( > 0.05), but postslaughter measurements revealed fewer livers with white spots in pigs from two-level pens. This result might be related to lower disease pressure in pens with lower stocking densities. In conclusion, a two-level pen seemed to have positive effects on the behavior of pigs (in terms of pen mate manipulation) and the increase in space allowance seemed to make the pigs calmer. The use of the second level by heavier pigs should, however, be further studied to investigate whether this extra space could be fully counted for the minimal space requirements.

Conformation Traits and Gaits in the Icelandic Horse are Associated with Genetic Variants in Myostatin (MSTN)

Many genes are known to have an influence on conformation and performance traits; however, the role of one gene, Myostatin (MSTN), has been highlighted in recent studies on horses. Myostatin acts as a repressor in the development and regulation of differentiation and proliferative growth of skeletal muscle. Several studies have examined the link between MSTN, conformation, and performance in racing breeds, but no studies have investigated the relationship in Icelandic horses. Icelandic horses, a highly unique breed, are known both for their robust and compact conformation as well as their additional gaits tölt and pace. Three SNPs (g.65868604G>T [PR8604], g.66493737C>T [PR3737], and g.66495826A>G [PR5826]) flanking or within equine MSTN were genotyped in 195 Icelandic horses. The SNPs and haplotypes were analyzed for association with official estimated breeding values (EBV) for conformation traits (n = 11) and gaits (n = 5). The EBV for neck, withers, and shoulders was significantly associated with both PR8604 and PR3737 (P < 0.05). PR8604 was also associated with EBV for total conformation (P = 0.05). These associations were all supported by the haplotype analysis. However, while SNP PR5826 showed a significant association with EBVs for leg stance and hooves (P < 0.05), haplotype analyses for these traits failed to fully support these associations. This study demonstrates the possible role of MSTN on both the form and function of horses from non-racing breeds. Further analysis of Icelandic horses as well as other non-racing breeds would be beneficial and likely help to completely understand the influence of MSTN on conformation and performance in horses.

Using an Inbred Horse Breed in a High Density Genome-Wide Scan for Genetic Risk Factors of Insect Bite Hypersensitivity (IBH)

While susceptibility to hypersensitive reactions is a common problem amongst humans and animals alike, the population structure of certain animal species and breeds provides a more advantageous route to better understanding the biology underpinning these conditions. The current study uses Exmoor ponies, a highly inbred breed of horse known to frequently suffer from insect bite hypersensitivity, to identify genomic regions associated with a type I and type IV hypersensitive reaction. A total of 110 cases and 170 controls were genotyped on the 670K Axiom Equine Genotyping Array. Quality control resulted in 452,457 SNPs and 268 individuals being tested for association. Genome-wide association analyses were performed using the GenABEL package in R and resulted in the identification of two regions of interest on Chromosome 8. The first region contained the most significant SNP identified, which was located in an intron of the DCC netrin 1 receptor gene. The second region identified contained multiple top SNPs and encompassed the PIGN, KIAA1468, TNFRSF11A, ZCCHC2, and PHLPP1 genes. Although additional studies will be needed to validate the importance of these regions in horses and the relevance of these regions in other species, the knowledge gained from the current study has the potential to be a step forward in unraveling the complex nature of hypersensitive reactions.

Differential Expression of Genes and DNA Methylation associated with Prenatal Protein Undernutrition by Albumen Removal in an avian model

Previously, long-term effects on body weight and reproductive performance have been demonstrated in the chicken model of prenatal protein undernutrition by albumen removal. Introduction of such persistent alterations in phenotype suggests stable changes in gene expression. Therefore, a genome-wide screening of the hepatic transcriptome by RNA-Seq was performed in adult hens. The albumen-deprived hens were created by partial removal of the albumen from eggs and replacement with saline early during embryonic development. Results were compared to sham-manipulated hens and non-manipulated hens. Grouping of the differentially expressed (DE) genes according to biological functions revealed the involvement of processes such as ‘embryonic and organismal development’ and ‘reproductive system development and function’. Molecular pathways that were altered were ‘amino acid metabolism’, ‘carbohydrate metabolism’ and ‘protein synthesis’. Three key central genes interacting with many DE genes were identified: UBC, NR3C1, and ELAVL1. The DNA methylation of 9 DE genes and 3 key central genes was examined by MeDIP-qPCR. The DNA methylation of a fragment (UBC_3) of the UBC gene was increased in the albumen-deprived hens compared to the non-manipulated hens. In conclusion, these results demonstrated that prenatal protein undernutrition by albumen removal leads to long-term alterations of the hepatic transcriptome in the chicken.

Genetic parameters and estimated breeding values of insect bite hypersensitivity in Belgian Warmblood horses

Genetic factors involved in susceptibility to insect bite hypersensitivity (IBH) in Belgian Warmblood horses (BWP) were investigated. Data relating to 3409 horses were collected using a questionnaire, administered to owners during sport competitions, BWP breeding days, breeder visits and after phone calls. Horses were classified as IBH-affected or unaffected, based on two 'disease classifiers': a lifetime record, based on owner information (life_status) and another based on whether or not the horse was showing clinical signs at the time of questioning (clin_status). IBH prevalence was 10% based on life_status, and 6.2% based on clin_status. The heritabilities estimated using threshold animal models varied from 0.65 to 0.78 on the underlying scale (0.18-0.26 on the observed scale). These research findings indicate that susceptibility to IBH is a heritable trait in BWP.

An intervention study demonstrates effects of MC4R genotype on boar taint and performances of growing-finishing pigs

The Asp298Asn polymorphism of the melanocortin-4 receptor (MC4R) in pigs is known to affect economically important traits such as growth rate and backfat thickness. We have assessed the possible use of this polymorphism as a molecular marker to perform genetic selection toward lower boar taint levels without compromising growth performance and carcass and meat quality in commercial boars and gilts. Homozygous boars and gilts of the AA genotype and GG genotype were compared in an intervention study with a 2 × 2 design to assess main effects and possible interactions between sex and genotype. The concentrations of the 3 boar taint compounds androstenone ( = 0.044), skatole ( = 0.049), and indole ( = 0.006) were significantly higher in fat of AA boars compared to GG boars. However, no effect on the sensory analysis of the fat samples could be observed. Between 20 and 115 kg BW, AA pigs showed higher ADFI than GG pigs ( < 0.001). An interaction between genotype and sex was observed for ADG ( = 0.044): AA boars had a significantly higher ADG than GG boars but there was no significant difference between the gilts. Daily lean meat gain tended to be higher in boars compared to gilts ( = 0.051), independent of genotype. Similarly, boars showed higher G:F compared to gilts ( < 0.001), without effect of genotype. Genotype and sex affected several carcass quality parameters but there was no interaction. Pigs of the AA genotype displayed a lower dressing percentage ( = 0.005), lower ham width ( = 0.024), lower muscle thickness ( = 0.011), and higher fat thickness ( < 0.001), resulting in a lower lean meat percentage ( < 0.001) in comparison with GG pigs. Gilts had a significantly higher dressing percentage ( < 0.001), higher muscle thickness ( < 0.001), higher ham width ( < 0.001), and lower ham angle ( < 0.001) compared to boars. Other than the boar taint compounds, meat quality was not affected by genotype. Pork of gilts was darker ( = 0.014) and less exudative during cooking ( < 0.001) and contained more intramuscular fat ( = 0.013). These results indicate that genetic selection against boar taint is possible using this marker. This will also result in lower feed intake and ADG and, consequently, better carcass quality.

On the use of the transmission disequilibrium test to detect pseudo-autosomal variants affecting traits with sex-limited expression

We herein describe the realization of a genome-wide association study for scrotal hernia and cryptorchidism in Norwegian and Belgian commercial pig populations. We have used the transmission disequilibrium test to avoid spurious associations due to population stratification. By doing so, we obtained genome-wide significant signals for both diseases with SNPs located in the pseudo-autosomal region in the vicinity of the pseudo-autosomal boundary. By further analyzing these signals, we demonstrate that the observed transmission disequilibria are artifactual. We determine that transmission bias at pseudo-autosomal markers will occur (i) when analyzing traits with sex-limited expression and (ii) when the allelic frequencies at the marker locus differ between X and Y chromosomes. We show that the bias is due to the fact that (i) sires will preferentially transmit the allele enriched on the Y (respectively X) chromosome to affected sons (respectively daughters) and (ii) dams will appear to preferentially transmit the allele enriched on the Y (respectively X) to affected sons (respectively daughters), as offspring inheriting the other allele are more likely to be non-informative. We define the conditions to mitigate these issues, namely by (i) extracting information from maternal meiosis only and (ii) ignoring trios for which sire and dam have the same heterozygous genotype. We show that by applying these rules to scrotal hernia and cryptorchidism, the pseudo-autosomal signals disappear, confirming their spurious nature.

Development of a quantitative method for the simultaneous analysis of the boar taint compounds androstenone, skatole and indole in porcine serum and plasma by means of ultra-high performance liquid chromatography coupled to high resolution mass spectrometry

Boar taint is an off-odour occurring while heating meat or fat from boars. A method detecting the three compounds (androstenone, skatole and indole) simultaneously in blood would offer substantial advantages since it would allow monitoring the impact of rearing strategies. Therefore, a UHPLC-HR-Orbitrap-MS analysis method is optimized and validated for the quantification of these compounds in plasma or serum. Sample pre-treatment involved an extraction with diethylether followed by a centrifugal filtration (30 kDa). Limits of detection and quantification varied between 0.5 and 1 μg L(-1) and 2 and 3 μg L(-1) for the three compounds, respectively. Besides, an excellent repeatability (RSD < 7.6%), within-laboratory reproducibility (RSD<10.5%), recovery (87-97%) and linearity (R(2)>0.99) were recorded. Correlations between serum/plasma and fat levels of the boar taint compounds were positive for skatole (r(serum) = 0.39 and r(plasma) = 0.84) and androstenone (r(serum) = 0.73-0.78 and r(plasma) = 0.32-0.80).

Art meets science: The Cosmopolitan Chicken Research Project

The Cosmopolitan Chicken Project is an artistic undertaking of renowned artist Koen Vanmechelen. In this project, the artist interbreeds domestic chickens from different countries aiming at the creation of a true Cosmopolitan Chicken as a symbol for global diversity. The unifying theme is the chicken and the egg, symbols that link scientific, political, philosophical and ethical issues. The Cosmopolitan Chicken Research Project is the scientific component of this artwork. Based on state of the art genomic techniques, the project studies the effect of the crossing of chickens on the genetic diversity. Also, this research is potentially applicable to the human population. The setup of the CC®P is quite different from traditional breeding experiments: starting from the crossbreed of two purebred chickens (Mechelse Koekoek x Poule de Bresse), every generation is crossed with a few animals from another breed. For 26 of these purebred and crossbred populations, genetic diversity was measured (1) under the assumption that populations were sufficiently large to maintain all informative SNP within a generation and (2) under the circumstances of the CCP breeding experiment. Under the first assumption, a steady increase in genetic diversity was witnessed over the consecutive generations, thus indeed indicating the creation of a "Cosmopolitan Chicken Genome". However, under the conditions of the CCP, which reflects the reality within the human population, diversity is seen to fluctuate within given boundaries instead of steadily increasing. A reflection on this might be that this is because, in humans, an evolutionary optimum in genetic diversity is reached. Key words.

Insect bite hypersensitivity in Belgian warmblood horses: prevalence and risk factors

Insect bite hypersensitivity (IBH) is an allergic reaction to the bites of certain Culicoides spp. or other insects. In this study, risk factors for IBH in Belgian warmblood horses stabled or grazing in Flanders (Belgium) were investigated. IBH records (n=3409) were collected in 2009 and 2011 using a questionnaire and face-to-face interviews. The classification of IBH-affected versus unaffected horses was based on the owner’s statement, and the reported IBH lifetime prevalence was 10%. Thirty eight percent of IBH affected horses had no clinical symptoms at the time of questioning. When only the presence or absence of clinical symptoms at the time of questioning was taken into account, the prevalence of IBH symptoms was 6.2%. Seventy percent of IBH-affected horses were treated with IBH measures to reduce clinical symptoms. Model selection was based on backwards elimination in a logistic regression framework starting with 17 factors. The age of the horse, vegetation of surrounding pasture and stud size were found to be significantly associated with the self-reported IBH status.

The MC4R c.893G>A mutation: a marker for growth and leanness associated with boar taint odour in Belgian pig breeds

Since surgical castration of male piglets without anaesthesia is under heavy societal pressure, finding a sustainable solution to reduce boar taint has become urgent. One way to circumvent this animal welfare violation is raising entire male pigs whilst selecting against the tainted phenotype through marker-assisted selection. Since slaughtering at a lower weight is often suggested to reduce boar taint, selection using a marker for that trait could be a promising strategy. Therefore, in this study a melanocortin 4 receptor (MC4R) mutation, frequently described in different pig breeds as marker for fat content, weight gain and feed intake, was examined in relation to boar taint in pig breeds used in Belgian pig farms. Although results suggest an association between this mutation and a boar taint odour score assigned by experts, no association was found between the mutation and the concentration of the individual chemical boar taint components androstenone, skatole and indole.

Chronic progressive lymphedema in the Belgian draft horse in Belgium: clinical phenotyping, prevalence and risk factor analysis

Chronic progressive lymphedema (CPL) is a well-known problem in the Belgian draft horse. However, tools for clinical diagnosis have not been standardized. Therefore, the exact disease prevalence and the possible CPL risk factors have not been investigated in the Belgian breed yet. The aim of the present study was to assess a uniform CPL clinical examination method to increase diagnostic objectivity as well as to determine the occurrence and severity of lesions. Using univariate and stepwise multivariable logistic regression models, the association of factors, such as age-gender interaction, coat color and season, with CPL occurrence was examined. In this study, it was demonstrated that CPL is highly prevalent in the Belgian draft horse: 60.66% of the horses of the total sample were affected, including a large proportion of young horses (< 3 years), whereas 85.86% was affected in a subset of older horses (≥ 3 years). In the latter, the lesions were more explicit as CPL is a chronic disorder. In some horses (i.e. 14% of the yearlings), mild clinical symptoms appeared at approximately the age of one, while distinct onset of the disease occurred from the age of three onwards. The factors age-gender interaction, coat color and season, which are significantly associated with CPL occurrence, were identified and quantified. In this study, a standardized diagnostic protocol is proposed that will enhance future data collection and furthermore will offer a foundation for quantitative genetic research. Ultimately, it will help to reduce CPL occurrence in the Belgian draft horse by means of selection.