A genome scan for positive selection in thoroughbred horses

Deciphering climate resilience in Indian cattle breeds by selection signature analyses

The signature of selection is a crucial concept in evolutionary biology that refers to the pattern of genetic variation which arises in a population due to natural selection. In the context of climate adaptation, the signature of selection can reveal the genetic basis of adaptive traits that enable organisms to survive and thrive in changing environmental conditions. Breeds living in diverse agroecological zones exhibit genetic "footprints" within their genomes that mirror the influence of climate-induced selective pressures, subsequently impacting phenotypic variance. It is assumed that the genomes of animals residing in these regions have been altered through selection for various climatic adaptations. These regions are known as signatures of selection and can be identified using various summary statistics. We examined genotypic data from eight different cattle breeds (Gir, Hariana, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, and Tharparkar) that are adapted to diverse regional climates. To identify selection signature regions in this investigation, we used four intra-population statistics: Tajima's D, CLR, iHS, and ROH. In this study, we utilized Bovine 50 K chip data and four genome scan techniques to assess the genetic regions of positive selection for high-temperature adaptation. We have also performed a genome-wide investigation of genetic diversity, inbreeding, and effective population size in our target dataset. We identified potential regions for selection that are likely to be caused by adverse climatic conditions. We observed many adaptation genes in several potential selection signature areas. These include genes like HSPB2, HSPB3, HSP20, HSP90AB1, HSF4, HSPA1B, CLPB, GAP43, MITF, and MCHR1 which have been reported in the cattle populations that live in varied climatic regions. The findings demonstrated that genes involved in disease resistance and thermotolerance were subjected to intense selection. The findings have implications for marker-assisted breeding, understanding the genetic landscape of climate-induced adaptation, putting breeding and conservation programs into action.

Mathematical modeling suggests cytotoxic T lymphocytes control growth of B16 tumor cells in collagin-fibrin gels by cytolytic and non-lytic mechanisms

Cytotoxic T lymphocytes (CTLs) are important in controlling some viral infections, and therapies involving transfer of large numbers of cancer-specific CTLs have been successfully used to treat several types of cancers in humans. While molecular mechanisms of how CTLs kill their targets are relatively well understood we still lack solid quantitative understanding of the kinetics and efficiency at which CTLs kill their targets in different conditions. Collagen-fibrin gel-based assays provide a tissue-like environment for the migration of CTLs, making them an attractive system to study the cytotoxicity in vitro. Budhu et al. [1] systematically varied the number of peptide (SIINFEKL)- pulsed B16 melanoma cells and SIINFEKL-specific CTLs (OT-1) and measured remaining targets at different times after target and CTL co-inoculation into collagen-fibrin gels. The authors proposed that their data were consistent with a simple model in which tumors grow exponentially and are killed by CTLs at a per capita rate proportional to the CTL density in the gel. By fitting several alternative mathematical models to these data we found that this simple "exponential-growth-mass-action-killing" model does not precisely fit the data. However, determining the best fit model proved difficult because the best performing model was dependent on the specific dataset chosen for the analysis. When considering all data that include biologically realistic CTL concentrations ( E ≤ 10 ⁷ cell/ml) the model in which tumors grow exponentially and CTLs suppress tumor's growth non-lytically and kill tumors according to the mass-action law (SiGMA model) fitted the data with best quality. Results of power analysis suggested that longer experiments (∼ 3 - 4 days) with 4 measurements of B16 tumor cell concentrations for a range of CTL concentrations would best allow to discriminate between alternative models. Taken together, our results suggest that interactions between tumors and CTLs in collagen-fibrin gels are more complex than a simple exponential-growth- mass-action killing model and provide support for the hypothesis that CTLs impact on tumors may go beyond direct cytotoxicity.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Genetics of Thoroughbred Racehorse Performance

Thoroughbred horses have been selected for racing performance for more than 400 years. Despite continued selection, race times have not improved significantly during the past 60 years, raising the question of whether genetic variation for racing performance still exists. Studies using phenotypes such as race time, money earned, and handicapping, however, demonstrate that there is extensive variation within these traits and that they are heritable. Even so, these are poor measures of racing success since Thoroughbreds race at different ages and distances and on different types of tracks, and some may not race at all. With the advent of genomic tools, DNA variants are being identified that contribute to racing success. Aside from strong associations for myostatin variants with best racing distance, weak to modest associations with racing phenotypes are reported for other genomic regions. These data suggest that diverse genetic strategies have contributed to producing a successful racehorse, and genetic variation contributing to athleticism remains important. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Role of genes related to performance and reproduction of Thoroughbreds in training and breeding - A review

Thoroughbreds have been selected for speed and stamina since the 1700s. This selection resulted in structural and functional system-wide adaptations that enhanced physiological characteristics for outstanding speed of 61-71 kph (38-44 mph) between 1,000 and 3,200 m (5 furlongs - 2 miles). At present, horseracing is still an economically important industrial sector, therefore intensive research is underway to explore genes that allow the utilisation of genetic abilities and are significant in breeding and training. This study aims to provide an overview of genetic research and its applicability related to Thoroughbreds.

Polymorphisms at Myostatin Gene (MSTN) and the Associations with Sport Performances in Anglo-Arabian Racehorses

One hundred and eighty Anglo-Arabian horses running 1239 races were sampled for the present study. DNA was extracted from the blood and myostatin gene, MSTN, was genotyped. Moreover, prizes won and places were achieved for the 1239 races to perform association analyses between the different genotypes and sport traits. Two SNPs already reported in previous studies regarding the Thoroughbred breed, rs69472472 and rs397152648, were revealed as polymorphic. The linkage disequilibrium analysis investigating the haplotype structure of MSTN did not evidence any association block. Polymorphism at SNP rs397152648, previously known as g.66493737 T>C, significantly influenced sport traits, with heterozygous horses TC showing better results than homozygotes TT. The portion of variance due to the random effect of the individual animal, and the other phenotypic effects of sex, percentage of Arabian blood and race distance, computed together with the genotype at MSTN in the statistical models, exerted a significant influence. Hence, this information is useful to improve knowledge of the genetic profile of Anglo-Arabian horses and a possible selection for better sport performance.

Variability of ACOX1 Gene Polymorphisms across Different Horse Breeds with Regard to Selection Pressure

Simple Summary

The genetic mechanisms occurring in organisms are shaped by selection pressure. Features that ought to be useful under given conditions leave their marks on the genome in the form of mutations, thereby creating different alleles. In this study, five different horse breeds were examined to find the connection between an individual’s lifestyle and the presence of the peroxisomal acyl-coenzyme A oxidase 1 (ACOX1) gene, which is necessary for some metabolic pathways. Results indicated that different ACOX1 gene alleles play various roles in primitive breeds and domesticated horses. This led to the conclusion that the DNA profile can be rated on the basis of adaptation to living conditions, opening the gate for further investigation.

Abstract

