First report of junctional epidermolysis bullosa (JEB) in the Italian draft horse.

BACKGROUND: Epitheliogenesis imperfecta in horses was first recognized at the beginning of the 20th century when it was proposed that the disease could have a genetic cause and an autosomal recessive inheritance pattern. Electron microscopy studies confirmed that the lesions were characterized by a defect in the lamina propria and the disease was therefore reclassified as epidermolysis bullosa. Molecular studies targeted two mutations affecting genes involved in dermal-epidermal junction: an insertion in LAMC2 in Belgians and other draft breeds and one large deletion in LAMA3 in American Saddlebred. CASE PRESENTATION: A mechanobullous disease was suspected in a newborn, Italian draft horse foal, which presented with multifocal to coalescing erosions and ulceration on the distal extremities. Histological examination of skin biopsies revealed a subepidermal cleft formation and transmission electron microscopy demonstrated that the lamina densa of the basement membrane remained attached to the dermis. According to clinical, histological and ultrastructural findings, a diagnosis of junctional epidermolysis bullosa (JEB) was made. Genetic tests confirmed the presence of 1368insC in LAMC2 in the foal and its relatives. CONCLUSION: This is the first report of JEB in Italy. The disease was characterized by typical macroscopic, histologic and ultrastructural findings. Genetic tests confirmed the presence of the 1368insC in LAMC2 in this case: further investigations are required to assess if the mutation could be present at a low frequency in the Italian draft horse population. Atypical breeding practices are responsible in this case and played a role as odds enhancer for unfavourable alleles. Identification of carriers is fundamental in order to prevent economic loss for the horse industry.

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.

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

BACKGROUND: Copy 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. RESULTS: Analysis 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. CONCLUSIONS: Comparative 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.

Y-Chromosome Haplotype Report among Eight Italian Horse Breeds.

Horse domestication and breed selection processes have profoundly influenced the development and transformation of human society and civilization over time. Therefore, their origin and history have always attracted much attention. In Italy, several local breeds have won prestigious awards thanks to their unique traits and socio-cultural peculiarities. Here, for the first time, we report the genetic variation of three loci of the male-specific region of the Y chromosome (MSY) of four local breeds and another one (Lipizzan, UNESCO) well-represented in the Italian Peninsula. The analysis also includes data from three Sardinian breeds and another forty-eight Eurasian and Mediterranean horse breeds retrieved from GenBank for comparison. Three haplotypes (HT1, HT2, and HT3) were found in Italian stallions, with different spatial distributions between breeds. HT1 (the ancestral haplotype) was frequent, especially in Bardigiano and Monterufolino, HT2 (Neapolitan/Oriental wave) was found in almost all local breeds, and HT3 (Thoroughbred wave) was detected in Maremmano and two Sardinian breeds (Sardinian Anglo-Arab and Sarcidano). This differential distribution is due to three paternal introgressions of imported stallions from foreign countries to improve local herds; however, further genetic analyses are essential to reconstruct the genetic history of native horse breeds, evaluate the impact of selection events, and enable conservation strategies.

Exploring the Italian equine gene pool via high-throughput genotyping.

Introduction: The Italian peninsula is in the center of the Mediterranean area, and historically it has been a hub for numerous human populations, cultures, and also animal species that enriched the hosted biodiversity. Horses are no exception to this phenomenon, with the peculiarity that the gene pool has been impacted by warfare and subsequent "colonization". In this study, using a comprehensive dataset for almost the entire Italian equine population, in addition to the most influential cosmopolitan breeds, we describe the current status of the modern Italian gene pool. Materials and Methods: The Italian dataset comprised 1,308 individuals and 22 breeds genotyped at a 70 k density that was merged with publicly available data to facilitate comparison with the global equine diversity. After quality control and supervised subsampling to ensure consistency among breeds, the merged dataset with the global equine diversity contained data for 1,333 individuals from 54 populations. Multidimensional scaling, admixture, gene flow, and effective population size were analyzed. Results and Discussion: The results show that some of the native Italian breeds preserve distinct gene pools, potentially because of adaptation to the different geographical contexts of the peninsula. Nevertheless, the comparison with international breeds highlights the presence of strong gene flow from renowned breeds into several Italian breeds, probably due to historical introgression. Coldblood breeds with stronger genetic identity were indeed well differentiated from warmblood breeds, which are highly admixed. Other breeds showed further peculiarities due to their breeding history. Finally, we observed some breeds that exist more on cultural, traditional, and geographical point of view than due to actual genetic distinctiveness.

