Whole genome detection of sequence and structural polymorphism in six diverse horses

EquCab_Finn: A new reference genome assembly for the domestic horse, Finnhorse

Finnhorse is Finland's native and national horse breed and it has genetic affinities to northern European and Asian horses. It has historical importance for agriculture, forest work and transport and as a war horse. Finnhorse has four breeding sections in the studbook and is under conservation and characterisation efforts. We sequenced and annotated the genome of a Finnhorse mare from the working horse section using PacBio and Omni-C data. This genome can complement the existing Thoroughbred reference genome (EquCab 3.0) and facilitate genetic studies of horses from northern Eurasia. We assembled 2.4 Gb of the genome with an N50 scaffold length of 83.8 Mb and the genome annotation resulted in a total of 19 748 protein coding genes of which 1200 were Finnhorse specific. The assembly has high quality and synteny with the current horse reference genome. We manually curated five genes of interest and deposited the final assembly in the European Nucleotide Archive under the accession no. PRJEB71364.

Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits

Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.

Whole-genome sequencing provides novel insights into the evolutionary history and genetic adaptation of reindeer populations in northern Eurasia

Domestic reindeer (Rangifer tarandus) play a vital role in the culture and livelihoods of indigenous people across northern Eurasia. These animals are well adapted to harsh environmental conditions, such as extreme cold, limited feed availability and long migration distances. Therefore, understanding the genomics of reindeer is crucial for improving their management, conservation and utilisation. In this study, we have generated a new genome assembly for the Fennoscandian domestic reindeer with high contiguity, making it the most complete reference genome for reindeer to date. The new genome assembly was utilised to explore genetic diversity, population structure and selective sweeps in Eurasian Rangifer tarandus populations which was based on the largest population genomic dataset for reindeer, encompassing 58 individuals from diverse populations. Phylogenetic analyses revealed distinct genetic clusters, with the Finnish wild forest reindeer (Rangifer tarandus fennicus) standing out as a unique subspecies. Divergence time estimates suggested a separation of ~ 52 thousand years ago (Kya) between the northern European Rangifer tarandus fennicus and Rangifer tarandus tarandus. Our study identified four main genetic clusters: Fennoscandian, the eastern/northern Russian and Alaskan group, the Finnish forest reindeer, and the Svalbard reindeer. Furthermore, two independent reindeer domestication processes were inferred, suggesting separate origins for the domestic Fennoscandian and eastern/northern Russian reindeer. Notably, shared genes under selection, including retroviral genes, point towards molecular domestication processes that aided adaptation of this species to diverse environments.

Genetic variation and domestication of horses revealed by 10 chromosome-level genomes and whole-genome resequencing

Understanding the genetic variations of the horse (Equus caballus) genome will improve breeding conservation and welfare. However, genetic variations in long segments, such as structural variants (SVs), remain understudied. We de novo assembled 10 chromosome-level three-dimensional horse genomes, each representing a distinct breed, and analysed horse SVs using a multi-assembly approach. Our findings suggest that SVs with the accumulation of mammalian-wide interspersed repeats related to long interspersed nuclear elements might be a horse-specific mechanism to modulate genome-wide gene regulatory networks. We found that olfactory receptors were commonly loss and accumulated deleterious mutations, but no purge of deleterious mutations occurred during horse domestication. We examined the potential effects of SVs on the spatial structure of chromatin via topologically associating domains (TADs). Breed-specific TADs were significantly enriched by breed-specific SVs. We identified 4199 unique breakpoint-resolved novel insertions across all chromosomes that account for 2.84 Mb sequences missing from the reference genome. Several novel insertions might have potential functional consequences, as 519 appeared to reside within 449 gene bodies. These genes are primarily involved in pathogen recognition, innate immune responses and drug metabolism. Moreover, 37 diverse horses were resequenced. Combining this with public data, we analysed 97 horses through a comparative population genomics approach to identify the genetic basis underlying breed characteristics using Thoroughbreds as a case study. We provide new scientific evidence for horse domestication, an understanding of the genetic mechanism underlying the phenotypic evolution of horses, and a comprehensive genetic variation resource for further genetic studies of horses.

Genes and sites under adaptation at the phylogenetic scale also exhibit adaptation at the population-genetic scale