The ACOX1 gene encodes peroxisomal acyl-coenzyme A oxidase 1, the first enzyme in the fatty acid β-oxidation pathway, which could be significant for organisms exposed to long periods of starvation and harsh living conditions. We hypothesized that variations within ACOX1, revealed by RNA Sequencing (RNA-Seq), might be based on adaptation to living conditions and had resulted from selection pressure. There were five different horse breeds used in this study, representing various utility types: Arabian, Thoroughbred, Polish Konik, draft horses, and Hucul. The single-nucleotide polymorphism (SNP) located in the ACOX1 (rs782885985) was used as a marker and was identified using the PCR restriction fragment length polymorphism method (PCR-RFLP). Results indicated extremely different genotype and allele distributions of the ACOX1 gene across breeds. A predominance of the G allele was exhibited in horses that had adapted to difficult environmental conditions, namely, Polish Konik and Huculs, which are considered to be primitive breeds. The prevalence of the T allele in Thoroughbreds indicated that ACOX1 is significant in energy metabolism during flat racing.

Analysis of Genetic Variability in the Argentine Polo Horse With a Panel of Microsatellite Markers

Argentine Polo (AP) is a young horse breed with a fast expansion because of an open registry policy and the application of assisted reproduction technologies. In the last years, AP showed a remarkable progress associated with the use of embryo transfer technology and intensive selection based on sport performance. However, these practices could have affected the genetic variability of the breed. To monitor these potential changes, a parentage panel of 11-15 microsatellites was investigated for changes in allele frequencies, heterozygosity, and exclusion probability over a 16 year period. Frequency of 36 alleles in 13 markers showed significant departures from expected changes because of random genetic drift. Six markers showed both allelic frequency increase and expected heterozygosity (He) reduction, suggesting the influence of selective breeding or hitchhiking effects. The average He decreased significantly although was lower than the observed heterozygosity, indicating a still low inbreeding level. Although the exclusion probability of the panel showed a trend to decrease, it is over 0.9995, the recommended value for equine parentage exclusion.

Genetic analysis of the endangered Cleveland Bay horse: A century of breeding characterised by pedigree and microsatellite data

The Cleveland Bay horse is one of the oldest equines in the United Kingdom, with pedigree data going back almost 300 years. The studbook is essentially closed and because of this, there are concerns about loss of genetic variation across generations. The breed is one of five equine breeds listed as “critical” (<300 registered adult breeding females) by the UK Rare Breeds Survival Trust in their annual Watchlist. Due to their critically endangered status, the current breadth of their genetic diversity is of concern, and assessment of this can lead to improved breed management strategies. Herein, both genealogical and molecular methods are combined in order to assess founder representation, lineage, and allelic diversity. Data from 15 microsatellite loci from a reference population of 402 individuals determined a loss of 91% and 48% of stallion and dam lines, respectively. Only 3 ancestors determine 50% of the genome in the living population, with 70% of maternal lineage being derived from 3 founder females, and all paternal lineages traced back to a single founder stallion. Methods and theory are described in detail in order to demonstrate the scope of this analysis for wider conservation strategies. We quantitatively demonstrate the critical nature of the genetic resources within the breed and offer a perspective on implementing this data in considered breed management strategies.

Variation in the SLC16A1 and the ACOX1 Genes Is Associated with Gallop Racing Performance in Arabian Horses

Arabian horses are not only one of the most ancient breeds in the world, but they are also one of the most appreciated racehorse breeds today. The breed generates attention for their phenomenal endurance ability and their capability for gallop racing. Consequently, genetic testing to select the best individuals is attracting ever increasing interests from the Arabian industry. As such, the aim of this study was to further investigate associations between performance and variation at candidate genes suspected of having a key role in Arabian gallop racing performance. Generalized linear models were fit to test associations between eight candidate gene variants and a variety of gallop racing performance traits in a sample of Arabian racehorses (n = 287). Two genes, solute carrier family 16 member 1 (SLC16A1) and acyl-CoA oxidase 1 (ACOX1), were significantly associated with multiple gallop racing performance traits, whereas another gene, actinin alpha 3 (ACTN3) was associated with best race distance. Previously established associations between these three genes and equine metabolism strongly suggest further investigation of these genes, and their relationship with Arabian horse performance is warranted.

A genome-wide scan for candidate lethal variants in Thoroughbred horses

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.

The effects of inbreeding on covering success, gestation length and foal sex ratio in Australian thoroughbred horses

Background

Horses produce only one foal from an eleven-month gestation period, making the maintenance of high reproductive rates essential. Genetic bottlenecks and inbreeding can increase the frequency of deleterious variants, resulting in reduced reproductive levels in a population. In this study we examined the influence of inbreeding levels on foaling rate, gestation length and secondary sex ratio in Australian Thoroughbred mares. We also investigated the genetic change in these traits throughout the history of the breed. Phenotypic data were obtained from 27,262 breeding records of Thoroughbred mares provided by three Australian stud farms. Inbreeding was estimated using the pedigree of each individual dating back to the foundation of the breed in the eighteenth century.

Results

While both gestation length and foaling rate were heritable, no measurable effect of inbreeding on either trait was found. However, we did find that the genetic value for both traits had decreased within recent generations. A number of environmental factors also had significant effects on foaling rate and gestation length. Secondary sex ratio had only an extremely small paternal heritable effect and was not susceptible to environmental influences.

Conclusions

In contrast to racing performance, inbreeding had no measurable effect on foaling rate or gestation length in Australian Thoroughbred horses. This could be because the level of inbreeding in the population examined is not high enough to show a discernible effect on reproductive traits. Populations that experience higher levels of inbreeding due to use of artificial reproductive technologies or extremely small population sizes may show a more pronounced reduction in natural foaling rate or gestation length. It is also possible that the intensive management techniques used in the Thoroughbred population masks any negative effects of inbreeding. The decrease in the genetic value of foaling rate is likely to be because horses with unfavourable genetic potential have not yet been selected out of the population. The change in genetic value of gestation length may be due to selective breeding favouring horses with shorter pregnancies. We also found that prioritising the mating of older mares, and avoiding out of season mating could lead to an increased breeding success.

Selection in Australian Thoroughbred horses acts on a locus associated with early two-year old speed

Thoroughbred horse racing is a global sport with major hubs in Europe, North America, Australasia and Japan. Regional preferences for certain traits have resulted in phenotypic variation that may result from adaptation to the local racing ecosystem. Here, we test the hypothesis that genes selected for regional phenotypic variation may be identified by analysis of selection signatures in pan-genomic SNP genotype data. Comparing Australian to non-Australian Thoroughbred horses (n = 99), the most highly differentiated loci in a composite selection signals (CSS) analysis were on ECA6 (34.75–34.85 Mb), ECA14 (33.2–33.52 Mb and 35.52–36.94 Mb) and ECA16 (24.28–26.52 Mb) in regions containing candidate genes for exercise adaptations including cardiac function (ARHGAP26, HBEGF, SRA1), synapse development and locomotion (APBB3, ATXN7, CLSTN3), stress response (NR3C1) and the skeletal muscle response to exercise (ARHGAP26, NDUFA2). In a genome-wide association study for field-measured speed in two-year-olds (n = 179) SNPs contained within the single association peak (33.2–35.6 Mb) overlapped with the ECA14 CSS signals and spanned a protocadherin gene cluster. Association tests using higher density SNP genotypes across the ECA14 locus identified a SNP within the PCDHGC5 gene associated with elite racing performance (n = 922). These results indicate that there may be differential selection for racing performance under racing and management conditions that are specific to certain geographic racing regions. In Australia breeders have principally selected horses for favourable genetic variants at loci containing genes that modulate behaviour, locomotion and skeletal muscle physiology that together appear to be contributing to early two-year-old speed.

