Population structure and genomic diversity of the Einsiedler horse

The breeding history of the Einsiedler horse is closely connected with the Benedictine cloister Einsiedeln. In the mid-nineteenth century, it was decided to use European Warmblood stallions for cross-breeding and to abandon the selection of stallions. Since that time, it has only been possible to trace back the origin of Einsiedler horses using maternal ancestry information. Here, we collected high-density genotype data for European Warmblood horses (Selle Français, Swiss Warmblood and Einsiedler) and Franches-Montagnes horses, the last native Swiss horse breed, to unravel the current population structure of the Einsiedler horse. Using commonly applied methods to ascertain fine-scale population structures, it was not possible to clearly differentiate the Einsiedler from other European Warmblood horses. However, by means of runs of homozygosity (ROH) we were able to detect breed-specific ROH islands for the Einsiedler horse, including genes involved in domestication and adaptation to high altitude. Therefore, future breeding activities should involve the screening of these breed-specific ROH segments, the revival of cryopreserved sperm and the selection of Einsiedler stallions.

An improved transmissibility model to detect transgenerational transmitted environmental effects

BACKGROUND

Evolutionary studies have reported that non-genetic information can be inherited across generations (epigenetic marks, microbiota, cultural inheritance). Non-genetic information is considered to be a key element to explain the adaptation of wild species to environmental constraints because it lies at the root of the transgenerational transmission of environmental effects. The "transmissibility model" was proposed several years ago to better predict the transmissible potential of each animal by taking these diverse sources of inheritance into account in a global transmissible potential. We propose to improve this model to account for the influence of the environment on the global transmissible potential as well. This extension of the transmissibility model is the "transmissibility model with environment" that considers a covariance between transmissibility samplings of animals sharing the same environment. The null hypothesis of "no transmitted environmental effect" can be tested by comparing the two models using a likelihood ratio test (LRT).

RESULTS

We performed simulations that mimicked an experimental design consisting of two lines of animals with one exposed to a particular environment at a given generation. This enabled us to evaluate the performances of the transmissibility model with environment so as to detect and quantify transgenerational transmitted environmental effects. The power and the realized type I error of the LRT were compared to those of a T-test comparing the phenotype of the two lines, three generations after the environmental exposure for different sets of parameters. The power of the LRT ranged from 45 to 94%, whereas that of the T-test was always lower than 26%. In addition, the realized type I error of the T-test was 15% and that of the LRT was 5%, as expected. Variances, the covariance between transmissibility samplings, and path coefficients of transmission estimated with the transmissibility model with environment were close to their true values for all sets of parameters.

CONCLUSIONS

The transmissibility model with environment is effective in modeling vertical transmission of environmental effects.

Genetic analysis of geometric morphometric 3D visuals of French jumping horses

BACKGROUND

For centuries, morphology has been the most commonly selected trait in horses. A 3D video recording enabled us to obtain the coordinates of 43 anatomical landmarks of 2089 jumping horses. Generalized Procrustes analysis provided centered and scaled coordinates that were independent of volume, i.e., centroid size. Genetic analysis of these coordinates (mixed model; 17,994 horses in the pedigree) allowed us to estimate a variance-covariance matrix. New phenotypes were then defined: the "summarized shapes". They were obtained by linear combinations of Procrustes coordinates with, as coefficients, the eigenvectors of the genetic variance-covariance matrix. These new phenotypes were used in genome-wide association analyses (GWAS) and multitrait genetic analysis that included judges' scores and competition results of the horses.

RESULTS

We defined ten shapes that represented 86% of the variance, with heritabilities ranging from 0.14 to 0.42. Only one of the shapes was found to be genetically correlated with competition success (rg = - 0.12, standard error = 0.07). Positive and negative genetic correlations between judges' scores and shapes were found. This means that the breeding objective defined by judges involves improvement of anatomical parts of the body that are negatively correlated with each other. Known single nucleotide polymorphisms (SNPs) on chromosomes 1 and 3 for height at withers were significant for centroid size but not for any of the shapes. As these SNPs were not associated with the shape that distinguished rectangular horses from square horses (with height at withers greater than body length), we hypothesize that these SNPs play a role in the overall development of horses, i.e. in height, width, and length but not in height at withers when standardized to unit centroid size. Several other SNPs were found significant for other shapes.

CONCLUSIONS

The main application of 3D morphometric analysis is the ability to define the estimated breeding value (EBV) of a sire based on the shape of its potential progeny, which is easier for breeders to visualize in a single synthetic image than a full description based on linear profiling. However, the acceptance of these new phenotypes by breeders and the complex nature of summarized shapes may be challenging. Due to the low genetic correlations of the summarized shapes with jumping performance, the methodology did not allow indirect performance selection criteria to be defined.

Copula miss-specification in REML multivariate genetic animal model estimation

BACKGROUND

In animal genetics, linear mixed models are used to deal with genetic and environmental effects. The variance and covariance terms of these models are usually estimated by restricted maximum likelihood (REML), which provides unbiased estimators. A strong hypothesis of REML estimation is the multi-normality of the response variables. However, in practice, even if the marginal distributions of each phenotype are normal, the multi-normality assumption may be violated by non-normality of the cross-sectional dependence structure, that is to say when the copula of the multivariate distribution is not Gaussian. This study uses simulations to evaluate the impact of copula miss-specification in a bivariate animal model on REML estimations of variance components.

RESULT

Bivariate phenotypes were simulated for populations undergoing selection, considering different copulas for the dependence structure between the error components. Two multi-trait situations were considered: two phenotypes were measured on the selection candidates, or only one phenotype was measured on the selection candidates. Three generations with random selection and five generations with truncation selection based on estimated breeding values were simulated. When selection was performed at random, no significant differences were observed between the REML estimations of variance components and the true parameters even for the non-Gaussian distributions. For the truncation selections, when two phenotypes were measured on candidates, biases were systematically observed in the variance components for high residual dependence in the case of non-Gaussian distributions, especially in the case of a heavy-tailed or asymmetric distribution when the two traits were measured. Conversely, when only one phenotype was measured on candidates, no difference was observed between the Gaussian and non-Gaussian distributions in REML estimations.