Genome-wide epigenetic modifications in sports horses during training as an adaptation phenomenon.

With his bicentennial breeding history based on athletic performance, the Thoroughbred horse can be considered the equine sport breed. Although genomic and transcriptomic tools and knowledge are at the state of the art in equine species, the epigenome and its modifications in response to environmental stimuli, such as training, are less studied. One of the major epigenetic modifications is cytosine methylation at 5' of DNA molecules. This crucial biochemical modification directly mediates biological processes and, to some extent, determines the organisms' phenotypic plasticity. Exercise indeed affects the epigenomic state, both in humans and in horses. In this study, we highlight, with a genome-wide analysis of methylation, how the adaptation to training in the Thoroughbred can modify the methylation pattern throughout the genome. Twenty untrained horses, kept under the same environmental conditions and sprint training regimen, were recruited, collecting peripheral blood at the start of the training and after 30 and 90 days. Extracted leukocyte DNA was analyzed with the methylation content sensitive enzyme ddRAD (MCSeEd) technique for the first time applied to animal cells. Approximately one thousand differently methylated genomic regions (DMRs) and nearby genes were called, revealing that methylation changes can be found in a large part of the genome and, therefore, referable to the physiological adaptation to training. Functional analysis via GO enrichment was also performed. We observed significant differences in methylation patterns throughout the training stages: we hypothesize that the methylation profile of some genes can be affected early by training, while others require a more persistent stimulus.

Unlocking Horse Y Chromosome Diversity.

The present equine genetic variation mirrors the deep influence of intensive breeding programs during the last 200 years. Here, we provide a comprehensive current state of knowledge on the trends and prospects on the variation in the equine male-specific region of the Y chromosome (MSY), which was assembled for the first time in 2018. In comparison with the other 12 mammalian species, horses are now the most represented, with 56 documented MSY genes. However, in contrast to the high variability in mitochondrial DNA observed in many horse breeds from different geographic areas, modern horse populations demonstrate extremely low genetic Y-chromosome diversity. The selective pressures employed by breeders using pedigree data (which are not always error-free) as a predictive tool represent the main cause of this lack of variation in the Y-chromosome. Nevertheless, the detailed phylogenies obtained by recent fine-scaled Y-chromosomal genotyping in many horse breeds worldwide have contributed to addressing the genealogical, forensic, and population questions leading to the reappraisal of the Y-chromosome as a powerful genetic marker to avoid the loss of biodiversity as a result of selective breeding practices, and to better understand the historical development of horse breeds.

Metabolic and Biomolecular Changes Induced by Incremental Long-Term Training in Young Thoroughbred Racehorses during First Workout Season.

Training has a huge effect on physiological homeostasis. The Thoroughbred racehorse is a valid animal model to investigate such changes for training schedule fine-tuning. As happens in human athletes, it is hypothesized that biochemical and immune response changes and related biomolecular variations could be induced by training programs. The aim of this study was to investigate, for the first time, the long-term metabolic and biomolecular modifications in young untrained Thoroughbred racehorses in the first 4-month timeframe training period. Twenty-nine clinically healthy, untrained, two-year-old Thoroughbred racehorses were followed during their incremental 4-month sprint exercise schedule. Blood collection was performed once a month, five times (T-30, T0, T30, T60, and T90). For each sample, lactate concentration, plasma cell volume (PCV), and hematobiochemical parameters (glucose, urea, creatinine, aspartate aminotransferase (AST), γ-glutamyltransferase (GGT), alkaline phosphatase (ALP), total bilirubin (Tbil), lactate dehydrogenase (LDH), creatine kinase (CK), cholesterol, triglycerides, albumin (Alb), total proteins (TPs), phosphorus (P), calcium (Ca2+), magnesium (Mg), sodium (Na+), potassium (K-), and chloride (Cl)) were determined. At T-30 and T90, serum protein electrophoresis (SPE), serum amyloid A (SAA), and real-time qPCR were performed on all samples to evaluate the expression of key genes and cytokines related to inflammatory and Th2 immunity responses: Interleukin-4 (IL-4), Interleukin-6 (IL-6), Interleukin-10 (IL-10), Interleukin-1ß (IL-1ß), Octamer-Binding Transcription Factor 1 (OCT1), B-cell lymphoma/leukemia 11A (BCL11A). Statistical analysis was performed (ANOVA and t test, p < 0.05). Significant modifications were identified compared with T-30 for PCV, glucose, triglycerides, cholesterol, lactate, urea, creatinine, Tbil, ALP, LDH, Na+, K-, Ca2+, SAA, TPs, SPE, IL-6, IL-4, Oct-1, and BCL11A. In conclusion, the first long-term training period was found to induce fundamental systemic changes in untrained Thoroughbreds.