Social Referencing in the Domestic Horse

Dogs and cats use human emotional information directed to an unfamiliar situation to guide their behavior, known as social referencing. It is not clear whether other domestic species show similar socio-cognitive abilities in interacting with humans. We investigated whether horses (n = 46) use human emotional information to adjust their behavior to a novel object and whether the behavior of horses differed depending on breed type. Horses were randomly assigned to one of two groups: an experimenter positioned in the middle of a test arena directed gaze and voice towards the novel object with either (a) a positive or (b) a negative emotional expression. The duration of subjects' position to the experimenter and the object in the arena, frequency of gazing behavior, and physical interactions (with either object or experimenter) were analyzed. Horses in the positive condition spent more time between the experimenter and object compared to horses in the negative condition, indicating less avoidance behavior towards the object. Horses in the negative condition gazed more often towards the object than horses in the positive condition, indicating increased vigilance behavior. Breed types differed in their behavior: thoroughbreds showed less human-directed behavior than warmbloods and ponies. Our results provide evidence that horses use emotional cues from humans to guide their behavior towards novel objects.

Genome-wide SNP analysis of Japanese Thoroughbred racehorses

The domestication process of plants and animals typically involves intense inbreeding and directional selection for various traits. Here, we genotyped 370 Japanese Thoroughbred horses using the recently developed 670k SNP array and performed various genome-wide analysis also using genotype data of other horse breeds. We identified a number of regions showing interesting patterns of polymorphisms. For instance, the region containing the MC1R locus associated with chestnut coat color may have been targeted by selection for a different mutation much earlier on than the recent selection for chestnut color. We also identified regions that show signatures of selection specific to Thoroughbreds. In addition, we found that intense inbreeding early in the history of the Thoroughbred breed and also before the formation of the breed has a significant impact on the genomic architecture of modern Thoroughbreds. Our study demonstrates that the horse 670k array can be utilized to gain important insight into the domestication process of horses and to understand the genetic basis of the phenotypic diversity in horses.

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).

Selection signatures in four German warmblood horse breeds: Tracing breeding history in the modern sport horse

The study of selection signatures helps to find genomic regions that have been under selective pressure and might host genes or variants that modulate important phenotypes. Such knowledge improves our understanding of how breeding programmes have shaped the genomes of livestock. In this study, 942 stallions were included from four, exemplarily chosen, German warmblood breeds with divergent historical and recent selection focus and different crossbreeding policies: Trakehner (N = 44), Holsteiner (N = 358), Hanoverian (N = 319) and Oldenburger (N = 221). Those breeds are nowadays bred for athletic performance and aptitude for show-jumping, dressage or eventing, with a particular focus of Holsteiner on the first discipline. Blood samples were collected during the health exams of the stallion preselections before licensing and were genotyped with the Illumina EquineSNP50 BeadChip. Autosomal markers were used for a multi-method search for signals of positive selection. Analyses within and across breeds were conducted by using the integrated Haplotype Score (iHS), cross-population Extended Haplotype Homozygosity (xpEHH) and Runs of Homozygosity (ROH). Oldenburger and Hanoverian showed very similar iHS signatures, but breed specificities were detected on multiple chromosomes with the xpEHH. The Trakehner clustered as a distinct group in a principal component analysis and also showed the highest number of ROHs, which reflects their historical bottleneck. Beside breed specific differences, we found shared selection signals in an across breed iHS analysis on chromosomes 1, 4 and 7. After investigation of these iHS signals and shared ROH for potential functional candidate genes and affected pathways including enrichment analyses, we suggest that genes affecting muscle functionality (TPM1, TMOD2-3, MYO5A, MYO5C), energy metabolism and growth (AEBP1, RALGAPA2, IGFBP1, IGFBP3-4), embryonic development (HOXB-complex) and fertility (THEGL, ZPBP1-2, TEX14, ZP1, SUN3 and CFAP61) have been targeted by selection in all breeds. Our findings also indicate selection pressure on KITLG, which is well-documented for influencing pigmentation.

Comparative population genomics unveils candidate genes for athletic performance in Hanoverians

Equine athletes have a genetic heritage that has been evolved for millions of years, which provides an opportunity to study the genetics of locomotion pattern and performance in mammals. The Hanoverian, a breed originating in Germany, is arguably among the most athletic of horse breeds, as well as possessing a balanced character and beautiful appearance. Here, we compared the whole genomes of Hanoverian with three other horse breeds (Akhal-Teke, Franches-Montagnes, and Standardbred), using the fixation index (Fst) and cross-population composite likelihood ratio (XP-CLR) methods for testing the multi-locus allele frequency differentiation between populations. We identified 299 and 485 positively selected genes using the Fst and XP-CLR methods, respectively. Further functional analyses showed that the ACTA1 gene is potentially involved in athletic performance in the Hanoverian breed, consistent with its role observed in human population. In addition, three other loci on chromosomes 1 and 20 were identified to be potentially involved in equine physical performance. The selected candidate genes identified in this study may be useful in current breeding efforts to develop improved breeds in regard to athletic performance.

Evaluation of recent changes in genetic variability in Japanese thoroughbred population based on a short tandem repeat parentage panel

The integrity of thoroughbreds is maintained under strict regulation involving DNA parentage testing, which is robust in a population with high genetic variability. The genetic variability of the thoroughbred population is possibly fluctuating because of selective breeding that has focused on adaptations for racing performance. To monitor genetic variability within the population and the effectiveness of short tandem repeat (STR) parentage testing, we investigated allele frequencies and the exclusion probability (PE) of 16-17 loci of a parentage panel in the Japanese thoroughbred population over 15 years. Expected heterozygosities (He) of 14 loci indicated a decreasing trend, and the average He of the population decreased significantly. Low genetic variability was possibly induced by a decrease in population size and a selective breeding bias. Four loci showed both a significant increase in allele frequency and a significant decrease in He; it is assumed that those loci were affected by positive selection for racing performance. There was a significant decrease in the PE because of the changes in genetic variability; however, it has remained over 0.99995. The current STR panel is still effective for parentage control, but it will be necessary to continuously monitor genetic variability, which has decreased over 15 years.

Identification and expression analysis of a novel miRNA derived from ERV-E1 LTR in Equus caballus