CONCLUSIONS

This study confirms that REML can be used by geneticists to evaluate breeding values in the multivariate case even if the multivariate phenotypes deviate from normality in the situation of random selection or if one trait is not measured for the candidate under selection. Nevertheless, when the two traits are measured, the violation of the normality assumption may lead to non-negligible biases in the REML estimations of the variance-covariance components.

Quality of breeding value predictions from longitudinal analyses, with application to residual feed intake in pigs

Background

An important goal in animal breeding is to improve longitudinal traits. The objective of this study was to explore for longitudinal residual feed intake (RFI) data, which estimated breeding value (EBV), or combination of EBV, to use in a breeding program. Linear combinations of EBV (summarized breeding values, SBV) or phenotypes (summarized phenotypes) derived from the eigenvectors of the genetic covariance matrix over time were considered, and the linear regression method (LR method) was used to facilitate the evaluation of their prediction accuracy.

Results

Weekly feed intake, average daily gain, metabolic body weight, and backfat thickness measured on 2435 growing French Large White pigs over a 10-week period were analysed using a random regression model. In this population, the 544 dams of the phenotyped animals were genotyped. These dams did not have own phenotypes. The quality of the predictions of SBV and breeding values from summarized phenotypes of these females was evaluated. On average, predictions of SBV at the time of selection were unbiased, slightly over-dispersed and less accurate than those obtained with additional phenotypic information. The use of genomic information did not improve the quality of predictions. The use of summarized instead of longitudinal phenotypes resulted in predictions of breeding values of similar quality.

Conclusions

For practical selection on longitudinal data, the results obtained with this specific design suggest that the use of summarized phenotypes could facilitate routine genetic evaluation of longitudinal traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00722-w.

Equine vitiligo-like depigmentation in grey horses is related to genes involved in immune response and tumor metastasis

Background

In horses, the autoimmune disease vitiligo is characterized by the loss of melanocytes and results in patchy depigmentation of the skin around the eyes, muzzle and the perianal region. Vitiligo-like depigmentation occurs predominantly in horses displaying the grey coat colour and is observed at a prevalence level of 26.0–67.0% in grey horses compared with only 0.8–3.5% in non-grey horses. While the polygenetic background of this complex disease is well documented in humans, the underlying candidate genes for this skin disorder in horses remain unknown. In this study we aim to perform a genome-wide association study (GWAS) for identifying putative candidate loci for vitiligo-like depigmentation in horses.

Methods

In the current study, we performed a GWAS analysis using high-density 670 k single nucleotide polymorphism (SNP) data from 152 Lipizzan and 104 Noriker horses, which were phenotyped for vitiligo-like depigmentation by visual inspection. After quality control 376,219 SNPs remained for analyses, the genome-wide Bonferroni corrected significance level was p < 1.33e-7.

Results

We identified seven candidate genes on four chromosomes (ECA1, ECA13, ECA17, ECA20) putatively involved in vitiligo pathogenesis in grey horses. The highlighted genes PHF11, SETDB2, CARHSP1 and LITAFD, are associated with the innate immune system, while the genes RCBTB1, LITAFD, NUBPL, PTP4A1, play a role in tumor suppression and metastasis. The antagonistic pathogenesis of vitiligo in relation to cancer specific enhanced cell motility and/or metastasis on typical melanoma predilection sites underlines a plausible involvement of RCBTB1, LITAFD, NUBPL, and PTP4A1.

Conclusions

The proposed candidate genes for equine vitiligo-like depigmentation, indicate an antagonistic relation between vitiligo and tumor metastasis in a horse population with higher incidence of melanoma. Further replication and expression studies should lead to a better understanding of this skin disorder in horses.

Genomic Correlations Between the Gaits of Young Horses Measured by Accelerometry and Functional Longevity in Jumping Competition

Functional longevity is essential for the well-being of horses and the satisfaction of riders. Conventional selection using longevity breeding values calculated from competition results is not efficient because it takes too long to obtain reliable information. Therefore, the objective was to identify early criteria for selection. We assessed two types of early criteria: gait traits of young horses and QTLs. Thus, our aim was to estimate the genetic correlation between gait traits and longevity and to perform a genome-wide association study (GWAS) for longevity. Measurements of gaits by accelerometry were recorded on 1,477 show jumping horses that were 4 to 5 years old. Gait analysis provided 9 principal components describing trot, canter, and walk. Longevity estimated breeding values (EBVs) for stallions were calculated using a survival analysis of more than 900,000 years of performances by 179,448 show jumping horses born from 1981 onwards. Longevity was measured as the number of years spent in competition. Model included region and month of birth, age at first competition, year, and performance level. Longevity EBVs were deregressed to obtain weighted pseudo-performances for 1,968 stallions. Genomic data were available for 3,658 jumping horses. Seventy-eight percent of the horses measured for gaits and twenty-five percent of those measured for longevity were genotyped. A GWAS of longevity revealed no significant QTLs. Genetic parameters between each of the 9 principal components of the gait variables and longevity were evaluated with a bi-trait animal linear mixed model using single-step GBLUP analysis with the relationship matrix constructed from genomic data and genealogy (24,448 ancestors over four generations). The heritability of the gait traits varied from 0.11 to 0.44. The third principal component for trot (high lateral activity) and the first principal component for canter (high dorsoventral activity and low stride frequency) were moderately genetically correlated with higher longevity: r_g = 0.38 (0.15) and 0.28 (0.13), respectively. Our study revealed that functional longevity is a polygenic trait with no major genes. We found new correlations between longevity and gait traits. Before using gait characteristics in a selection plan, these correlations need to be understood better at the biomechanical level.

Genomic analysis of gaits and racing performance of the French trotter

