Archaeological and genetic evidence concerning the time and mode of wild horse (Equus ferus) domestication is still debated. High levels of genetic diversity in horse mtDNA have been detected when analyzing the control region; recurrent mutations, however, tend to blur the structure of the phylogenetic tree. Here, we brought the horse mtDNA phylogeny to the highest level of molecular resolution by analyzing 83 mitochondrial genomes from modern horses across Asia, Europe, the Middle East, and the Americas. Our data reveal 18 major haplogroups (A–R) with radiation times that are mostly confined to the Neolithic and later periods and place the root of the phylogeny corresponding to the Ancestral Mare Mitogenome at ∼130–160 thousand years ago. All haplogroups were detected in modern horses from Asia, but F was only found inE. przewalskii—the only remaining wild horse. Therefore, a wide range of matrilineal lineages from the extinctE. ferusunderwent domestication in the Eurasian steppes during the Eneolithic period and were transmitted to modernE. caballusbreeds. Importantly, now that the major horse haplogroups have been defined, each with diagnostic mutational motifs (in both the coding and control regions), these haplotypes could be easily used to (i) classify well-preserved ancient remains, (ii) (re)assess the haplogroup variation of modern breeds, including Thoroughbreds, and (iii) evaluate the possible role of mtDNA backgrounds in racehorse performance.

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F_ST-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse. A breed of the horse typically consists of individuals sharing very similar aesthetic and performance traits. However, a great deal of variation in traits exists between breeds. The range of variation observed among breeds can be illustrated by the size difference between the Miniature horse (0.74 m and 100 kg) and draft horse (1.8 m and 900 kg), or by comparing the optimum racing distance of the Quarter Horse (1/4 mile) to that of the Arabian (100 miles or more). In this study, we exploited the breed structure of the horse to identify regions of the genome that are significantly different between breeds and therefore may harbor genes and genetic variants targeted by selective breeding. This work resulted in the identification of variants in the Paint and Quarter Horse significantly associated with altered muscle fiber type proportions favorable for increased sprinting ability. A strong signature of selection was also identified in breeds that perform alternative gaits, and several genomic regions identified are hypothesized to be involved in the determination of size. This study has demonstrated the utility of this approach for studying the equine genome and is the first to show a functional consequence of selective breeding in the horse.

Adequate stress response is a critical factor during athlete horses' training and is central to our capacity to obtain better performances while safeguarding animal welfare.

Bovine milk is important for human nutrition, but its fat content is often criticized as a risk factor in cardiovascular disease. Selective breeding programs could be used to alter the fatty acid (FA) composition of bovine milk to improve the healthiness of dairy products for human consumption. Here, we performed a genome-wide association study (GWAS) on bovine milk to identify genomic regions or specific genes associated with FA profile and to investigate genetic differences between the Italian Simmental (IS) and Italian Holstein (IH) breeds. To achieve this, we first characterized milk samples from 416 IS cows and 436 IH cows for their fat profile by gas chromatography. Subjects were genotyped with single nucleotide polymorphism array and a single-marker regression model for GWAS was performed. Our findings confirm previously reported quantitative trait loci strongly associated with bovine milk fat composition. More specifically, our GWAS results revealed significant signals on chromosomes Bos taurus autosome 19 and 26 for milk FA. Further analysis using a gene-centric approach and pathway meta-analysis identified not only some well-known genes underlying quantitative trait loci for milk FA components, such as FASN, SCD, and DGAT1, but also other significant candidate genes, including some with functional roles in pathways related to "Lipid metabolism." Highlighted genes related to FA profile include ECI2, PCYT2, DCXR, G6PC3, PYCR1, and ALG12 in IS, and CYP17A1, ACO2, PI4K2A, GOT1, GPT, NT5C2, PDE6G, POLR3H, and COX15 in IH. Overall, the breed-specific association outcomes reflect differences in the genetic backgrounds of the IS and IH breeds and their selective breeding histories.

Milk extracellular vesicles (mEVs) seem to be one of the main maternal messages delivery systems. Extracellular vesicles (EVs) are micro/nano-sized membrane-bound structures enclosing signaling molecules and thus acting as signal mediators between distant cells and/or tissues, exerting biological effects such as immune modulation and pro-regenerative activity. Milk is also a unique, scalable, and reliable source of EVs. Our aim was to characterize the RNA content of cow, donkey, and goat mEVs through transcriptomic analysis of mRNA and small RNA libraries. Over 10,000 transcripts and 2000 small RNAs were expressed in mEVs of each species. Among the most represented transcripts, 110 mRNAs were common between the species with cow acting as the most divergent. The most represented small RNA class was miRNA in all the species, with 10 shared miRNAs having high impact on the immune regulatory function. Functional analysis for the most abundant mRNAs shows epigenetic functions such as histone modification, telomere maintenance, and chromatin remodeling for cow; lipid catabolism, oxidative stress, and vitamin metabolism for donkey; and terms related to chemokine receptor interaction, leukocytes migration, and transcriptional regulation in response to stress for goat. For miRNA targets, shared terms emerged as the main functions for all the species: immunity modulation, protein synthesis, cellular cycle regulation, transmembrane exchanges, and ion channels. Moreover, donkey and goat showed additional terms related to epigenetic modification and DNA maintenance. Our results showed a potential mEVs immune regulatory purpose through their RNA cargo, although in vivo validation studies are necessary.

The horse is an optimal model organism for studying the genomic response to exercise-induced stress, due to its natural aptitude for athletic performance and the relative homogeneity of its genetic and environmental backgrounds. Here, we applied RNA-sequencing analysis through the use of SOLiD technology in an experimental framework centered on exercise-induced stress during endurance races in equine athletes. We monitored the transcriptional landscape by comparing gene expression levels between animals at rest and after competition. Overall, we observed a shift from coding to non-coding regions, suggesting that the stress response involves the differential expression of not annotated regions. Notably, we observed significant post-race increases of reads that correspond to repeats, especially the intergenic and intronic L1 and L2 transposable elements. We also observed increased expression of the antisense strands compared to the sense strands in intronic and regulatory regions (1 kb up- and downstream) of the genes, suggesting that antisense transcription could be one of the main mechanisms for transposon regulation in the horse under stress conditions. We identified a large number of transcripts corresponding to intergenic and intronic regions putatively associated with new transcriptional elements. Gene expression and pathway analysis allowed us to identify several biological processes and molecular functions that may be involved with exercise-induced stress. Ontology clustering reflected mechanisms that are already known to be stress activated (e.g., chemokine-type cytokines, Toll-like receptors, and kinases), as well as “nucleic acid binding” and “signal transduction activity” functions. There was also a general and transient decrease in the global rates of protein synthesis, which would be expected after strenuous global stress. In sum, our network analysis points toward the involvement of specific gene clusters in equine exercise-induced stress, including those involved in inflammation, cell signaling, and immune interactions.