Genetic Parameters and Inbreeding Effect of Morphological Traits in Sardinian Anglo Arab Horse.

The purpose of this study was to estimate the heritability and genetic correlations of four biometric measurements and an overall score (OS) in the Sardinian Anglo-Arab horse (SAA); moreover, the effect of inbreeding on these traits was investigated. A dataset with 43,624 horses (27,052 females and 16,572 males) was provided by the Agricultural Research Agency of Sardinia (AGRIS). Cannon bone circumference (BC), chest girth (CG), shoulder length (SL), and withers height (WH) were measured on 6033 SAA horses born in Sardinia between 1967 and 2005; beside the measurements, an overall score (OS) was taken comparing the morphology of each horse to an "ideal type" that is scored out of 100. The mean value is 20.5 cm for BC, 185.9 cm for CG, 67.6 cm for SL, 160.8 cm for WH, and 73.2 for the OS. The heritability estimates ranged from 0.78 to 0.23. The results allow to foresee high genetic progress through the breeding programs. The most affected trait by the inbreeding rate seems to only be the withers height.

Gallop Racing Shifts Mature mRNA towards Introns: Does Exercise-Induced Stress Enhance Genome Plasticity?

Physical exercise is universally recognized as stressful. Among the "sport species", the horse is probably the most appropriate model for investigating the genomic response to stress due to the homogeneity of its genetic background. The aim of this work is to dissect the whole transcription modulation in Peripheral Blood Mononuclear Cells (PBMCs) after exercise with a time course framework focusing on unexplored regions related to introns and intergenic portions. PBMCs NGS from five 3 year old Sardinian Anglo-Arab racehorses collected at rest and after a 2000 m race was performed. Apart from differential gene expression ascertainment between the two time points the complexity of transcription for alternative transcripts was identified. Interestingly, we noted a transcription shift from the coding to the non-coding regions. We further investigated the possible causes of this phenomenon focusing on genomic repeats, using a differential expression approach and finding a strong general up-regulation of repetitive elements such as LINE. Since their modulation is also associated with the "exonization", the recruitment of repeats that act with regulatory functions, suggesting that there might be an active regulation of this transcriptional shift. Thanks to an innovative bioinformatic approach, our study could represent a model for the transcriptomic investigation of stress.

Mitochondrial DNA Survey Reveals the Lack of Accuracy in Maremmano Horse Studbook Records.

* Correspondence: andrea [...].

Elo Method and Race Traits: A New Integrated System for Sport Horse Genetic Evaluation.

This first survey on Sardinian Anglo-Arab horse (SAA) race traits highlights important aspects for the breeding purpose of this population. The heritability of the race traits were estimated through a trivariate model; the estimates were 0.39, 0.33, and 0.30 for the number of placings, total earnings and Elo rating, respectively. The genetic progress could be improved by using an MT genetic evaluation of stallions and mares, combining information from competition traits.

Genetic Variability and Population Structure in the Sardinian Anglo-Arab Horse.