The aim was to disentangle gait characteristics from other qualities needed for racing performances with a genomic analysis of French trotters (FT). A sample of 1,390 horses were recruited, from which 46% were genotyped with Illumina chip of 54,602 SNPs, 49% with Affymetrix chip of 670,806 SNPs and 586 had a completed questionnaire on trotting technique. Racing performances cover the period 1996 to 2018. There were 252,368 FT-born; 96,617 qualified and 83,962 which participated in a race. After quality control, 377,611 SNPs were retained and imputed. Questionnaire described trotting technique over 13 questions which were summarized, after principal component analysis in 3 traits: pacer, heavy trot/gallop and other defects. GWAS and genomic evaluation were performed using single-step approach. We found 25 QTL for racing performances and 9 for trotting technique. Only DMRT3 mutation was significant for both traits. To tend to pace avoid the defect at gallop and lead to a better early career for earnings, less percentage of disqualified races at all ages and more harness than under saddle career. This is the portrait of AA genotype at DMRT3. We found 5 other QTL, not linked to gait traits, which might improve selection of genetically independent performance traits of earnings per races and percentage of finished races. For only earnings at different ages and in under saddle or harness races, genomic evaluation remains the best way to predict performances.

Inclusive inheritance for residual feed intake in pigs and rabbits

Non‐genetic information (epigenetic, microbiota, behaviour) that results in different phenotypes in animals can be transmitted from one generation to the next and thus is potentially involved in the inheritance of traits. However, in livestock species, animals are selected based on genetic inheritance only. The objective of the present study was to determine whether non‐genetic inherited effects play a role in the inheritance of residual feed intake (RFI) in two species: pigs and rabbits. If so, the path coefficients of the information transmitted from sire and dam to offspring would differ from the expected transmission factor of 0.5 that occurs if inherited information is of genetic origin only. Two pigs (pig1, pig2) and two rabbits (rabbit1, rabbit2) datasets were used in this study (1,603, 3,901, 5,213 and 4,584 records, respectively). The test of the path coefficients to 0.5 was performed for each dataset using likelihood ratio tests (null model: transmissibility model with both path coefficients equal to 0.5, full model: unconstrained transmissibility model). The path coefficients differed significantly from 0.5 for one of the pig datasets (pig2). Although not significant, we observed, as a general trend, that sire path coefficients of transmission were lower than dam path coefficients in three of the datasets (0.46 vs 0.53 for pig1, 0.39 vs 0.44 for pig2 and 0.38 vs 0.50 for rabbit1). These results suggest that phenomena other than genetic sources of inheritance explain the phenotypic resemblance between relatives for RFI, with a higher transmission from the dam's side than from the sire's side.

Accelerometers Provide Early Genetic Selection Criteria for Jumping Horses

The aim of this study was to evaluate the genetic component of the locomotor jumping ability, via a wearable accelerometer sensor, and to estimate the genetic correlation with performance in competition, to introduce such criteria in selection schema. A sample of 1,056 young 3-year-old horses were equipped with a 3-dimensional accelerometer during a free jumping test, in regular breeding shows from 2015 to 2017. Seven variables were extracted from the dorso-ventral acceleration curve for the last three jumps over a double bar obstacle of 1.15 m for the front pole and 1.20 m for the back pole with a 1.20 m spread. Variables were the peaks of forelimbs, hindlimbs, and landing acceleration, the duration between peaks at take-off, the peak of forelimb acceleration and start of jump, jump duration and duration between the beginning of the impact of forelimbs and the peak at landing. During breeding shows, judges scored balance, strength, style, and reactivity for free jumping and jumping tests under saddle. Jumping competition results were recorded by logarithm of the sum of points earned in each competition. All horses in official competitions were included, i.e., 160,257 horses born in 1997 with a total of 649,491 annual performances. An animal mixed model with complete pedigree over four generations (353,236 horses) were used with fixed effects of jumping test location and date, morning/afternoon, gender, month of birth, rank of jump for accelerometric data, effect of year of competition, combined with age and gender for competition results. As a result, jump duration was the most heritable and repeatable for jump variables: h ² = 0.16 (0.06), r = 0.52 (0.02), while accelerations were moderately heritable (h ² = 0.05-0.09, r = 0.39-0.51). Judgement scores were heritable: 0.21 (0.07)-0.33 (0.09) and were highly correlated. Scores during free jumping were genetically correlated to jump duration: 0.71 (0.15)-0.88 (0.16). Both jump duration and judgement scores were genetically correlated to competition performance: 0.59 (0.13) for jump duration, from 0.60 (0.11) to 0.77 (0.12) for scores. Jump duration and judgement scores can be used as early selection criteria. The advantage of the accelerometric measurement is its objectivity and the ease of recording.

A Unified Model for Inclusive Inheritance in Livestock Species

For years, animal selection in livestock species has been performed by selecting animals based on genetic inheritance. However, evolutionary studies have reported that nongenetic information that drives natural selection can also be inherited across generations (epigenetic, microbiota, environmental inheritance). In response to this finding, the concept of inclusive heritability, which combines all sources of information inherited across generations, was developed. To better predict the transmissible potential of each animal by taking into account these diverse sources of inheritance and improve selection in livestock species, we propose the "transmissibility model." Similarly to the animal model, this model uses pedigree and phenotypic information to estimate variance components and predict the transmissible potential of an individual, but differs by estimating the path coefficients of inherited information from parent to offspring instead of using a set value of 0.5 for both the sire and the dam (additive genetic relationship matrix). We demonstrated the structural identifiability of the transmissibility model, and performed a practical identifiability and power study of the model. We also performed simulations to compare the performances of the animal and transmissibility models for estimating the covariances between relatives and predicting the transmissible potential under different combinations of sources of inheritance. The transmissibility model provided similar results to the animal model when inheritance was of genetic origin only, but outperformed the animal model for estimating the covariances between relatives and predicting the transmissible potential when the proportion of inheritance of nongenetic origin was high or when the sire and dam path coefficients were very different.

Genome-Wide Homozygosity Patterns and Evidence for Selection in a Set of European and Near Eastern Horse Breeds