Horses (Equus caballus) have been domesticated and bred to enhance speed, strength, and agility. Members of the Equus caballus Endogenous Retrovirus (EqERV) family affect several of these abilities in horses. EqERV elements have been integrated in the horse genome during evolution and generate repeat elements such as long terminal repeats (LTRs). LTR sequences are involved in retrovirus replication and play an essential function in post-transcriptional control mechanisms, such as by providing binding sites for microRNAs (miRNAs) or generating miRNA precursors. In this study, we identified a novel miRNA derived from EqERV-E1 LTR using various bioinformatics tools. To examine the relationship between EqERV-E1 LTR and similar elements, we used BLAST2seq and phylogenetic analysis. LTR sequences were located in the untranslated region (UTR) of mRNAs and also formed the stem-loop secondary structure. The sequence was registered in the DDBJ database as LTR derived miRNA under the accession number corresponding to LC383797 (referred to eca-miR-1804). Quantitative polymerase chain reaction (qPCR) to confirm the expression of eca-miR-1804 and the similar miR-1255a, showed an almost identical expression pattern in eight different equine tissues. Therefore, these data imply that the LTR could function as an miRNA, which is expressed in the examined equine tissues. In addition, the current study provides inputs for additional functional studies concerning the LTR of other EqERV families.

A genome-wide association study for harness racing success in the Norwegian-Swedish coldblooded trotter reveals genes for learning and energy metabolism

Background

Although harness racing is of high economic importance to the global equine industry, significant genomic resources have yet to be applied to mapping harness racing success. To identify genomic regions associated with harness racing success, the current study performs genome-wide association analyses with three racing performance traits in the Norwegian-Swedish Coldblooded Trotter using the 670 K Axiom Equine Genotyping Array.

Results

Following quality control, 613 horses and 359,635 SNPs were retained for further analysis. After strict Bonferroni correction, nine genome-wide significant SNPs were identified for career earnings. No genome-wide significant SNPs were identified for number of gallops or best km time. However, four suggestive genome-wide significant SNPs were identified for number of gallops, while 19 were identified for best km time. Multiple genes related to intelligence, energy metabolism, and immune function were identified as potential candidate genes for harness racing success.

Conclusions

Apart from the physiological requirements needed for a harness racing horse to be successful, the results of the current study also advocate learning ability and memory as important elements for harness racing success. Further exploration into the mental capacity required for a horse to achieve racing success is likely warranted.

Electronic supplementary material

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

Integrated analysis of microRNA and mRNA expressions in peripheral blood leukocytes of Warmblood horses before and after exercise

Exercise capacity is a valuable trait in horses, and it has been used as a horse selection criterion. Although exercise affects molecular homeostasis and adaptation in horses, the mechanisms underlying these effects are not fully described. This study was carried out to identify changes in the blood profiles of microRNAs (miRNAs) and mRNAs induced by exercise in horse leukocytes. Total RNAs isolated from the peripheral blood leukocytes of four Warmblood horses before and after exercise were subjected to next-generation sequencing (NGS) and microarray analyses to determine the miRNA and mRNA expression profiles, respectively. The expressions of 6 miRNAs, including 4 known and 2 novel miRNAs, were altered by exercise. The predicted target genes of the differentially expressed miRNAs identified by NGS were matched to the exercise-induced mRNAs determined by microarray analysis. Five genes (LOC100050849, LOC100054517, KHDRBS3, LOC100053996, and LOC100062720) from the microarray analysis were matched to the predicted target genes of the 6 miRNAs. The subset of mRNAs and miRNAs affected by exercise in peripheral blood leukocytes may be useful in elucidating the molecular mechanisms of exercise-associated physiology in horses.

Identification on novel locus of dairy traits of Kazakh horse in Xinjiang

The utility of high-density single nucleotide polymorphism (SNP) data help to accurately identify genomic regions that have undergone positive selection. In this study, the Affymetrix Equine 670 K high-density SNP array was used to genotype Kazakh and Yili horse population. After quality control, 370,227 autosomal SNPs were used to detect selection signatures by using global fixation index (F_ST) and cross-population extended haplotype homozygosity (XP-EHH). The database of Ensemble, Genecards, and NCBI were used to make gene annotation and functional analysis. The results showed that there were 134 candidate SNPs overlapped between F_ST and XP-EHH in Kazakh horse. We also discovered some potential selective sweep regions associated with milk trait, including NUMB, LGALS2, ADCY8, SLC25A30, and CA8 genes. New findings from this research have potential value for milk traits selecting in horse.

A potential regulatory region near the EDN3 gene may control both harness racing performance and coat color variation in horses

The Swedish-Norwegian Coldblooded trotter and the heavier North-Swedish draught horse both descend from the North-Swedish horse, but the Coldblooded trotters have been selected for racing performance while the North-Swedish draught horse is mainly used for agricultural and forestry work. By comparing the genomes of Coldblooded trotters, North-Swedish draught horses and Standardbreds for a large number of single-nucleotide polymorphisms (SNPs), the aim of the study was to identify genetic regions that may be under selection for racing performance. We hypothesized that the selection for racing performance, in combination with unauthorized crossbreeding of Coldblooded trotters and Standardbreds, has created regions in the genome where the Coldblooded trotters and Standardbreds are similar, but differ from the North-Swedish draught horse. A fixation index (Fst) analysis was performed and sliding window Delta Fst values were calculated across the three breeds. Five windows, where the average Fst between Coldblooded trotters and Standardbreds was low and the average Fst between Coldblooded trotters and North-Swedish draught horses was high, were selected for further investigation. Associations between the most highly ranked SNPs and harness racing performance were analyzed in 400 raced Coldblooded trotters with race records. One SNP showed a significant association with racing performance, with the CC genotype appearing to be negatively associated. The SNP identified was genotyped in 1915 horses of 18 different breeds. The frequency of the TT genotype was high in breeds typically used for racing and show jumping while the frequency of the CC genotype was high in most pony breeds and draught horses. The closest gene in this region was the Endothelin3 gene (EDN3), a gene mainly involved in melanocyte and enteric neuron development. Both functional genetic and physiological studies are needed to fully understand the possible impacts of the gene on racing performance.

Founder-specific inbreeding depression affects racing performance in Thoroughbred horses

The Thoroughbred horse has played an important role in both sporting and economic aspects of society since the establishment of the breed in the 1700s. The extensive pedigree and phenotypic information available for the Thoroughbred horse population provides a unique opportunity to examine the effects of 300 years of selective breeding on genetic load. By analysing the relationship between inbreeding and racing performance of 135,572 individuals, we found that selective breeding has not efficiently alleviated the Australian Thoroughbred population of its genetic load. However, we found evidence for purging in the population that might have improved racing performance over time. Over 80% of inbreeding in the contemporary population is accounted for by a small number of ancestors from the foundation of the breed. Inbreeding to these ancestors has variable effects on fitness, demonstrating that an understanding of the distribution of genetic load is important in improving the phenotypic value of a population in the future. Our findings hold value not only for Thoroughbred and other domestic breeds, but also for small and endangered populations where such comprehensive information is not available.

Analysis of cross-population differentiation between Thoroughbred and Jeju horses

Objective