The Sardinian Anglo Arab (SAA) is a famous horse breed in Italy, with a significant historical background in the island of Sardinia. The aim of the study is to perform an evaluation of genetic variability in SAA using pedigree and mitochondrial data. In the current population, pedigree completeness was observed to be close to 100%, while the inbreeding coefficient and the average relatedness were lower than 3%. The ratio of effective founders/numbers of ancestors was 3.68 for the whole pedigree. The effective population size (Ne) computed by an individual increase in inbreeding (Ne_1) was 456.86, the Ne on equivalent generations (Ne_2) was 184.75, and this value slightly increased to 209.31 when computed by log-regression on equivalent generations (Ne_3). These results suggest the presence of crossbreeding and bottleneck phenomena, and they were compared with other Italian horses (reported in bibliography) to present the SAA among the Italian horse breeds scenario. Furthermore, the noteworthy mitochondrial variability reflects the use of a considerable number of founder mares; the contribution of L lineage was very important, probably because of the re-colonization from the Iberian Peninsula after the Last Glacial Maximum.

A Genetic Window on Sardinian Native Horse Breeds through Uniparental Molecular Systems.

Sardinia, an island located to the west of Italy in the Mediterranean Sea, boasts three native horse breeds: Giara, Sarcidano, and Sardinian Anglo-Arab. Here, we have investigated for the first time three loci of the non-recombining region of the Y chromosome (NRY) in 34 stallions from these breeds and performed a phylogenetic analysis of the maternal relationships among 178 previously published mitochondrial control regions. We found that the current NRY diversity of Sardinian horse breeds is linked to three haplotypes (HT), all identified within Sarcidano. Each breed showed a typical HT: HT1 (ancestral) was the most represented in Sarcidano, HT2 (Neapolitan/Oriental wave) in Giara, and HT3 (Thoroughbred wave) in Sardinian Anglo-Arab. The specificity of each haplotype suggests the influence of independent breeding strategies and the effect of genetic drift in each Sardinian population. The female counterpart, extended to 178 horses, showed a low genetic variability and a common maternal origin for Giara and Sarcidano. The higher variability of the Sardinian Anglo-Arab indicates multiple mare lineages in its current population. Further genetic analyses will be crucial to understand the paternal history of male horses, preserve the endangered mares' and stallions' lineages, and improve the enhancement of autochthonous genetic resources on this island.

Exercise-induced up-regulation of MMP-1 and IL-8 genes in endurance horses.

BACKGROUND: The stress response is a critical factor in the training of equine athletes; it is important for performance and for protection of the animal against physio-pathological disorders.In this study, the molecular mechanisms involved in the response to acute and strenuous exercise were investigated using peripheral blood mononuclear cells (PBMCs). RESULTS: Quantitative real-time PCR (qRT-PCR) was used to detect modifications in transcription levels of the genes for matrix metalloproteinase-1 (MMP-1) and interleukin 8 (IL-8), which were derived from previous genome-wide expression analysis. Significant up-regulation of these two genes was found in 10 horses that had completed a race of 90-120 km in a time-course experimental design. CONCLUSION: These results suggest that MMP-1 and IL-8 are both involved in the exercise-induced stress response, and this represents a starting point from which to understand the adaptive responses to this phenomenon.

Equine Adipose-Derived Mesenchymal Stromal Cells Release Extracellular Vesicles Enclosing Different Subsets of Small RNAs.

BACKGROUND: Equine adipose-derived mesenchymal stromal cells (e-AdMSC) exhibit attractive proregenerative properties strongly related to the delivery of extracellular vesicles (EVs) that enclose different kinds of molecules including RNAs. In this study, we investigated small RNA content of EVs produced by e-AdMSC with the aim of speculating on their possible biological role. METHODS: EVs were obtained by ultracentrifugation of the conditioned medium of e-AdMSC of 4 subjects. Transmission electron microscopy and scanning electron microscopy were performed to assess their size and nanostructure. RNA was isolated, enriched for small RNAs (<200 nt), and sequenced by Illumina technology. After bioinformatic analysis with state-of-the-art pipelines for short sequences, mapped reads were used to describe EV RNA cargo, reporting classes, and abundances. Enrichment analyses were performed to infer involved pathways and functional categories. RESULTS: Electron microscopy showed the presence of vesicles ranging in size from 30 to 300 nm and expressing typical markers. RNA analysis revealed that ribosomal RNA was the most abundant fraction, followed by small nucleolar RNAs (snoRNAs, 13.67%). Miscellaneous RNA (misc_RNA) reached 4.57% of the total where Y RNA, RNaseP, and vault RNA represented the main categories. miRNAs were sequenced at a lower level (3.51%) as well as protein-coding genes (1.33%). Pathway analyses on the protein-coding fraction revealed a significant enrichment for the "ribosome" pathway followed by "oxidative phosphorylation." Gene Ontology analysis showed enrichment for terms like "extracellular exosome," "organelle envelope," "RNA binding," and "small molecule metabolic process." The miRNA target pathway analysis revealed the presence of "signaling pathways regulating pluripotency of stem cells" coherent with the source of the samples. CONCLUSION: We herein demonstrated that e-AdMSC release EVs enclosing different subsets of small RNAs that potentially regulate a number of biological processes. These findings shed light on the role of EVs in the context of MSC biology.