Intensive artificial and natural selection have shaped substantial variation among European horse breeds. Whereas most equine selection signature studies employ divergent genetic population structures in order to derive specific inter-breed targets of selection, we screened a total of 1476 horses originating from 12 breeds for the loss of genetic diversity by runs of homozygosity (ROH) utilizing a 670,000 single nucleotide polymorphism (SNP) genotyping array. Overlapping homozygous regions (ROH islands) indicating signatures of selection were identified by breed and similarities/dissimilarities between populations were evaluated. In the entire dataset, 180 ROH islands were identified, whilst 100 islands were breed specific, all other overlapped in 36 genomic regions with at least one ROH island of another breed. Furthermore, two ROH hot spots were determined at horse chromosome 3 (ECA3) and ECA11. Besides the confirmation of previously documented target genes involved in selection for coat color (MC1R, STX17, ASIP), body size (LCORL/NCAPG, ZFAT, LASP1, HMGA2), racing ability (PPARGC1A), behavioral traits (GRIN2B, NTM/OPCML) and gait patterns (DMRT3), several putative target genes related to embryonic morphogenesis (HOXB), energy metabolism (IGFBP-1, IGFBP-3), hair follicle morphogenesis (KRT25, KRT27, INTU) and autophagy (RALB) were highlighted. Furthermore, genes were pinpointed which might be involved in environmental adaptation of specific habitats (UVSSA, STXBP4, COX11, HLF, MMD).

89 Effect of embryo transfer and recipient breed on offspring performance in equine show jumping

Maternal preconceptional, gestational, and postnatal environments alter offspring phenotype in humans and animals. In particular, embryo technologies were shown in horses and other species to affect offspring growth and metabolism. Quantitative genetic calculations in a database containing genealogy and a measured trait (performance) of individuals enable geneticists to estimate the respective contribution of (A) genetic aptitude for the trait (inherited from all parents), (B) genetic maternal aptitude (additional contribution of the mother and her ancestors, acting through pregnancy, lactation, and nursing), and (C) the permanent maternal environment (similarities between the offspring of a dam that are not attributed to genetics) to the variability of this trait in this population. Additionally, fixed effects can be introduced and estimated to take into account additional factors acting on individual performance. The aim of this project was to apply this method to estimate the respective contributions of A, B, and C as well as the fixed effects of maternal parity (D), age at foaling (E), foaling interval (F), embryo transfer (G), and breed of recipient mare (for ET recipients) on the adult offspring sport performance in show jumping horses. The data came from the French Horse Database established in 1976 for genealogy and performance; calculation was done using WOMBAT® software. The performance of 207,984 horses (4.5 years per horse), born to 90,438 dams (of which 4187 were ET recipients, the others being bred by AI or natural service), were studied. Performance was expressed as follows. For each event, points were distributed depending on the horse’s rank within the event and the event’s difficulty. The logarithm of the sum of points of the year has a Gaussian distribution and is then normalised to mean=0 and standard deviation=1. The genetic maternal aptitude and permanent maternal environment showed a small but significant contribution to the variance of performance (1.4 and 1.5%, respectively) compared to a direct genetic contribution of 24.4%. Within each age class of dam, offspring of primiparous mares have a handicap (−0.04 to −0.06) compared to multiparous mares. Within multiparous mares, a foaling interval of 1 year only provided a small advantage (+0.01, nonsignificant) compared to longer intervals (&gt;2 years). The effect of embryo transfer was strong and significant (+0.24), which could reflect the better management of the horses produced by embryo transfer because of their genetic and financial value rather than a physiologic effect. The performance of offspring born to draft and warmblood recipients (+0.27v. +0.24, nonsignificant) were not different. In conclusion, although jumping performance occurs several years after birth, maternal effects are still measurable and embryo transfer is associated with a strongly increased performance independent of genetic value, probably due to better management of these horses. No effect of recipient mare’s breed was found.

Novel insights into Sabino1 and splashed white coat color patterns in horses

Within the framework of genome‐wide analyses using the novel Axiom^® genotyping array, we investigated the distribution of two previously described coat color patterns, namely sabino1 (SBI), associated with the KIT gene (KI16+1037A), and splashed white, associated with the PAX3 gene (ECA6:g.11429753C>T; PAX3^C70Y), including a total of 899 horses originating from eight different breeds (Achal Theke, Purebred Arabian, Partbred Arabian, Anglo‐Arabian, Shagya Arabian, Haflinger, Lipizzan and Noriker). Based on the data we collected we were able to demonstrate that, besides Quarter horses, the PAX3^C70Y allele is also present in Noriker (seven out of 189) and Lipizzan (three out of 329) horses. The SB1 allele was present in three breeds (Haflinger, 14 out of 98; Noriker, four out of 189; Lipizzan one out of 329). Furthermore, we examined the phenotypes of SB1‐ and PAX3^C70Y‐carrier horses for their characteristic white spotting patterns. None of the SB1/sb1‐carrier horses met the criteria defining the Sabino1 pattern according to current applied protocols. From 10 heterozygous PAX3^C70Y‐carrier horses, two had nearly a splashed white phenotype. The results of this large‐scale experiment on the genetic association of white spotting patterns in horses underline the influence of gene interactions and population differences on complex traits such as Sabino1 and splashed white.

Human papillomavirus does not play a role in the Barrett esophagus: a French cohort

The role of human papillomavirus (HPV) in Barrett's esophagus (BE) has been examined but remains unclear. The purpose of the study is to dispute the connection between HPV and BE in a prospective case-control study. Biopsies were performed above and inside the Barrett's segment for BE patients and in the distal third of the esophagus for control patients for histological interpretation and for virological analysis. Biopsies for virological analysis were placed in a virus transport medium and immediately frozen in liquid nitrogen. Virological analysis involved real-time PCR using the SyBr® green protocol with modified SPF10 general primers. A total of 180 patients (119 control and 61 BE, respectively) were included. In BE patients, 31, 18, and 12 patients had, respectively, no dysplasia, low-grade dysplasia, and high grade dysplasia. Overall, nine were found to be HPV positive: five were control patients and four BE patients. HPV positive status was not associated with BE. No factors were associated with HPV, in particular the degree of BE dysplasia. HPV infection appears unlikely to be significant in the etiology of BE compared with control patients. (ClinicalTrials.gov, Number NCT02549053).