This study was intended to identify genes positively selected in Thoroughbred horses (THBs) that potentially contribute to their running performances.

Methods

The genomes of THB and Jeju horses (JH, Korean native horse) were compared to identify genes positively selected in THB. We performed cross-population extended haplotype homozygosity (XP-EHH) and cross-population composite likelihood ratio test (XP-CLR) statistical methods for our analysis using whole genome resequencing data of 14 THB and 6 JH.

Results

We identified 98 (XP-EHH) and 200 (XP-CLR) genes that are under positive selection in THB. Gene enrichment analysis identified 72 gene ontology biological process (GO BP) terms. The genes and GO BP terms explained some of THB’s characteristics such as immunity, energy metabolism and eye size and function related to running performances. GO BP terms that play key roles in several cell signaling mechanisms, which affected ocular size and visual functions were identified. GO BP term Eye photoreceptor cell differentiation is among the terms annotated presumed to affect eye size.

Conclusion

Our analysis revealed some positively selected candidate genes in THB related to their racing performances. The genes detected are related to the immunity, ocular size and function, and energy metabolism.

Identification and Expression Analyses of Equine Endogenous Retroviruses in Horses

Endogenous retroviruses (ERVs) have been integrated into vertebrate genomes and have momentously affected host organisms. Horses (Equus caballus) have been domesticated and selected for elite racing ability over centuries. ERVs played an important role in the evolutionary diversification of the horse genome. In the present study, we identified six equine ERV families (EqERVs-E1, I1, M2, P1, S1, and Y4), their full-length viral open reading frames (ORFs), and elucidated their phylogenetic relationships. The divergence time of EqERV families assuming an evolutionary rate of 0.2%/Myr indicated that EqERV-S3 (75.4 million years ago; mya) on chromosome 10 is an old EqERV family and EqERV-P5 (1.2 Mya) on chromosome 12 is a young member. During the evolutionary diversification of horses, the EqERV-I family diverged 1.7 Mya to 38.7 Mya. Reverse transcription quantitative real-time PCR (RT-qPCR) amplification of EqERV pol genes showed greater expression in the cerebellum of the Jeju horse than the Thoroughbred horse. These results could contribute further dynamic studies for horse genome in relation to EqERV gene function.

Foundations of performance – factors that contribute to excellence in equine exercise

Horses are renowned for their incredible capacity for a range of athletic activities, and participation in athletic events arguably represents the most critical strut of the equine industry. Successful performance is typically a primary focus during participation in competitive athletic events, and relies upon a variety of innate physiological and structural factors of the athlete. However, a wide range of external factors also influence performance, and many of these can be readily manipulated. Therefore, thorough assessment of the individual’s inherent capacity for a specific athletic discipline must be combined with optimisation of external factors including nutrition and training to promote excellent performance. Recent progress in methods of athlete selection and monitoring of training responses are assisting continued improvements in equine performance.

Molecular analysis of alternative transcripts of equine AXL receptor tyrosine kinase gene

OBJECTIVE

Since athletic performance is a most importance trait in horses, most research focused on physiological and physical studies of horse athletic abilities. In contrast, the molecular analysis as well as the regulatory pathway studies remain insufficient for evaluation and prediction of horse athletic abilities. In our previous study, we identified AXL receptor tyrosine kinase (AXL) gene which was expressed as alternative spliced isoforms in skeletal muscle during exercise. In the present study, we validated two AXL alternative splicing transcripts (named as AXLa for long form and AXLb for short form) in equine skeletal muscle to gain insight(s) into the role of each alternative transcript during exercise.

METHODS

We validated two isoforms of AXL transcripts in horse tissues by reverse transcriptase polymerase chain reaction (RT-PCR), and then cloned the transcripts to confirm the alternative locus and its sequences. Additionally, we examined the expression patterns of AXLa and AXLb transcripts in horse tissues by quantitative RT-PCR (qRT-PCR).

RESULTS

Both of AXLa and AXLb transcripts were expressed in horse skeletal muscle and the expression levels were significantly increased after exercise. The sequencing analysis showed that there was an alternative splicing event at exon 11 between AXLa and AXLb transcripts. 3-dimentional (3D) prediction of the alternative protein structures revealed that the structural distance of the connective region between fibronectin type 3 (FN3) and immunoglobin (Ig) domain was different between two alternative isoforms.

CONCLUSION

It is assumed that the expression patterns of AXLa and AXLb transcripts would be involved in regulation of exercise-induced stress in horse muscle possibly through an NF-κB signaling pathway. Further study is necessary to uncover biological function(s) and significance of the alternative splicing isoforms in race horse skeletal muscle.

Association of low race performance with mtDNA haplogroup L3b of Australian thoroughbred horses

Mitochondrial DNA (mtDNA) encodes the genes for respiratory chain sub-units that determine the efficiency of oxidative phosphorylation in mitochondria. The aim of this study was to determine if there were any haplogroups and variants in mtDNA that could be associated with athletic performance of Thoroughbred horses. The whole mitochondrial genomes of 53 maternally unrelated Australian Thoroughbred horses were sequenced and an association study was performed with the competition histories of 1123 horses within their maternal lineages. A horse mtDNA phylogenetic tree was constructed based on a total of 195 sequences (including 142 from previous reports). The association analysis showed that the sample groups with poor racing performance history were enriched in haplogroup L3b (p = .0003) and its sub-haplogroup L3b1a (p = .0007), while those that had elite performance appeared to be not significantly associated with haplogroups G2 and L3a1a1a (p > .05). Haplogroup L3b and L3b1a bear two and five specific variants of which variant T1458C (site 345 in 16s rRNA) is the only potential functional variant. Furthermore, secondary reconstruction of 16s RNA showed considerable differences between two types of 16s RNA molecules (with and without T1458C), indicating a potential functional effect. The results suggested that haplogroup L3b, could have a negative association with elite performance. The T1458C mutation harboured in haplogroup L3b could have a functional effect that is related to poor athletic performance.

Transcriptomic evaluation of bovine blastocysts obtained from peri-pubertal oocyte donors

Assisted reproduction technologies (ART) and high selection pressure in the dairy industry are leading towards the use of younger females for reproduction, thereby reducing the interval between generations. This situation may have a negative impact on embryo quality, thus reducing the success rate of the procedures. This study aimed to document the effects of oocyte donor age on embryo quality, at the transcriptomic level, in order to characterize the effects of using young females for reproduction purpose. Young Holstein heifers (n = 10) were used at three different ages for ovarian stimulation protocols and oocyte collections (at 8, 11 and 14 months). All of the oocytes were fertilized in vitro with the semen of one adult bull, generating three lots of embryos per animal. Each animal was its own control for the evaluation of the effects of age. The EmbryoGENE platform was used for the assessment of gene expression patterns at the blastocyst stage. Embryos from animals at 8 vs 14 months and at 11 vs 14 months were used for microarray hybridization. Validation was done by performing RT-qPCR on seven candidate genes. Age-related contrast analysis (8 vs 14 mo and 11 vs 14 mo) identified 242 differentially expressed genes (DEGs) for the first contrast, and 296 for the second. The analysis of the molecular and biological functions of the DEGs suggests a metabolic cause to explain the differences that are observed between embryos from immature and adult subjects. The mTOR and PPAR signaling pathways, as well as the NRF2-mediated oxidative stress response pathways were among the gene expression pathways affected by donor age. In conclusion, the main differences between embryos produced at peri-pubertal ages are related to metabolic conditions resulting in a higher impact of in vitro conditions on blastocyts from younger heifers.

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH ) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED ). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome-wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.