Exercise induced stress in horses: selection of the most stable reference genes for quantitative RT-PCR normalization.

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. In order to investigate the molecular mechanisms underlying this process, several studies have been conducted that take advantage of microarray and quantitative real-time PCR (qRT-PCR) technologies to analyse the expression of candidate genes involved in the cellular stress response. Appropriate application of qRT-PCR, however, requires the use of reference genes whose level of expression is not affected by the test, by general physiological conditions or by inter-individual variability. RESULTS: The expression of nine potential reference genes was evaluated in lymphocytes of ten endurance horses during strenuous exercise. These genes were tested by qRT-PCR and ranked according to the stability of their expression using three different methods (implemented in geNorm, NormFinder and BestKeeper). Succinate dehydrogenase complex subunit A (SDHA) and hypoxanthine phosphoribosyltransferase (HPRT) always ranked as the two most stably expressed genes. On the other hand, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), transferrin receptor (TFRC) and ribosomal protein L32 (RPL32) were constantly classified as the less reliable controls. CONCLUSION: This study underlines the importance of a careful selection of reference genes for qRT-PCR studies of exercise induced stress in horses. Our results, based on different algorithms and analytical procedures, clearly indicate SDHA and HPRT as the most stable reference genes of our pool.

Circulating miRNAs as Putative Biomarkers of Exercise Adaptation in Endurance Horses.

Endurance exercise induces metabolic adaptations and has recently been reported associated with the modulation of a particular class of small noncoding RNAs, microRNAs, that act as post-transcriptional regulators of gene expression. Released into body fluids, they termed circulating miRNAs, and they have been recognized as more effective and accurate biomarkers than classical serum markers. This study examined serum profile of miRNAs through massive parallel sequencing in response to prolonged endurance exercise in samples obtained from four competitive Arabian horses before and 2 h after the end of competition. MicroRNA identification, differential gene expression (DGE) analysis and a protein-protein interaction (PPI) network showing significantly enriched pathways of target gene clusters, were assessed and explored. Our results show modulation of more than 100 miRNAs probably arising from tissues involved in exercise responses and indicating the modulation of correlated processes as muscle remodeling, immune and inflammatory responses. Circulating miRNA high-throughput sequencing is a promising approach for sports medicine for the discovery of putative biomarkers for predicting risks related to prolonged activity and monitoring metabolic adaptations.

Analysis of peripheral blood mononuclear cells gene expression in endurance horses by cDNA-AFLP technique.

The knowledge of molecular mechanisms of stress response in athlete horses can allow us to plan an appropriate and high-grade training to obtain better performance and to preserve horse welfare. It is well known that excessive muscular exercise can lead to a number of responses which may be associated with modification of the mRNA levels for a number of metabolic genes such as those involved in the immune response. In the present study cDNA-AFLP technique was applied to Arab endurance horses under stressing conditions to visualise variations of transcriptional profiles; 49 transcript derived fragments (TDFs), differentially expressed, were cloned and sequenced. Four of these showed high sequence similarity with genes probably involved in exercise-induced stress response and resulted to be not sequenced in the horse. Their modulation was confirmed by RT-PCR and the full-length transcripts were isolated by RACE-PCR. The mRNAs sequences obtained were included in the GenBank database as Equus caballus interleukin 8 (IL8), E. caballus retinoblastoma binding protein 6 mRNA (RBBP6), E. caballus eukaryotic translation initiation factor 4 gamma 3 (eIF4G3) and E. caballus heat shock protein 90 (Hsp90). The expression pattern of these genes was verified in other endurance horses under stressing conditions, strengthening the hypothesis of their real involvement in exercise stress-induced response.