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 in E. przewalskii-the only remaining wild horse. Therefore, a wide range of matrilineal lineages from the extinct E. ferus underwent domestication in the Eurasian steppes during the Eneolithic period and were transmitted to modern E. caballus breeds. 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.horse mitochondrial genome | mtDNA haplogroups | origin of Equus caballus | Przewalski's horse | animal domestication

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

Background: 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.

The application of molecular diagnostic techniques along with nucleotide sequence determination to permit contemporary phylogenetic analysis of European field isolates of equine infectious anemia virus (EIAV) has not been widely reported. As a result, of extensive testing instigated following the 2006 outbreak of equine infectious anemia in Italy, 24 farms with a history of exposure to this disease were included in this study. New PCR-based methods were developed, which, especially in the case of DNA preparations from peripheral blood cells, showed excellent correlation with OIE-approved agar gel immunodiffusion (AGID) tests for identifying EIAV-infected animals. In contrast, the OIE-recommended oligonucleotide primers for EIAV failed to react with any of the Italian isolates. Similar results were also obtained with samples from four Romanian farms. In addition, for the first time complete characterization of gag genes from five Italian isolates and one Romanian isolate has been achieved, along with acquisition of extensive sequence information (86% of the total gag gene) from four additional EIAV isolates (one Italian and three Romanian). Furthermore, in another 23 cases we accomplished partial characterization of gag gene sequences in the region encoding the viral matrix protein. Analysis of this information suggested that most Italian isolates were geographically restricted, somewhat reminiscent of the "clades" described for human immunodeficiency virus type 1 (HIV-1). Collectively this represents the most comprehensive genetic study of European EIAV isolates conducted to date.

BackgroundThe climatic and cultural diversity of the Italian Peninsula triggered, over time, the development of a great variety of horse breeds, whose origin and history are still unclear. To clarify this issue, analyses on phenotypic traits and genealogical data were recently coupled with molecular screening.MethodologyTo provide a comprehensive overview of the horse genetic variability in Italy, we produced and phylogenetically analyzed 407 mitochondrial DNA (mtDNA) control-region sequences from ten of the most important Italian riding horse and pony breeds: Bardigiano, Esperia, Giara, Lipizzan, Maremmano, Monterufolino, Murgese, Sarcidano, Sardinian Anglo-Arab, and Tolfetano. A collection of 36 Arabian horses was also evaluated to assess the genetic consequences of their common use for the improvement of some local breeds.ConclusionsIn Italian horses, all previously described domestic mtDNA haplogroups were detected as well as a high haplotype diversity. These findings indicate that the ancestral local mares harbored an extensive genetic diversity. Moreover, the limited haplotype sharing (11%) with the Arabian horse reveals that its impact on the autochthonous mitochondrial gene pools during the final establishment of pure breeds was marginal, if any. The only significant signs of genetic structure and differentiation were detected in the geographically most isolated contexts (i.e. Monterufolino and Sardinian breeds). Such a geographic effect was also confirmed in a wider breed setting, where the Italian pool stands in an intermediate position together with most of the other Mediterranean stocks. However, some notable exceptions and peculiar genetic proximities lend genetic support to historical theories about the origin of specific Italian breeds.

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.

A 15-year-old Belgian gelding was referred for fever, depression, and respiratory distress. Lung biopsy revealed interstitial fibrosis consistent with chronic interstitial pneumonia. Equid herpesvirus 5 (EHV-5) DNA was detected by polymerase chain reaction (PCR) in bronchoalveolar lavage and biopsy specimens. A presumptive diagnosis of equine multinodular pulmonary fibrosis (EMPF) was made, and the horse was administered a systemic treatment with corticosteroids and antiviral drugs. Despite initial clinical improvement, 4 weeks later, the condition of the horse rapidly deteriorated, and the animal was euthanized. Postmortem examination confirmed the presumptive diagnosis of EMPF. The EHV-5 DNA load in different tissues was estimated using a quantitative real-time PCR. Lung had a remarkable viral load, higher than in other organs, especially within the pulmonary fibrotic nodules, and a linkage between high viral burden and the most severely affected tissues was observed. The results suggest that the quantitative real-time PCR is a useful tool to quantify the EHV-5 load in different organs and to understand the relationship between EHV-5 and EMPF. The bronchoalveolar lavage was determined to be a good clinical sample to estimate the EHV-5 load in lung.

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.