Michela Felicetti

Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication

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.

Genetic Diversity in the Modern Horse Illustrated from Genome-Wide SNP Data

Horses were domesticated from the Eurasian steppes 5,000–6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. FST calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.

Molecular Detection, Epidemiology, and Genetic Characterization of Novel European Field Isolates of Equine Infectious Anemia Virus

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

Analysis of copy number variants by three detection algorithms and their association with body size in horses

BackgroundCopy number variants (CNVs) have been shown to play an important role in genetic diversity of mammals and in the development of many complex phenotypic traits. The aim of this study was to perform a standard comparative evaluation of CNVs in horses using three different CNV detection programs and to identify genomic regions associated with body size in horses.ResultsAnalysis was performed using the Illumina Equine SNP50 genotyping beadchip for 854 horses. CNVs were detected by three different algorithms, CNVPartition, PennCNV and QuantiSNP. Comparative analysis revealed 50 CNVs that affected 153 different genes mainly involved in sensory perception, signal transduction and cellular components. Genome-wide association analysis for body size showed highly significant deleted regions on ECA1, ECA8 and ECA9. Homologous regions to the detected CNVs on ECA1 and ECA9 have also been shown to be correlated with human height.ConclusionsComparative analysis of CNV detection algorithms was useful to increase the specificity of CNV detection but had certain limitations dependent on the detection tool. GWAS revealed genome-wide associated CNVs for body size in horses.

RNA Sequencing of the Exercise Transcriptome in Equine Athletes

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.

Quantification of Equid herpesvirus 5 DNA in clinical and necropsy specimens collected from a horse with equine multinodular pulmonary fibrosis

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.

Athletic humans and horses: Comparative analysis of interleukin-6 (IL-6) and IL-6 receptor (IL-6R) expression in peripheral blood mononuclear cells in trained and untrained subjects at rest

BackgroundHorses and humans share a natural proclivity for athletic performance. In this respect, horses can be considered a reference species in studies designed to optimize physical training and disease prevention. In both species, interleukin-6 (IL-6) plays a major role in regulating the inflammatory process induced during exercise as part of an integrated metabolic regulatory network. The aim of this study was to compare IL-6 and IL-6 receptor (IL-6R) mRNA expression in peripheral blood mononuclear cells (PBMCs) in trained and untrained humans and horses.ResultsNine highly trained male swimmers (training volume: 21.6 ± 1.7 h/wk in 10-12 sessions) were compared with two age-matched control groups represented by eight lightly trained runners (training volume: 6.4 ± 2.6 h/wk in 3-5 sessions) and nine untrained subjects. In addition, eight trained horses (training volume: 8.0 ± 2.1 h/wk in 3-4 sessions) were compared with eight age-matched sedentary mares. In humans, IL-6 mRNA levels in PBMCs determined by quantitative reverse transcription-polymerase chain reaction were significantly higher in highly trained subjects, whereas IL-6R expression did not differ among groups. In horses, transcripts of both IL-6 and IL-6R were significantly up-regulated in the trained group.ConclusionsUp-regulation of IL-6R expression in PBMCs in horses could reflect a mechanism that maintains an adequate anti-inflammatory environment at rest through ubiquitous production of anti-inflammatory cytokines throughout the body. These findings suggest that the system that controls the inflammatory response in horses is better adapted to respond to exercise than that in humans.

Genetic diversity in the Maremmano horse and its relationship with other European horse breeds

The Maremmano is an Italian warmblood horse breed from central Italy. We characterized the genetic diversity and the degree of admixture in Maremmano in comparison to 14 other European horse breeds using 30 microsatellites. Between-breed diversity explained about 9 per cent of the total genetic diversity. Cluster analysis, genetic distances and genetic differentiation coefficients showed a close relationship of Maremmano with Hanoverian and Lusitano in accordance with breed history.

Microarray analysis after strenuous exercise in peripheral blood mononuclear cells of endurance horses

It is known that moderate physical activity may have beneficial effects on health, whereas strenuous effort induces a state resembling inflammation. The molecular mechanisms underlying the cellular response to exercise remain unclear, although it is clear that the immune system plays a key role. It has been hypothesized that the physio-pathological condition that develops in athletes subjected to heavy training is caused by derangement of cellular immune regulation. The purpose of the present study was to obtain information on endurance horse gene transcription under strenuous conditions and to identify candidate genes causing immune system derangement. We performed a wide gene expression scan, using microarray technology, on peripheral blood mononuclear cells of ten horses chosen from high-level participants in national and international endurance races. The use of three different timepoints revealed changes in gene expression when post-effort samples (T1, taken immediately after the race; and T2, taken 24 h after the race) were compared with basal sample (T0, at rest). Statistical analysis showed no differences in gene expression between T0 and T2 samples, indicating complete restoration of homeostasis by 24 h after racing, whereas T1 showed strong modulation of expression, affecting 132 genes (97 upregulated, 35 downregulated). Ingenuity pathway analysis revealed that the main mechanisms and biofunctions involved were significantly associated with immunological and inflammatory responses. Real-time PCR was performed on 26 gene products to validate the array data.

Transcription of LINE-derived sequences in exercise-induced stress in horses

A large proportion of mammalian genomes is represented by transposable elements (TE), most of them being long interspersed nuclear elements 1 (LINE-1 or L1). An increased expression of LINE-1 elements may play an important role in cellular stress-related conditions exerting drastic effects on the mammalian transcriptome. To understand the impact of TE on the known horse transcriptome, we masked the horse EST database, pointing out that the amount is consistent with other major vertebrates. A previously developed transcript-derived fragments (TDFs) dataset, deriving from exercise-stimulated horse peripheral blood mononuclear cells (PBMCs), was found to be enriched with L1 (26.8% in terms of bp). We investigated the involvement of TDFs in exercise-induced stress through bioinformatics and gene expression analysis. Results indicate that LINE-derived sequences are not only highly but also differentially expressed during physical effort, hinting at interesting scenarios in the regulation of gene expression in relation to exercise.