Endurance Exercise Ability in the Horse: A Trait with Complex Polygenic Determinism

Endurance horses are able to run at more than 20 km/h for 160 km (in bouts of 30-40 km). This level of performance is based on intense aerobic metabolism, effective body heat dissipation and the ability to endure painful exercise. The known heritabilities of endurance performance and exercise-related physiological traits in Arabian horses suggest that adaptation to extreme endurance exercise is influenced by genetic factors. The objective of the present genome-wide association study (GWAS) was to identify single nucleotide polymorphisms (SNPs) related to endurance racing performance in 597 Arabian horses. The performance traits studied were the total race distance, average race speed and finishing status (qualified, eliminated or retired). We used three mixed models that included a fixed allele or genotype effect and a random, polygenic effect. Quantile-quantile plots were acceptable, and the regression coefficients for actual vs. expected log10p-values ranged from 0.865 to 1.055. The GWAS revealed five significant quantitative trait loci (QTL) corresponding to 6 SNPs on chromosomes 6, 1, 7, 16, and 29 (two SNPs) with corrected p-values from 1.7 × 10-6 to 1.8 × 10-5. Annotation of these 5 QTL revealed two genes: sortilin-related VPS10-domain-containing receptor 3 (SORCS3) on chromosome 1 is involved in protein trafficking, and solute carrier family 39 member 12 (SLC39A12) on chromosome 29 is active in zinc transport and cell homeostasis. These two coding genes could be involved in neuronal tissues (CNS). The other QTL on chromosomes 6, 7, and 16 may be involved in the regulation of the gene expression through non-coding RNAs, CpG islands and transcription factor binding sites. On chromosome 6, a new candidate equine long non-coding RNA (KCNQ1OT1 ortholog: opposite antisense transcript 1 of potassium voltage-gated channel subfamily Q member 1 gene) was predicted in silico and validated by RT-qPCR in primary cultures of equine myoblasts and fibroblasts. This lncRNA could be one element of the cardiac rhythm regulation. Our GWAS revealed that equine performance during endurance races is a complex polygenic trait, and is partially governed by at least 5 QTL: two coding genes involved in neuronal tissues and three other loci with many regulatory functions such as slowing down heart rate.

Diagnostic value of ultrasonography in elbow trauma in children: Prospective study of 34 cases

INTRODUCTION

Among the various elbow injuries in children that initially have normal radiographs, a certain number of occult fractures are only diagnosed correctly after the fact, during a follow-up visit.

PURPOSE

This study evaluated the diagnostic contribution of ultrasonography in the treatment of acute elbow injuries in children and the strategic and economic impact of using this tool alongside radiography.

MATERIALS AND METHODS

During this prospective study performed between January 1 and April 1 2014, elbow ultrasonography was performed within 6 days in all children under 15 years of age with a suspected occult fracture. The ultrasonography exam looked for lipohemarthrosis, the posterior fat pad sign and cortical disruption. If no fracture was visible on ultrasonography, a removable splint was given to the patient to relieve pain, and no radiological or clinical follow-up was scheduled. The patients were contacted again at least 15 days later to determine whether an undetected fracture was present. Lastly, we evaluated the cost of treatment with and without ultrasonography in the cases where no fracture was diagnosed.

RESULTS

In 13 cases, ultrasonography revealed lipohemarthrosis and a fat pat sign, with cortical disruption also present in 11 of these cases. In two cases, the diagnosis was made based solely on the presence of lipohemarthrosis and a fat pat sign. There were seven lateral condyle fractures, two medial epicondyle fractures and two supracondylar fractures. Among the 21 patients with normal ultrasonography, no fracture was diagnosed later on. In patients without a fracture, using ultrasonography resulted in a cost savings of €29.10 per patient versus not using it.

CONCLUSION

In our study, ultrasonography is a sensitive examination for the diagnosis of occult elbow fractures in children. When the radiography and ultrasonography are both normal, the possibility of fracture can be rule out definitively, which reduces the need for immobilization, follow-up and treatment costs. The findings of this preliminary study should be validated with a larger prospective study.

Genome Wide Association Study Identifies New Loci Associated with Undesired Coat Color Phenotypes in Saanen Goats

This paper reports a quantitative genetics and genomic analysis of undesirable coat color patterns in goats. Two undesirable coat colors have routinely been recorded for the past 15 years in French Saanen goats. One fifth of Saanen females have been phenotyped “pink” (8.0%) or “pink neck” (11.5%) and consequently have not been included in the breeding program as elite animals. Heritability of the binary “pink” and “pink neck” phenotype, estimated from 103,443 females was 0.26 for “pink” and 0.21 for “pink neck”. Genome wide association studies (using haplotypes or single SNPs) were implemented using a daughter design of 810 Saanen goats sired by 9 Artificial Insemination bucks genotyped with the goatSNP50 chip. A highly significant signal (-log10pvalue = 10.2) was associated with the “pink neck” phenotype on chromosome 11, suggesting the presence of a major gene. Highly significant signals for the “pink” phenotype were found on chromosomes 5 and 13 (-log10p values of 7.2 and, 7.7 respectively). The most significant SNP on chromosome 13 was in the ASIP gene region, well known for its association with coat color phenotypes. Nine significant signals were also found for both traits. The highest signal for each trait was detected by both single SNP and haplotype approaches, whereas the smaller signals were not consistently detected by the two methods. Altogether these results demonstrated a strong genetic control of the “pink” and “pink neck” phenotypes in French Saanen goats suggesting that SNP information could be used to identify and remove undesired colored animals from the breeding program.

Does heterozygosity at the DMRT3 gene make French trotters better racers?

BACKGROUND

