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.

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.

BAMBI: A new method for automated assessment of bidirectional early-life interaction between maternal behavior and pup vocalization in mouse dam-pup dyads

Vital early-life dyadic interaction in mice requires a pup to signal its needs adequately, and a dam to recognize and respond to the pup's cues accurately and timely. Previous research might have missed important biological and/or environmental elements of this complex bidirectional interaction, because it often focused on one dyadic member only. In laboratory rodents, the Pup Retrieval Test (PRT) is the leading procedure to assess pup-directed maternal care. The present study describes BAMBI (Bidirectional Automated Mother-pup Behavioral Interaction test), a novel automated PRT methodology based on synchronous video recording of maternal behavior and audio recording of pup vocalizations, which allows to assess bidirectional dam-pup dyadic interaction. We were able to estimate pup retrieval and pup vocalization parameters accurately in 156 pups from 29 dams on postnatal days (PND) 5, 7, 9, 11, and 13. Moreover, we showed an association between number of emitted USVs and retrieval success, indicating dyadic interdependency and bidirectionality. BAMBI is a promising new automated home-cage behavioral method that can be applied to both basic and preclinical studies investigating complex phenotypes related to early-life social development.

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.

Automated procedure to assess pup retrieval in laboratory mice

All mammalian mothers form some sort of caring bond with their infants that is crucial to the development of their offspring. The Pup Retrieval Test (PRT) is the leading procedure to assess pup-directed maternal care in laboratory rodents, used in a wide range of basic and preclinical research applications. Most PRT protocols require manual scoring, which is prone to bias and spatial and temporal inaccuracies. This study proposes a novel procedure using machine learning algorithms to enable reliable assessment of PRT performance. Automated tracking of a dam and one pup was established in DeepLabCut and was combined with automated behavioral classification of "maternal approach", "carrying" and "digging" in Simple Behavioral Analysis (SimBA). Our automated procedure estimated retrieval success with an accuracy of 86.7%, whereas accuracies of "approach", "carry" and "digging" were estimated at respectively 99.3%, 98.6% and 85.0%. We provide an open-source, step-by-step protocol for automated PRT assessment, which aims to increase reproducibility and reliability, and can be easily shared and distributed.

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.