Diversity and Evolutionary Analysis of Iron-Containing (Type-III) Alcohol Dehydrogenases in Eukaryotes

BACKGROUND

Alcohol dehydrogenase (ADH) activity is widely distributed in the three domains of life. Currently, there are three non-homologous NAD(P)+-dependent ADH families reported: Type I ADH comprises Zn-dependent ADHs; type II ADH comprises short-chain ADHs described first in Drosophila; and, type III ADH comprises iron-containing ADHs (FeADHs). These three families arose independently throughout evolution and possess different structures and mechanisms of reaction. While types I and II ADHs have been extensively studied, analyses about the evolution and diversity of (type III) FeADHs have not been published yet. Therefore in this work, a phylogenetic analysis of FeADHs was performed to get insights into the evolution of this protein family, as well as explore the diversity of FeADHs in eukaryotes.

PRINCIPAL FINDINGS

Results showed that FeADHs from eukaryotes are distributed in thirteen protein subfamilies, eight of them possessing protein sequences distributed in the three domains of life. Interestingly, none of these protein subfamilies possess protein sequences found simultaneously in animals, plants and fungi. Many FeADHs are activated by or contain Fe2+, but many others bind to a variety of metals, or even lack of metal cofactor. Animal FeADHs are found in just one protein subfamily, the hydroxyacid-oxoacid transhydrogenase (HOT) subfamily, which includes protein sequences widely distributed in fungi, but not in plants), and in several taxa from lower eukaryotes, bacteria and archaea. Fungi FeADHs are found mainly in two subfamilies: HOT and maleylacetate reductase (MAR), but some can be found also in other three different protein subfamilies. Plant FeADHs are found only in chlorophyta but not in higher plants, and are distributed in three different protein subfamilies.

CONCLUSIONS/SIGNIFICANCE

FeADHs are a diverse and ancient protein family that shares a common 3D scaffold with a patchy distribution in eukaryotes. The majority of sequenced FeADHs from eukaryotes are distributed in just two subfamilies, HOT and MAR (found mainly in animals and fungi). These two subfamilies comprise almost 85% of all sequenced FeADHs in eukaryotes.

Glucose-independent Acetate Metabolism Promotes Melanoma Cell Survival and Tumor Growth

Tumors rely on multiple nutrients to meet cellular bioenergetics and macromolecular synthesis demands of rapidly dividing cells. Although the role of glucose and glutamine in cancer metabolism is well understood, the relative contribution of acetate metabolism remains to be clarified. We show that glutamine supplementation is not sufficient to prevent loss of cell viability in a subset of glucose-deprived melanoma cells, but synergizes with acetate to support cell survival. Glucose-deprived melanoma cells depend on both oxidative phosphorylation and acetate metabolism for cell survival. Acetate supplementation significantly contributed to maintenance of ATP levels in glucose-starved cells. Unlike acetate, short chain fatty acids such as butyrate and propionate failed to prevent loss of cell viability from glucose deprivation. In vivo studies revealed that in addition to nucleo-cytoplasmic acetate assimilating enzyme ACSS2, mitochondrial ACSS1 was critical for melanoma tumor growth in mice. Our data indicate that acetate metabolism may be a potential therapeutic target for BRAF mutant melanoma.

Thermal manipulation of the chicken embryo triggers differential gene expression in response to a later heat challenge

Background

Meat type chickens have limited capacities to cope with high environmental temperatures, this sometimes leading to mortality on farms and subsequent economic losses. A strategy to alleviate this problem is to enhance adaptive capacities to face heat exposure using thermal manipulation (TM) during embryogenesis. This strategy was shown to improve thermotolerance during their life span. The aim of this study was to determine the effects of TM (39.5 °C, 12 h/24 vs 37.8 °C from d7 to d16 of embryogenesis) and of a subsequent heat challenge (32 °C for 5 h) applied on d34 on gene expression in the Pectoralis major muscle (PM). A chicken gene expression microarray (8 × 60 K) was used to compare muscle gene expression profiles of Control (C characterized by relatively high body temperatures, Tb) and TM chickens (characterized by a relatively low Tb) reared at 21 °C and at 32 °C (CHC and TMHC, respectively) in a dye-swap design with four comparisons and 8 broilers per treatment. Real-time quantitative PCR (RT-qPCR) was subsequently performed to validate differential expression in each comparison. Gene ontology, clustering and network building strategies were then used to identify pathways affected by TM and heat challenge.

Results

Among the genes differentially expressed (DE) in the PM (1.5 % of total probes), 28 were found to be differentially expressed between C and TM, 128 between CHC and C, and 759 between TMHC and TM. No DE gene was found between TMHC and CHC broilers. The majority of DE genes analyzed by RT-qPCR were validated. In the TM/C comparison, DE genes were involved in energy metabolism and mitochondrial function, cell proliferation, vascularization and muscle growth; when comparing heat-exposed chickens to their own controls, TM broilers developed more specific pathways than C, especially involving genes related to metabolism, stress response, vascularization, anti-apoptotic and epigenetic processes.

Conclusions

This study improved the understanding of the long-term effects of TM on PM muscle. TM broilers displaying low Tb may have lower metabolic intensity in the muscle, resulting in decreased metabolic heat production, whereas modifications in vascularization may enhance heat loss. These specific changes could in part explain the better adaptation of TM broilers to heat.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2661-y) contains supplementary material, which is available to authorized users.

Transcriptomic profile adaptations following exposure of equine satellite cells to nutriactive phytochemical gamma-oryzanol

BACKGROUND

Adult skeletal muscle myogenesis depends on the activation of satellite cells that have the potential to differentiate into new fibers. Gamma-oryzanol (GO), a commercially available nutriactive phytochemical, has gained global interest on account of its muscle-building and regenerating effects. Here, we investigated GO for its potential influence on myogenesis, using equine satellite cell culture model, since the horse is a unique animal, bred and exercised for competitive sport. To our knowledge, this is the first report where the global gene expression in cultured equine satellite cells has been described.

METHODS

Equine satellite cells were isolated from semitendinosus muscle and cultured until the second day of differentiation. Differentiating cells were incubated with GO for the next 24 h. Subsequently, total RNA from GO-treated and control cells was isolated, amplified, labeled, and hybridized to two-color Horse Gene Expression Microarray slides. Quantitative PCR was used for the validation of microarray data.

RESULTS