Recently, a mutation was discovered in the DMRT3 gene that controls pacing in horses. The mutant allele A is fixed in the American Standardbred trotter breed, while in the French trotter breed, the frequency of the wild-type allele C is still 24%. This study aimed at measuring the effect of DMRT3 genotypes on the performance of French trotters and explaining why the polymorphism still occurs in this breed. Using a mixed animal model, genetic parameters and environmental effects on performance traits were estimated from data on 173 176 French trotter races. The effect of the DMRT3 gene was then estimated by the effect of genotype at the highly linked SNP BIEC2-620109 (C-C, A-T) for 630 horses. A selection scheme that included qualification and racing performances was modeled to (1) verify if the observed superiority of heterozygous CT horses at this SNP could be explained only by selection and (2) understand why allele C has not disappeared in French trotters.

RESULTS

Heritability of racing performance traits was high for qualification test (0.56), moderate for annual earnings per finished race (0.26 to 0.31) and low for proportion of disqualified races (0.06 to 0.09). Genotype CC was always unfavorable compared to genotype TT for qualification: the probability to be qualified was 20% for CC vs. 48% for TT and earnings were -0.96 σy lower for CC than for TT. Genotype CT was also unfavorable for qualification (40%) and earnings at 3 years (-0.21 σy), but favorable for earnings at ages greater than 5 years: +0.41 σy (P = 7.10(-4)). Selection on qualification could not explain more than 19% of the difference between genotypes CC and CT in earnings at ages greater than 5 years. Only a scenario for which genotype CT has a favorable effect on the performance of horses older than 5 years could explain that the polymorphism at the DMRT3 gene still exists in the French trotter breed.

CONCLUSIONS

The use of mature horses in the French racing circuit can explain that the CA genotype is still present in the French trotter horses.

A comparison of methods for whole-genome QTL mapping using dense markers in four livestock species

Background

With dense genotyping, many choices exist for methods to detect quantitative trait loci (QTL) in livestock populations. However, no across-species study has been conducted on the performance of different methods using real data. We compared three methods that correct for relatedness either implicitly or explicitly: linkage and linkage disequilibrium haplotype-based analysis (LDLA), efficient mixed-model association (EMMA) analysis, and Bayesian whole-genome regression (BayesC). We analyzed one chromosome in each of five datasets (dairy cattle, beef cattle, sheep, horses, and pigs) using real genotypes based on dense single nucleotide polymorphisms and phenotypes. The P values corrected for multiple testing or Bayes factors greater than 150 were considered to be significant. To complete the real data study, we also simulated quantitative trait loci (QTL) for the same datasets based on the real genotypes. Several scenarios were chosen, with different QTL effects and linkage disequilibrium patterns. A pseudo-null statistical distribution was chosen to make the significance thresholds comparable across methods.

Results

For the real data, the three methods generally agreed within 1 or 2 cM for the locations of QTL regions and disagreed when no signals were significant (e.g. in pigs). For certain datasets, LDLA had more significant signals than EMMA or BayesC, but they were concentrated around the same peaks. Therefore, the three methods detected approximately the same number of QTL regions. For the simulated data, LDLA was slightly less powerful and accurate than either EMMA or BayesC but this depended strongly on how thresholds were set in the simulations.

Conclusions

All three methods performed similarly for real and simulated data. No method was clearly superior across all datasets or for any particular dataset. For computational efficiency and ease of interpretation, EMMA is recommended, but using more than one method is suggested.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0087-7) contains supplementary material, which is available to authorized users.

Estimation of genetic parameters for resistance to gastro-intestinal nematodes in pure blood Arabian horses

Equine internal parasites, mostly cyathostomins, affect both horse welfare and performance. The appearance of anthelmintic-resistant parasites creates a pressing need for optimising drenching schemes. This optimization may be achieved by identifying genetic markers associated with host susceptibility to infection and then to drench carriers of these markers. The aim of our study was to characterise the genetics of horse resistance to strongyle infection by estimating heritability of this trait in an Arabian pure blood population. A population of 789 Arabian pure blood horses from the Michałów stud farm, Poland were measured for strongyle egg excretion twice a year, over 8 years. Low repeatability values were found for faecal egg counts. Our analyses showed that less than 10% of the observed variation for strongyle faecal egg counts in this population had a genetic origin. However, additional analyses highlighted an age-dependent increase in heritability which was 0.04 (±0.02) in young horses (up to 3 years of age) but 0.21 (±0.04) in older ones. These results suggest that a significant part of the inter-individual variation has a genetic origin. This paves the way to a genomic dissection of horse-nematode interactions which might provide predictive markers of susceptibility, allowing individualised drenching schemes.

Genome-wide association study for jumping performances in French sport horses

A genome-wide association study was performed to identify single nucleotide polymorphisms (SNPs) associated with jumping performances of warmbloods in France. The 999 horses included in the study for jumping performances were sport horses [mostly Selle Français (68%), Anglo-Arabians (13%) and horses from the other European studbooks]. Horses were genotyped using the Illumina EquineSNP50 BeadChip. Of the 54,602 SNPs available on this chip, 44,424 were retained after quality testing. Phenotypes were obtained by deregressing official breeding values for jumping competitions to use all available information, that is, the performances of each horse as well as those of its relatives. Two models were used to test the effects of the genotypes on deregressed phenotypes: a single-marker mixed model and a haplotype-based mixed model (significant: P < 1E-05; suggestive: P < 1E-04). Both models included a polygenic effect to take into account familial structures. For jumping performances, one suggestive quantitative trait locus (QTL) located on chromosome 1 (BIEC2_31196 and BIEC2_31198) was detected with both models. This QTL explains 0.7% of the phenotypic variance. RYR2, a gene encoding a major calcium channel in cardiac muscle in humans and mice, is located 0.55 Mb from this potential QTL.

Is the use of formulae a reliable way to predict the accuracy of genomic selection?