Adaptation in protein-coding sequences can be detected from multiple sequence alignments across species or alternatively by leveraging polymorphism data within a population. Across species, quantification of the adaptive rate relies on phylogenetic codon models, classically formulated in terms of the ratio of nonsynonymous over synonymous substitution rates. Evidence of an accelerated nonsynonymous substitution rate is considered a signature of pervasive adaptation. However, because of the background of purifying selection, these models are potentially limited in their sensitivity. Recent developments have led to more sophisticated mutation-selection codon models aimed at making a more detailed quantitative assessment of the interplay between mutation, purifying, and positive selection. In this study, we conducted a large-scale exome-wide analysis of placental mammals with mutation-selection models, assessing their performance at detecting proteins and sites under adaptation. Importantly, mutation-selection codon models are based on a population-genetic formalism and thus are directly comparable to the McDonald and Kreitman test at the population level to quantify adaptation. Taking advantage of this relationship between phylogenetic and population genetics analyses, we integrated divergence and polymorphism data across the entire exome for 29 populations across 7 genera and showed that proteins and sites detected to be under adaptation at the phylogenetic scale are also under adaptation at the population-genetic scale. Altogether, our exome-wide analysis shows that phylogenetic mutation-selection codon models and the population-genetic test of adaptation can be reconciled and are congruent, paving the way for integrative models and analyses across individuals and populations.

Short Insertion and Deletion Discoveries via Whole-Genome Sequencing of 101 Thoroughbred Racehorses

Thoroughbreds are some of the most famous racehorses worldwide and are currently animals of high economic value. To understand genomic variability in Thoroughbreds, we identified genome-wide insertions and deletions (INDELs) and obtained their allele frequencies in this study. INDELs were obtained from whole-genome sequencing data of 101 Thoroughbred racehorses by mapping sequence reads to the horse reference genome. By integrating individual data, 1,453,349 and 113,047 INDELs were identified in the autosomal (1–31) and X chromosomes, respectively, while 18 INDELs were identified on the mitochondrial genome, totaling 1,566,414 INDELs. Of those, 779,457 loci (49.8%) were novel INDELs, while 786,957 loci (50.2%) were already registered in Ensembl. The sizes of diallelic INDELs ranged from −286 to +476, and the majority, 717,736 (52.14%) and 220,672 (16.03%), were 1-bp and 2-bp variants, respectively. Numerous INDELs were found to have lower frequencies of alternative (Alt) alleles. Many rare variants with low Alt allele frequencies (<0.5%) were also detected. In addition, 5955 loci were genotyped as having a minor allele frequency of 0.5 and being heterogeneous genotypes in all the horses. While short-read sequencing and its mapping to reference genome is a simple way of detecting variants, fake variants may be detected. Therefore, our data help to identify true variants in Thoroughbred horses. The INDEL database we constructed will provide useful information for genetic studies and industrial applications in Thoroughbred horses, including a gene-editing test for gene-doping control and a parentage test using INDELs for horse registration and identification.

Selection signatures for heat tolerance in Brazilian horse breeds

Since domestication, horse breeds have adapted to their environments and differentiated from one another. This paper uses two methods to detect selection signatures in 23 horse breeds, eight of which are Brazilian (610 animals), both cold-blooded and warm-blooded, from temperate and tropical regions. These animals were genotyped using the GGP Equine BeadChip and we analysed the data by Principal Component Analysis (PCA). The samples were separated into groups based on their geographical area of origin and PCA results studied. The genomic regions under selection were detected by hapFLK and PCAdapt methodologies, identifying six regions under selection with at least one Brazilian horse breed. These regions contain genes associated with heat tolerance, skin colour, body size, energy production/metabolism, genes involved in protein degradation/turnover/DNA repair, genes reducing the impact of oxidative stress/cellular repair, and transcriptional regulation. This work confirmed LCORL and NCAPG gene regions in previous studies associated with body size on Equine Chromosome Autosome 3 (ECA3). On the same ECA3, a region implicating genes linked to coat colour was identified, also previously related to heat stress. Regions with genes coding heat shock proteins were found on ECA1 and 2, and many candidate genes for oxidation-reduction which are a natural response to heat stress. However, a larger sample size and whole-genome SNPs are needed to understand better and identify new candidate regions as well as their functional relation with heat tolerance.