Our results revealed 58 genes with changed expression in GO-treated vs. control cells. Analysis of expression changes suggests that various processes are reinforced by GO in differentiating equine satellite cells, including inhibition of myoblast differentiation, increased proliferation and differentiation, stress response, and increased myogenic lineage commitment.

CONCLUSIONS

The present study may confirm putative muscle-enhancing abilities of GO; however, the collective role of GO in skeletal myogenesis remains equivocal. The diversity of these changes is likely due to heterogenous growth rate of cells in primary culture. Genes identified in our study, modulated by the presence of GO, may become potential targets of future research investigating impact of this supplement in skeletal muscle on proteomic and biochemical level.

Interleukin-22-Induced Antimicrobial Phospholipase A2 Group IIA Mediates Protective Innate Immunity of Nonhematopoietic Cells against Listeria monocytogenes

Listeria monocytogenes is a bacterial pathogen which establishes intracellular parasitism in various cells, including macrophages and nonhematopoietic cells, such as hepatocytes. It has been reported that several proinflammatory cytokines have pivotal roles in innate protection against L. monocytogenes infection. We found that a proinflammatory cytokine, interleukin 22 (IL-22), was expressed by CD3(+) CD4(+) T cells at an early stage of L. monocytogenes infection in mice. To assess the influence of IL-22 on L. monocytogenes infection in hepatocytes, cells of a human hepatocellular carcinoma line, HepG2, were treated with IL-22 before L. monocytogenes infection in vitro. Gene expression analysis of the IL-22-treated HepG2 cells identified phospholipase A2 group IIA (PLA2G2A) as an upregulated antimicrobial molecule. Addition of recombinant PLA2G2A to the HepG2 culture significantly suppressed L. monocytogenes infection. Culture supernatant of the IL-22-treated HepG2 cells contained bactericidal activity against L. monocytogenes, and the activity was abrogated by a specific PLA2G2A inhibitor, demonstrating that HepG2 cells secreted PLA2G2A, which killed extracellular L. monocytogenes. Furthermore, colocalization of PLA2G2A and L. monocytogenes was detected in the IL-22-treated infected HepG2 cells, which suggests involvement of PLA2G2A in the mechanism of intracellular killing of L. monocytogenes by HepG2 cells. These results suggest that IL-22 induced at an early stage of L. monocytogenes infection enhances innate immunity against L. monocytogenes in the liver by stimulating hepatocytes to produce an antimicrobial molecule, PLA2G2A.

Advanced Running Performance by Genetic Predisposition in Male Dummerstorf Marathon Mice (DUhTP) Reveals Higher Sterol Regulatory Element-Binding Protein (SREBP) Related mRNA Expression in the Liver and Higher Serum Levels of Progesterone

Long-term-selected DUhTP mice represent a non-inbred model for inborn physical high-performance without previous training. Abundance of hepatic mRNA in 70-day male DUhTP and control mice was analyzed using the Affymetrix mouse array 430A 2.0. Differential expression analysis with PLIER corrected data was performed using AltAnalyze. Searching for over-representation in biochemical pathways revealed cholesterol metabolism being most prominently affected in DUhTP compared to unselected control mice. Furthermore, pathway analysis by AltAnalyze plus PathVisio indicated significant induction of glycolysis, fatty acid synthesis and cholesterol biosynthesis in the liver of DUhTP mice versus unselected control mice. In contrast, gluconeogenesis was partially inactivated as judged from the analysis of hepatic mRNA transcript abundance in DUhTP mice. Analysis of mRNA transcripts related to steroid hormone metabolism inferred elevated synthesis of progesterone and reduced levels of sex steroids. Abundance of steroid delta isomerase-5 mRNA (Hsd3b5, FC 4.97) was increased and steroid 17-alpha-monooxygenase mRNA (Cyp17a1, FC -11.6) was massively diminished in the liver of DUhTP mice. Assessment of steroid profiles by LC-MS revealed increased levels of progesterone and decreased levels of sex steroids in serum from DUhTP mice versus controls. Analysis of hepatic mRNA transcript abundance indicates that sterol regulatory element-binding protein-1 (SREBP-1) may play a major role in metabolic pathway activation in the marathon mouse model DUhTP. Thus, results from bioinformatics modeling of hepatic mRNA transcript abundance correlated with direct steroid analysis by mass spectrometry and further indicated functions of SREBP-1 and steroid hormones for endurance performance in DUhTP mice.

Skeletal muscle adaptations and muscle genomics of performance horses

Skeletal muscles in horses are characterised by specific adaptations, which are the result of the natural evolution of the horse as a grazing animal, centuries of selective breeding and the adaptability of this tissue in response to training. These adaptations include an increased muscle mass relative to body weight, a great locomotor efficiency based upon an admirable muscle-tendon architectural design and an adaptable fibre-type composition with intrinsic shortening velocities greater than would be predicted from an animal of comparable body size. Furthermore, equine skeletal muscles have a high mitochondrial volume that permits a higher whole animal aerobic capacity, as well as large intramuscular stores of energy substrates (glycogen in particular). Finally, high buffer and lactate transport capacities preserve muscles against fatigue during anaerobic exercise. Many of these adaptations can improve with training. The publication of the equine genome sequence in 2009 has provided a major advance towards an improved understanding of equine muscle physiology. Equine muscle genomics studies have revealed a number of genes associated with elite physical performance and have also identified changes in structural and metabolic genes following exercise and training. Genes involved in muscle growth, muscle contraction and specific metabolic pathways have been found to be functionally relevant for the early performance evaluation of elite athletic horses. The candidate genes discussed in this review are important for a healthy individual to improve performance. However, muscle performance limiting conditions are widespread in horses and many of these conditions are also genetically influenced.

A Genome-wide Scan for Selective Sweeps in Racing Horses

Using next-generation sequencing, we conducted a genome-wide scan of selective sweeps associated with selection toward genetic improvement in Thoroughbreds. We investigated potential phenotypic consequence of putative candidate loci by candidate gene association mapping for the finishing time in 240 Thoroughbred horses. We found a significant association with the trait for Ral GApase alpha 2 (RALGAP2) that regulates a variety of cellular processes of signal trafficking. Neighboring genes around RALGAP2 included insulinoma-associated 1 (INSM1), pallid (PLDN), and Ras and Rab interactor 2 (RIN2) genes have similar roles in signal trafficking, suggesting that a co-evolving gene cluster located on the chromosome 22 is under strong artificial selection in racehorses.

Runs of homozygosity reveal signatures of positive selection for reproduction traits in breed and non-breed horses

Background

Modern horses represent heterogeneous populations specifically selected for appearance and performance. Genomic regions under high selective pressure show characteristic runs of homozygosity (ROH) which represent a low genetic diversity. This study aims at detecting the number and functional distribution of ROHs in different horse populations using next generation sequencing data.

Methods