We studied four formulae used to predict the accuracy of genomic selection prior to genotyping. The objectives of our study were to investigate the impact of the parameters of each formula on the values of accuracy calculated using these formulae, and to check whether the accuracies reported in the literature are in agreement with the formulae. First, we computed the marginal distribution of accuracy (by integration) for each parameter of all four formulae: heritability h(2) , reference population size T, number of markers M and number of effective segments in the genome Me . Then, we collected 145 accuracies and corresponding parameters reported in 13 publications on genomic selection (mainly in dairy cattle), and performed analysis of variance to test the differences between observed and predicted accuracy with effects of formulae and parameters. The variation of accuracy for different values of each parameter indicated that two parameters, T and Me, had a significant impact and that considerable differences existed between the formulae (mean accuracies differed by up to 0.20 point). The results of our meta-analysis showed a big formula effect on the accuracies predicted using each formula, and also a significant effect of the value obtained for Me calculated from Ne (effective population size). Each formula can therefore be demonstrated to be optimal depending on the assumption used for Me . In conclusion, no rules can be applied to predict the reliability of genomic selection using these formulae.

Computation of deregressed proofs for genomic selection when own phenotypes exist with an application in French show-jumping horses

Genomic evaluations often use as pseudo-phenotypes corrected means of progeny performances, like daughter yield deviations (DYD) in dairy species. In horse breeding, own performances are also available and performances from other relatives (as half sibs) may play an important part in the EBV because the number of progeny remains low, even for stallions. The first step for genomic selection in horses is therefore to generate pseudo-phenotypes for genomic analysis when parental or own information is considered. This work presents an easy method to compute deregressed EBV from regular pedigree-based genetic evaluations (EBV, reliabilities) to be used in genomic evaluations. The proposed methodology builds deregressed proofs so that they combine own performances (from genotyped individuals) and performances of relatives (outside of the genotyped sample). An application to show jumping horse data is presented. A sample of 908 stallions specialized in show jumping [71% Selle Français (SF), 17% foreign sport horses (FH), 13% Anglo Arab (AA)] were genotyped. Genotyping was performed using the Illumina Equine SNP50 BeadChip, and after quality tests, 44,444 SNP were retained. Two methods were used for genomic evaluation: GBLUP and BayesCπ, and 6 validation data sets were compared, chosen according to breeds SF + FH + AA or SF + FH, family structure (more than 3 half sibs), reliability of sires (>0.97) or sons (>0.72). In spite of a favorable genetic structure [linkage disequilibrium equal to 0.24 at 50 kb pairs], results showed low advantage of genomic evaluation. On the validation sample SF + FH + AA, the correlation between deregressed proofs and GBLUP or BayesCπ predictions was 0.39, 0.37, 0.51 according to the different validation data sets, compared with 0.36, 0.33, 0.53 obtained with BLUP predictions. Correlations were much lower on the SF + FH sample. Research is pursued to understand this low advantage of genomic selection and to improve the methodology for genomic evaluation in this context, which is less favorable than dairy cattle breeding.

Statistical distributions of test statistics used for quantitative trait association mapping in structured populations

Background

Spurious associations between single nucleotide polymorphisms and phenotypes are a major issue in genome-wide association studies and have led to underestimation of type 1 error rate and overestimation of the number of quantitative trait loci found. Many authors have investigated the influence of population structure on the robustness of methods by simulation. This paper is aimed at developing further the algebraic formalization of power and type 1 error rate for some of the classical statistical methods used: simple regression, two approximate methods of mixed models involving the effect of a single nucleotide polymorphism (SNP) and a random polygenic effect (GRAMMAR and FASTA) and the transmission/disequilibrium test for quantitative traits and nuclear families. Analytical formulae were derived using matrix algebra for the first and second moments of the statistical tests, assuming a true mixed model with a polygenic effect and SNP effects.

Results

The expectation and variance of the test statistics and their marginal expectations and variances according to the distribution of genotypes and estimators of variance components are given as a function of the relationship matrix and of the heritability of the polygenic effect. These formulae were used to compute type 1 error rate and power for any kind of relationship matrix between phenotyped and genotyped individuals for any level of heritability. For the regression method, type 1 error rate increased with the variability of relationships and with heritability, but decreased with the GRAMMAR method and was not affected with the FASTA and quantitative transmission/disequilibrium test methods.

Conclusions

The formulae can be easily used to provide the correct threshold of type 1 error rate and to calculate the power when designing experiments or data collection protocols. The results concerning the efficacy of each method agree with simulation results in the literature but were generalized in this work. The power of the GRAMMAR method was equal to the power of the FASTA method at the same type 1 error rate. The power of the quantitative transmission/disequilibrium test was low. In conclusion, the FASTA method, which is very close to the full mixed model, is recommended in association mapping studies.

Genome-wide association studies for osteochondrosis in French Trotter horses

A genome-wide association study for osteochondrosis (OC) in French Trotter horses was carried out to detect QTL using genotype data from the Illumina EquineSNP50 BeadChip assay. Analysis data came from 161 sire families of French Trotter horses with 525 progeny and family sizes ranging from 1 to 20. Genotypes were available for progeny (n = 525) and sires with at least 2 progeny (n = 98). Radiographic data were obtained from progeny using at least 10 views to reveal OC. All radiographic findings were described by at least 2 veterinary experts in equine orthopedics, and severity indices (scores) were assigned based on the size and location of the lesion. Traits used were a global score, the sum of all severity scores lesions (GM, quantitative measurement), and the presence or absence of OC on the fetlock (FM), hock (HM), and other sites (other). Data were analyzed using 2 mixed models including fixed effects, polygenic effects, and SNP or haplotype cluster effects. By combining results with both methods at moderate evidence of association threshold P < 5 × 10(-5), this genome-wide association study displayed 1 region for GM on the Equus caballus chromosome (ECA) 13, 2 for HM on ECA 3 and 14, and 1 for other on ECA 15. One region on ECA 3 for HM represented the most significant hit (P = 3 × 10(-6)). By comparing QTL between traits at a decreased threshold (P < 5 × 10(-4)), the 4 QTL detected for GM were associated to a QTL detected for FM or HM but never both. Another interesting result was that no QTL were found in common between HM and FM.

Genetic analysis of the longevity of French sport horses in jumping competition