GWAS Identifies a Region Containing the SALL1 Gene in Variation of Pigmentation Intensity Within the Chestnut Coat Color of Horses

Chestnut coat color in horses is determined by a missense mutation within the MC1R gene. However, the intensity of the chestnut color can vary widely within individuals possessing this genotype. Here, we investigated this variation using standardized photographs of 96 horses. Each horse was ranked lightest to darkest within the cohort for phenotype by 3 blinded observers. A genome-wide association study utilizing the relative shade ranking as the phenotype and using 268 487 single-nucleotide polymorphisms (SNPs) genotyped using the Affymetrix Equine 670k array identified a single significantly associated region on chromosome 3 (P = 2.934 × 10-8). Analysis of whole-genome sequences for horses spanning the diverse range of chestnut color identified candidate SNPs within the coding sequence of the only gene in the region: SALL1. The function of SALL1 is largely unknown, though it is predicted to interact with the Hermansky-Pudlak Syndrome type 1 (HPS1) protein, which causes partial albinism in humans. However, with only one study suggesting a circumstantial influence of the SALL1 protein on pigmentation, additional work is needed to confirm this new coat color locus in larger populations and investigate the function of this protein for impacts on equine health.

Rare and common variant discovery by whole-genome sequencing of 101 Thoroughbred racehorses

The Thoroughbred breed was formed by crossing Oriental horse breeds and British native horses and is currently used in horseracing worldwide. In this study, we constructed a single-nucleotide variant (SNV) database using data from 101 Thoroughbred racehorses. Whole genome sequencing (WGS) revealed 11,570,312 and 602,756 SNVs in autosomal (1–31) and X chromosomes, respectively, yielding a total of 12,173,068 SNVs. About 6.9% of identified SNVs were rare variants observed only in one allele in 101 horses. The number of SNVs detected in individual horses ranged from 4.8 to 5.3 million. Individual horses had a maximum of 25,554 rare variants; several of these were functional variants, such as non-synonymous substitutions, start-gained, start-lost, stop-gained, and stop-lost variants. Therefore, these rare variants may affect differences in traits and phenotypes among individuals. When observing the distribution of rare variants among horses, one breeding stallion had a smaller number of rare variants compared to other horses, suggesting that the frequency of rare variants in the Japanese Thoroughbred population increases through breeding. In addition, our variant database may provide useful basic information for industrial applications, such as the detection of genetically modified racehorses in gene-doping control and pedigree-registration of racehorses using SNVs as markers.

A QTL for conformation of back and croup influences lateral gait quality in Icelandic horses

BACKGROUND

The back plays a vital role in horse locomotion, where the spine functions as a spring during the stride cycle. A complex interaction between the spine and the muscles of the back contribute to locomotion soundness, gait ability, and performance of riding and racehorses. Conformation is commonly used to select horses for breeding and performance in multiple horse breeds, where the back and croup conformation plays a significant role. The conformation of back and croup plays an important role on riding ability in Icelandic horses. However, the genes behind this trait are still unknown. Therefore, the aim of this study was to identify genomic regions associated with conformation of back and croup in Icelandic horses and to investigate their effects on riding ability. One hundred seventy-seven assessed Icelandic horses were included in the study. A genome-wide association analysis was performed using the 670 K+ Axiom Equine Genotyping Array, and the effects of different haplotypes in the top associated region were estimated for riding ability and additional conformation traits assessed during breeding field tests.

RESULTS

A suggestive quantitative trait loci (QTL) for the score of back and croup was detected on Equus caballus (ECA) 22 (p-value = 2.67 × 10- 7). Haplotype analysis revealed two opposite haplotypes, which resulted in higher and lower scores of the back and croup, respectively (p-value < 0.001). Horses with the favorable haplotype were more inclined to have a well-balanced backline with an uphill conformation and had, on average, higher scores for the lateral gaits tölt (p-value = 0.02) and pace (p-value = 0.004). This genomic region harbors three genes: C20orf85, ANKRD60 and LOC100056167. ANKRD60 is associated with body height in humans. C20orf85 and ANKRD60 are potentially linked to adolescent idiopathic scoliosis in humans.

CONCLUSIONS