Next generation sequencing was performed for two Sorraia, one Dülmen Horse, one Arabian, one Saxon-Thuringian Heavy Warmblood, one Thoroughbred and four Hanoverian. After quality control reads were mapped to the reference genome EquCab2.70. ROH detection was performed using PLINK, version 1.07 for a trimmed dataset with 11,325,777 SNPs and a mean read depth of 12. Stretches with homozygous genotypes of >40 kb as well as >400 kb were defined as ROHs. SNPs within consensus ROHs were tested for neutrality. Functional classification was done for genes annotated within ROHs using PANTHER gene list analysis and functional variants were tested for their distribution among breed or non-breed groups.

Results

ROH detection was performed using whole genome sequences of ten horses of six populations representing various breed types and non-breed horses. In total, an average number of 3492 ROHs were detected in windows of a minimum of 50 consecutive homozygous SNPs and an average number of 292 ROHs in windows of 500 consecutive homozygous SNPs. Functional analyses of private ROHs in each horse revealed a high frequency of genes affecting cellular, metabolic, developmental, immune system and reproduction processes. In non-breed horses, 198 ROHs in 50-SNP windows and seven ROHs in 500-SNP windows showed an enrichment of genes involved in reproduction, embryonic development, energy metabolism, muscle and cardiac development whereas all seven breed horses revealed only three common ROHs in 50-SNP windows harboring the fertility-related gene YES1. In the Hanoverian, a total of 18 private ROHs could be shown to be located in the region of genes potentially involved in neurologic control, signaling, glycogen balance and reproduction. Comparative analysis of homozygous stretches common in all ten horses displayed three ROHs which were all located in the region of KITLG, the ligand of KIT known to be involved in melanogenesis, haematopoiesis and gametogenesis.

Conclusions

The results of this study give a comprehensive insight into the frequency and number of ROHs in various horses and their potential influence on population diversity and selection pressures. Comparisons of breed and non-breed horses suggest a significant artificial as well as natural selection pressure on reproduction performance in all types of horse populations.

Electronic supplementary material

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

Fifteen years of genomewide scans for selection: trends, lessons and unaddressed genetic sources of complication

Genomewide scans for natural selection (GWSS) have become increasingly common over the last 15 years due to increased availability of genome-scale genetic data. Here, we report a representative survey of GWSS from 1999 to present and find that (i) between 1999 and 2009, 35 of 49 (71%) GWSS focused on human, while from 2010 to present, only 38 of 83 (46%) of GWSS focused on human, indicating increased focus on nonmodel organisms; (ii) the large majority of GWSS incorporate interpopulation or interspecific comparisons using, for example F(ST), cross-population extended haplotype homozygosity or the ratio of nonsynonymous to synonymous substitutions; (iii) most GWSS focus on detection of directional selection rather than other modes such as balancing selection; and (iv) in human GWSS, there is a clear shift after 2004 from microsatellite markers to dense SNP data. A survey of GWSS meant to identify loci positively selected in response to severe hypoxic conditions support an approach to GWSS in which a list of a priori candidate genes based on potential selective pressures are used to filter the list of significant hits a posteriori. We also discuss four frequently ignored determinants of genomic heterogeneity that complicate GWSS: mutation, recombination, selection and the genetic architecture of adaptive traits. We recommend that GWSS methodology should better incorporate aspects of genomewide heterogeneity using empirical estimates of relevant parameters and/or realistic, whole-chromosome simulations to improve interpretation of GWSS results. Finally, we argue that knowledge of potential selective agents improves interpretation of GWSS results and that new methods focused on correlations between environmental variables and genetic variation can help automate this approach.

Genome-wide selection signatures in Pinzgau cattle

The aim of this study was to identify the evidence of recent selection based on estimation of the integrated Haplotype Score (iHS), population differentiation index (FST) and characterize affected regions near QTL associated with traits under strong selection in Pinzgau cattle. In total 21 Austrian and 19 Slovak purebreed bulls genotyped with Illumina bovineHD and  bovineSNP50 BeadChip were used to identify genomic regions under selection. Only autosomal loci with call rate higher than 90%, minor allele frequency higher than 0.01 and Hardy-Weinberg equlibrium limit of 0.001 were included in the subsequent analyses of selection sweeps presence. The final dataset was consisted from 30538 SNPs with 81.86 kb average adjacent SNPs spacing. The iHS score were averaged into non-overlapping 500 kb segments across the genome. The FST values were also plotted against genome position based on sliding windows approach and averaged over 8 consecutive SNPs. Based on integrated Haplotype Score evaluation only 7 regions with iHS score higher than 1.7 was found. The average iHS score observed for each adjacent syntenic regions indicated slight effect of recent selection in analysed group of Pinzgau bulls. The level of genetic differentiation between Austrian and Slovak bulls estimated based on FST index was low. Only 24% of FST values calculated for each SNP was greather than 0.01. By using sliding windows approach was found that 5% of analysed windows had higher value than 0.01. Our results indicated use of similar selection scheme in breeding programs of Slovak and Austrian Pinzgau bulls. The evidence for genome-wide association between signatures of selection and regions affecting complex traits such as milk production was insignificant, because the loci in segments identified as affected by selection were very distant from each other. Identification of genomic regions that may be under pressure of selection for phenotypic traits to better understanding of the relationship between genotype and phenotype is one of the challenges for livestock genetics.

Detection of selection signatures in dairy and beef cattle using high-density genomic information

Background

Artificial selection for economically important traits in cattle is expected to have left distinctive selection signatures on the genome. Access to high-density genotypes facilitates the accurate identification of genomic regions that have undergone positive selection. These findings help to better elucidate the mechanisms of selection and to identify candidate genes of interest to breeding programs.

Results

Information on 705 243 autosomal single nucleotide polymorphisms (SNPs) in 3122 dairy and beef male animals from seven cattle breeds (Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental) were used to detect selection signatures by applying two complementary methods, integrated haplotype score (iHS) and global fixation index (F_ST). To control for false positive results, we used false discovery rate (FDR) adjustment to calculate adjusted iHS within each breed and the genome-wide significance level was about 0.003. Using the iHS method, 83, 92, 91, 101, 85, 101 and 86 significant genomic regions were detected for Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental cattle, respectively. None of these regions was common to all seven breeds. Using the F_ST approach, 704 individual SNPs were detected across breeds. Annotation of the regions of the genome that showed selection signatures revealed several interesting candidate genes i.e. DGAT1, ABCG2, MSTN, CAPN3, FABP3, CHCHD7, PLAG1, JAZF1, PRKG2, ACTC1, TBC1D1, GHR, BMP2, TSG1, LYN, KIT and MC1R that play a role in milk production, reproduction, body size, muscle formation or coat color. Fifty-seven common candidate genes were found by both the iHS and global F_ST methods across the seven breeds. Moreover, many novel genomic regions and genes were detected within the regions that showed selection signatures; for some candidate genes, signatures of positive selection exist in the human genome. Multilevel bioinformatic analyses of the detected candidate genes suggested that the PPAR pathway may have been subjected to positive selection.

Conclusions

This study provides a high-resolution bovine genomic map of positive selection signatures that are either specific to one breed or common to a subset of the seven breeds analyzed. Our results will contribute to the detection of functional candidate genes that have undergone positive selection in future studies.

Electronic supplementary material

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