The longevity of sport horses is an economic and ethical issue. The aim of this study was to analyze the duration of the competitive life of jumping horses in France to assess the potential for genetic evaluation and to propose rules of management for sport horses. Data included lifetimes spent in jumping competitions for the 209,296 horses born from 1968 onward with performances between 1972 and 2008; the data set contained 22% right-censored records. Longevity was measured in years. Discrete survival analysis included fixed effects of region of birth, month of birth, year of recording, age at first competition, interaction between sex and level of jumping performance as measured by the logarithm of earnings adjusted for sex, age, year, and random sire and maternal grand-sire effects. There were 16,668 sires and maternal grand-sires. All fixed effects were highly significant (P < 0.001). Management of the sports career had an important effect on longevity: against common belief, the younger the horse started competing, the longer it stayed in competition. For horses that started competing at an age of 6 yr, the risk of culling was 1.33-fold that of horses having started at 4 yr of age. The less success in competition, the greater was the chance for leaving competition, especially for horses without earnings. For a gelding without earnings, the risk of culling was 1.40-fold that of an average-performance gelding and 2.57-fold that of a top-rated gelding (performance at least 2 SD above the mean). Mares always had greater relative risk than geldings or stallions because they may be culled from competition to be used for breeding. The risk of culling for females was 1.45-fold that of a gelding with the same performance. The heritability of the length of competitive life was 0.10. Breeding values were predicted for sires, and 3,303 sires showed an accuracy greater than 0.60. Among these sires, 262 were used for breeding in 2008.

Validation of models for analysis of ranks in horse breeding evaluation

BACKGROUND

Ranks have been used as phenotypes in the genetic evaluation of horses for a long time through the use of earnings, normal score or raw ranks. A model, ("underlying model" of an unobservable underlying variable responsible for ranks) exists. Recently, a full Bayesian analysis using this model was developed. In addition, in reality, competitions are structured into categories according to the technical level of difficulty linked to the technical ability of horses (horses considered to be the "best" meet their peers). The aim of this article was to validate the underlying model through simulations and to propose a more appropriate model with a mixture distribution of horses in the case of a structured competition. The simulations involved 1000 horses with 10 to 50 performances per horse and 4 to 20 horses per event with unstructured and structured competitions.

RESULTS

The underlying model responsible for ranks performed well with unstructured competitions by drawing liabilities in the Gibbs sampler according to the following rule: the liability of each horse must be drawn in the interval formed by the liabilities of horses ranked before and after the particular horse. The estimated repeatability was the simulated one (0.25) and regression between estimated competing ability of horses and true ability was close to 1. Underestimations of repeatability (0.07 to 0.22) were obtained with other traditional criteria (normal score or raw ranks), but in the case of a structured competition, repeatability was underestimated (0.18 to 0.22). Our results show that the effect of an event, or category of event, is irrelevant in such a situation because ranks are independent of such an effect. The proposed mixture model pools horses according to their participation in different categories of competition during the period observed. This last model gave better results (repeatability 0.25), in particular, it provided an improved estimation of average values of competing ability of the horses in the different categories of events.

CONCLUSIONS

The underlying model was validated. A correct drawing of liabilities for the Gibbs sampler was provided. For a structured competition, the mixture model with a group effect assigned to horses gave the best results.

Genetic diversity of a large set of horse breeds raised in France assessed by microsatellite polymorphism

The genetic diversity and structure of horses raised in France were investigated using 11 microsatellite markers and 1679 animals belonging to 34 breeds. Between-breed differences explained about ten per cent of the total genetic diversity (Fst = 0.099). Values of expected heterozygosity ranged from 0.43 to 0.79 depending on the breed. According to genetic relationships, multivariate and structure analyses, breeds could be classified into four genetic differentiated groups: warm-blooded, draught, Nordic and pony breeds. Using complementary maximisation of diversity and aggregate diversity approaches, we conclude that particular efforts should be made to conserve five local breeds, namely the Boulonnais, Landais, Merens, Poitevin and Pottok breeds.

[Bacterial DNA alteration by plasma generated atomic nitrogen. Real-time PCR detection contribution]

AIMS

Bacterial sterilization by the technology of plasma in post-discharge shows a growing interest. The main appeal of this new process resides in its action at dry and low temperature (60 degrees C). This technology would be therefore useful for the complex medical equipment, sensitive to the oxidization, humidity and/or requiring a temperature lower than 60 degrees C. The objective of this survey is to demonstrate the activity of an atomic flux emanating a plasma of pure molecular nitrogen on the bacterial DNA: does the plasma of nitrogen damage the genetic material?

MATERIALS AND METHODS

The bacteria tested (Bacillus stearothermophilus, Staphylococcus aureus MRSA and Helicobacter pylori) are cultivated on suitable agar, and the bacterial DNA is extracted from every CFU by the technique of the columns (High Pure PCR Template, Roche). Every quantity of DNA extracted is diluted in 1 ml pure water. Then, 50 microl of each of these solutions of DNA are laid down in sterile Nunclon's plates holes, which undergo an advanced emptiness cycle during 60 minutes. The DNA residues will be then introduced during 40 minutes in a plasma sterilization surrounding wall (Plasmalyse), Satelec) where the debit of nitrogen, the pressure and the temperature are adjusted respectively to 1 L x min(-1), 5 Torrs and 60 degrees C. The DNA so ''plasmalysé'', is recovered then by aspiration in 500 mul pure water and processed to undergo an amplification/detection by Real-Time PCR (LightCycler2.0, Roche). The DNA ''plasmalysé'' will be compared to the intact DNA control(1), to the DNA control control(2) having undergone the cycle of emptiness solely, as well as to the DNA control(3) solely heated to 60 degrees C during 40 min.

RESULTS

The amplification curves demonstrated that the only advanced emptiness and the only heat don't have any activity on the bacterial DNA. On the other hand, the DNA ''plasmalysé'' shows a deterioration of the amplified sequences.

CONCLUSION

The genomic bacterial DNA, once extracted, is damaged by the gaseous flux of nitrogen plasma. A new sterilization process of the medical material will presumably impose itself in medium-term.