Our results show that the detected QTL for conformation of back and croup is of importance for quality of lateral gaits in Icelandic horses. These findings could result in a genetic test to aid in the selection of breeding horses, thus they are of major interest for horse breeders. The results may also offer a gateway to comparative functional genomics by potentially linking both motor laterality and back inclination in horses with scoliosis in humans.

Cross-matching of allogeneic mesenchymal stromal cells eliminates recipient immune targeting

Allogeneic mesenchymal stromal cells (MSCs) have been used clinically for decades, without cross-matching, on the assumption that they are immune-privileged. In the equine model, we demonstrate innate and adaptive immune responses after repeated intra-articular injection with major histocompatibility complex (MHC) mismatched allogeneic MSCs, but not MHC matched allogeneic or autologous MSCs. We document increased peri-articular edema and synovial effusion, increased synovial cytokine and chemokine concentrations, and development of donor-specific antibodies in mismatched recipients compared with recipients receiving matched allogeneic or autologous MSCs. Importantly, in matched allogeneic and autologous recipients, but not mismatched allogeneic recipients, there was increased stromal derived factor-1 along with increased MSC concentrations in synovial fluid. Until immune recognition of MSCs can be avoided, repeated clinical use of MSCs should be limited to autologous or cross-matched allogeneic MSCs. When non-cross-matched allogeneic MSCs are used in single MSC dose applications, presensitization against donor MHC should be assessed.

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.

Metabogenomics reveals four candidate regions involved in the pathophysiology of Equine Metabolic Syndrome

An analogous condition to human metabolic syndrome, Equine Metabolic Syndrome (EMS) is defined by several clinical signs including obesity, hyperinsulinemia, and peripheral insulin dysregulation (ID). Affected horses may also exhibit hypertension, hyperlipemia and systemic inflammation. Measures of ID typically comprise the gold-standard for diagnosis in veterinary care. Yet, the dynamic nature of insulin homeostasis and complex procedures of typical assays make accurate quantification of ID and EMS challenging. This work aimed to investigate new strategies for identification of biochemical markers and correlated genes in EMS. To quantify EMS risk within this population, we utilized a composite score derived from nine common diagnostic variables. We applied a global liquid chromatography/mass spectroscopy approach (HPLC/MS) to whole plasma collected from 49 Arabian horses, resulting in 3392 high-confidence features and identification of putative metabolites in public databases. We performed a genome wide association analysis with genotypes from the 670k Affymetrix Equine SNP array utilizing EMS-correlated metabolites as phenotypes. We discovered four metabolite features significantly correlated with EMS score (P < 1.474 × 10^-5). GWAs for these features results (P = 6.787 × 10^-7, Bonferroni) identified four unique candidate regions (r² > 0.4) containing 63 genes. Significant genomic markers capture 43.52% of the variation in the original EMS score phenotype. The identified genomic loci provide insight into the pathways controlling variation in EMS and the origin of genetic predisposition to the condition. Rapid, feasible and accurate diagnostic tools derived from metabogenomics can be translated into measurable benefits in the timeline and quality of preventative management practices to preserve health in horses.

Genome Diversity and the Origin of the Arabian Horse

The Arabian horse, one of the world’s oldest breeds of any domesticated animal, is characterized by natural beauty, graceful movement, athletic endurance, and, as a result of its development in the arid Middle East, the ability to thrive in a hot, dry environment. Here we studied 378 Arabian horses from 12 countries using equine single nucleotide polymorphism (SNP) arrays and whole-genome re-sequencing to examine hypotheses about genomic diversity, population structure, and the relationship of the Arabian to other horse breeds. We identified a high degree of genetic variation and complex ancestry in Arabian horses from the Middle East region. Also, contrary to popular belief, we could detect no significant genomic contribution of the Arabian breed to the Thoroughbred racehorse, including Y chromosome ancestry. However, we found strong evidence for recent interbreeding of Thoroughbreds with Arabians used for flat-racing competitions. Genetic signatures suggestive of selective sweeps across the Arabian breed contain candidate genes for combating oxidative damage during exercise, and within the “Straight Egyptian” subgroup, for facial morphology. Overall, our data support an origin of the Arabian horse in the Middle East, no evidence for reduced global genetic diversity across the breed, and unique genetic adaptations for both physiology and conformation.