Genome sequence, comparative analysis, and population genetics of the domestic horse

Multiomic analysis of different horse breeds reveals that gut microbial butyrate enhances racehorse athletic performance

Gut microbes play a vital role in host physiology, but whether specific bacterial functions contribute to the exceptional athletic performance of racehorses needs to be better understood. Here, we identify an association of gut butyrate-producing bacteria with athletic performance in racehorses (Thoroughbred horse). Butyrate-producing bacteria and microbial butyrate synthesis genes were significantly enriched in the racehorse gut, and the GC-MS results confirmed this conclusion. Using a mouse model, we demonstrated that sodium butyrate is sufficient to increase treadmill run time performance. We also show that butyrate improves the host response to exercise, significantly altering muscle fibre type in skeletal muscle, and increasing muscle mitochondrial function and activity. In addition, in-depth analysis of the published data showed that the gene for the synthesis of butyrate was also significantly enriched in the gut microbes of human athletes. Overall, our study indicates that gut microbial butyrate improves run time via the gut-muscle axis, providing novel insights into gut microbial functions and paving the way for improving athletic performance by targeted gut microbiome manipulation.

Distinct cellular and molecular responses to infection in three target cell types from horses, a species naturally susceptible to Ross River virus

Our current understanding of the pathogenesis of Ross River virus (RRV) infection has been derived from murine models, which do not reproduce clinical disease as experienced by infected humans and horses. This prompted us to establish more relevant host model systems to study host-virus interactions using ex vivo peripheral blood mononuclear cells (PBMCs) and in vitro primary synovial fibroblast and epidermal keratinocyte cultures. Transcriptomic analysis revealed that the expression of the transmembrane protein matrix remodelling associated 8 (mxra8), recently found to be responsible for RRV cell entry, was downregulated in all cell types when infected with RRV, compared to mock-infected controls. Potent antiviral and inflammatory responses were generated by both synovial fibroblasts and epidermal keratinocytes upon RRV infection. Upregulation of multiple genes, inducible by double-stranded RNA, together with upregulation of toll-like receptor (TLR) tlr-3, but not tlr-7, 8 and 9, suggests possible abortive replication of RRV in these cell types and potent antiviral mechanisms. This was corroborated by virus growth kinetic studies which indicated inefficient RRV replication in synovial fibroblasts and epidermal keratinocytes. Cellular metabolic flux studies on PBMCs and synovial fibroblasts showed that RRV infected cells had reduced mitochondrial function. In addition, compared to PBMCs of seronegative horses, an enhanced antiviral state and reduced inflammation related gene expression was seen in PBMCs of seropositive horses infected with RRV. Thus, despite potent antiviral and inflammatory responses via the interferon pathway exhibited in all cell types, restricting virus growth, mitochondria capacity and function of infected cells remained negatively impacted.

Running a genetic stop sign accelerates oxygen metabolism and energy production in horses

Horses are among nature's greatest athletes, yet the ancestral molecular adaptations fueling their energy demands are poorly understood. Within a clinically important pathway regulating redox and metabolic homeostasis (NRF2/KEAP1), we discovered an ancient mutation-conserved in all extant equids-that increases mitochondrial respiration while decreasing tissue-damaging oxidative stress. This mutation is a de novo premature opal stop codon in KEAP1 that is translationally recoded into a cysteine through previously unknown mechanisms, producing an R15C mutation in KEAP1 that is more sensitive to electrophiles and reactive oxygen species. This recoding enables increased NRF2 activity, which enhances mitochondrial adenosine 5'-triphosphate production and cellular resistance to oxidative damage. Our study illustrates how recoding of a de novo stop codon, a strategy thought restricted to viruses, can facilitate adaptation in vertebrates.

CENP-A/CENP-B uncoupling in the evolutionary reshuffling of centromeres in equids

BACKGROUND

While CENP-A is the epigenetic determinant of the centromeric function, the role of CENP-B, a centromeric protein binding a specific DNA sequence, the CENP-B-box, remains elusive. In the few mammalian species analyzed so far, the CENP-B box is contained in the major satellite repeat that is present at all centromeres, with the exception of the Y chromosome. We previously demonstrated that, in the genus Equus, numerous centromeres lack any satellite repeat.

RESULTS

In four Equus species, CENP-B is expressed but does not bind the majority of satellite-based centromeres, or the satellite-free ones, while it is localized at several ancestral, now-inactive, centromeres. Centromeres lacking CENP-B are functional and recruit normal amounts of CENP-A and CENP-C. The absence of CENP-B is related to the lack of CENP-B boxes rather than to peculiar features of the protein itself. CENP-B boxes are present in a previously undescribed repeat which is not the major satellite bound by CENP-A. Comparative sequence analysis suggests that this satellite was centromeric in the equid ancestor, lost centromeric function during evolution, and gave rise to a shorter CENP-A bound repeat not containing the CENP-B box but enriched in dyad symmetries.

CONCLUSIONS

We propose that the uncoupling between CENP-B and CENP-A may have played a role in the extensive evolutionary reshuffling of equid centromeres. This study provides new insights into the complexity of centromere organization in a largely biodiverse world where the majority of mammalian species still have to be studied.

The Evolution and Role of Molecular Tools in Measuring Diversity and Genomic Selection in Livestock Populations (Traditional and Up-to-Date Insights): A Comprehensive Exploration

Distinctive molecular approaches and tools, particularly high-throughput SNP genotyping, have been applied to determine and discover SNPs, potential genes of interest, indicators of evolutionary selection, genetic abnormalities, molecular indicators, and loci associated with quantitative traits (QTLs) in various livestock species. These methods have also been used to obtain whole-genome sequencing (WGS) data, enabling the implementation of genomic selection. Genomic selection allows for selection decisions based on genomic-estimated breeding values (GEBV). The estimation of GEBV relies on the calculation of SNP effects using prediction equations derived from a subset of individuals in the reference population who possess both SNP genotypes and phenotypes for target traits. Compared to traditional methods, modern genomic selection methods offer advantages for sex-limited traits, low heritability traits, late-measured traits, and the potential to increase genetic gain by reducing generation intervals. The current availability of high-density genotyping and next-generation sequencing data allow for genome-wide scans for selection. This investigation provides an overview of the essential role of advanced molecular tools in studying genetic diversity and implementing genomic selection. It also highlights the significance of adaptive selection in light of new high-throughput genomic technologies and the establishment of selective comparisons between different genomes. Moreover, this investigation presents candidate genes and QTLs associated with various traits in different livestock species, such as body conformation, meat production and quality, carcass characteristics and composition, milk yield and composition, fertility, fiber production and characteristics, and disease resistance.

Advancing equine genomics: the development of a high density Axiom_Ashwa SNP chip for Indian horses and ponies

The unique horse and pony breeds of India are declining at an alarming rate. These horses have been integral to the Indian culture and customs for centuries and represent a valuable genetic resource. It is imperative to harness the potential of this equine genetic resource that urgently needs conservation. The study highlights the design and development of a high density SNP array, the Axiom_Ashwa to aid in the genetic analysis and conservation efforts for Indian horse and pony breeds. With 613,950 SNPs, this chip offers extensive genome coverage having an average inter-marker distance of 4 kb. The Axiom_Ashwa has been validated on a larger set of diverse indigenous samples as well as Thoroughbreds, demonstrating a high call rate of 99.4% and robustness for genotyping indigenous breeds. Linkage disequilibrium (LD) analysis showed higher average LD in Indian breeds compared to exotic breeds, suggesting a limited effective population size and recent bottlenecks. Phylogenetic and population stratification analyses using PCA and DAPC clearly distinguished horses, ponies and Thoroughbreds, confirming the efficacy of the Axiom_Ashwa chip. These findings underscore the urgent need for conservation efforts for Indian horse breeds, which have experienced significant drop in population size. The Axiom_Ashwa SNP chip offers advantages such as cost-effectiveness and high throughput, providing a more accurate genetic representation of Indian horses.

Immune horses rapidly increase antileukoproteinase and lack type I interferon secretion during mucosal innate immune responses against equine herpesvirus type 1

Equine herpesvirus type 1 (EHV-1) is one of the most prevalent respiratory pathogens in horses with a high impact on animal health worldwide. Entry of the virus into epithelial cells of the upper respiratory tract and rapid local viral replication is followed by infection of local lymphoid tissues leading to cell-associated viremia and disease progression. Pre-existing mucosal immunity has previously been shown to reduce viral shedding and prevent viremia, consequently limiting severe disease manifestations. Here, nasopharyngeal transcriptomic profiling was used to identify differentially expressed genes following EHV-1 challenge in horses with different EHV-1 immune statuses. Immune horses (n = 4) did neither develop clinical disease nor viremia and did not shed virus after experimental infection, while non-immune horses (n = 4) did all the above. RNA sequencing was performed on nasopharyngeal samples pre- and 24 hours post-infection (24hpi). At 24hpi, 109 and 44 genes were upregulated in immune horses and non-immune horses, respectively, and three genes were explored in further detail. Antileukoproteinase (SLPI) gene expression increased 2.1-fold within 24 hours in immune horses in concert with protein secretion. Interferon (IFN)-induced proteins with tetratricopeptide repeats 2 (IFIT2) and 3 (IFIT3) were upregulated in non-immune horses, corresponding with nasal IFN-α secretion and viral replication. By contrast, neither IFIT expression nor IFN-α secretion was induced by EHV-1 infection of immune horses. Transcriptomic profiling offered a tool to identify, for the first time, the role of SLPI in innate immunity against EHV-1, and further emphasized the central role of the type I IFN response in the anti-viral defense of non-immune horses.

IMPORTANCE

Equine herpesvirus type 1 (EHV-1) remains a considerable concern in the equine industry, with yearly outbreaks resulting in morbidity, mortality, and economic losses. In addition to its importance in equine health, EHV-1 is a respiratory pathogen and an alphaherpesvirus, and it may serve as a model for other viruses with similar pathogenicity or phylogeny. Large animal models allow the collection of high-volume samples longitudinally, permitting in-depth investigation of immunological processes. This study was performed on bio-banked nasopharyngeal samples from an EHV-1 infection experiment, where clinical outcomes had previously been determined. Matched nucleic acid and protein samples throughout infection permitted longitudinal quantification of the protein or related proteins of selected differentially expressed genes detected during the transcriptomic screen. The results of this manuscript identified novel innate immune pathways of the upper respiratory tract during the first 24 hours of EHV-1 infection, offering a first look at the components of early mucosal immunity that are indicative of protection.

Genome-wide population affinities and signatures of adaptation in hydruntines, sussemiones and Asian wild asses

The extremely rich palaeontological record of the horse family, also known as equids, has provided many examples of macroevolutionary change over the last ~55 Mya. This family is also one of the most documented at the palaeogenomic level, with hundreds of ancient genomes sequenced. While these data have advanced understanding of the domestication history of horses and donkeys, the palaeogenomic record of other equids remains limited. In this study, we have generated genome-wide data for 25 ancient equid specimens spanning over 44 Ky and spread across Anatolia, the Caucasus, Central Asia and Mongolia. Our dataset includes the genomes from two extinct species, the European wild ass, Equus hydruntinus, and the sussemione Equus ovodovi. We document, for the first time, the presence of sussemiones in Mongolia and their survival around ~3.9 Kya, a finding that should be considered when discussing the timing of the first arrival of the domestic horse in the region. We also identify strong spatial differentiation within the historical ecological range of Asian wild asses, Equus hemionus, and incomplete reproductive isolation in several groups yet considered as different species. Finally, we find common selection signatures at ANTXR2 gene in European, Asian and African wild asses. This locus, which encodes a receptor for bacterial toxins, shows no selection signal in E. ovodovi, but a 5.4-kb deletion within intron 7. Whether such genetic modifications played any role in the sussemione extinction remains unknown.

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.

Congenital heart defects in Arabian horses and the prospects of genetic testing: A review

Congenital heart defects (CHDs) can have profound and potentially life-threatening consequences on horses' health and performance capability. While CHDs are rare in the general horse population, the Arabian breed is disproportionately overrepresented and thus is widely suspected to be genetically predisposed. This review discusses the most common CHDs in Arabian horses, including ventricular septal defect (VSD), tetralogy of Fallot (TOF), patent duct arteriosus (PDA), tricuspid valve atresia (TVA) and atrial septal defect (ASD). This review also explores how future research into the genetic factors that likely underpin many CHDs can revolutionise the way these disorders are managed in Arabian horses.

Naturally occurring horse model of miscarriage reveals temporal relationship between chromosomal aberration type and point of lethality

Chromosomal abnormalities are a common cause of human miscarriage but rarely reported in any other species. As a result, there are currently inadequate animal models available to study this condition. Horses present one potential model since mares receive intense gynecological care. This allowed us to investigate the prevalence of chromosomal copy number aberrations in 256 products of conception (POC) in a naturally occurring model of pregnancy loss (PL). Triploidy (three haploid sets of chromosomes) was the most common aberration, found in 42% of POCs following PL over the embryonic period. Over the same period, trisomies and monosomies were identified in 11.6% of POCs and subchromosomal aberrations in 4.2%. Whole and subchromosomal aberrations involved 17 autosomes, with chromosomes 3, 4, and 20 having the highest number of aberrations. Triploid fetuses had clear gross developmental anomalies of the brain. Collectively, data demonstrate that alterations in chromosome number contribute to PL similarly in women and mares, with triploidy the dominant ploidy type over the key period of organogenesis. These findings, along with highly conserved synteny between human and horse chromosomes, similar gestation lengths, and the shared single greatest risk for PL being advancing maternal age, provide strong evidence for the first animal model to truly recapitulate many key features of human miscarriage arising due to chromosomal aberrations, with shared benefits for humans and equids.

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus)

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

Vertebrate centromeres in mitosis are functionally bipartite structures stabilized by cohesin

Centromeres are scaffolds for the assembly of kinetochores that ensure chromosome segregation during cell division. How vertebrate centromeres obtain a three-dimensional structure to accomplish their primary function is unclear. Using super-resolution imaging, capture-C, and polymer modeling, we show that vertebrate centromeres are partitioned by condensins into two subdomains during mitosis. The bipartite structure is found in human, mouse, and chicken cells and is therefore a fundamental feature of vertebrate centromeres. Super-resolution imaging and electron tomography reveal that bipartite centromeres assemble bipartite kinetochores, with each subdomain binding a distinct microtubule bundle. Cohesin links the centromere subdomains, limiting their separation in response to spindle forces and avoiding merotelic kinetochore-spindle attachments. Lagging chromosomes during cancer cell divisions frequently have merotelic attachments in which the centromere subdomains are separated and bioriented. Our work reveals a fundamental aspect of vertebrate centromere biology with implications for understanding the mechanisms that guarantee faithful chromosome segregation.

Deregulation of Metalloproteinase Expression in Gray Horse Melanoma Ex Vivo and In Vitro

The ability of human melanoma cells to switch from an epithelial to a mesenchymal phenotype contributes to the metastatic potential of disease. Metalloproteinases (MPs) are crucially involved in this process by promoting the detachment of tumor cells from the primary lesion and their migration to the vasculature. In gray horse melanoma, epithelial-mesenchymal transition (EMT) is poorly understood, prompting us to address MP expression in lesions versus intact skin by transcriptome analyses and the immunofluorescence staining (IF) of gray horse tumor tissue and primary melanoma cells. RNAseq revealed the deregulation of several MPs in gray horse melanoma and, notably, a 125-fold upregulation of matrix metalloproteinase 1 (MMP1) that was further confirmed by RT-qPCR from additional tumor material. The IF staining of melanoma tissue versus intact skin for MMP1 and tumor marker S100 revealed MMP1 expression in all lesions. The co-expression of S100 was observed at different extents, with some tumors scoring S100-negative. The IF staining of primary tumor cells explanted from the tumors for MMP1 showed that the metalloproteinase is uniformly expressed in the cytoplasm of 100% of tumor cells. Overall, the presented data point to MP expression being deregulated in gray horse melanoma, and suggest that MMP1 has an active role in gray horse melanoma by driving EMT-mediated tumor cell dissemination via the degradation of the extracellular matrix. Whilst S100 is considered a reliable tumor marker in human MM, gray horse melanomas do not seem to regularly express this protein.

Predicted genetic burden and frequency of phenotype-associated variants in the horse

Disease-causing variants have been identified for less than 20% of suspected equine genetic diseases. Whole genome sequencing (WGS) allows rapid identification of rare disease causal variants. However, interpreting the clinical variant consequence is confounded by the number of predicted deleterious variants that healthy individuals carry (predicted genetic burden). Estimation of the predicted genetic burden and baseline frequencies of known deleterious or phenotype associated variants within and across the major horse breeds have not been performed. We used WGS of 605 horses across 48 breeds to identify 32,818,945 variants, demonstrate a high predicted genetic burden (median 730 variants/horse, interquartile range: 613-829), show breed differences in predicted genetic burden across 12 target breeds, and estimate the high frequencies of some previously reported disease variants. This large-scale variant catalog for a major and highly athletic domestic animal species will enhance its ability to serve as a model for human phenotypes and improves our ability to discover the bases for important equine phenotypes.

Population genomics provide insight into ancestral relationships and diversity of the feral horses of Theodore Roosevelt National Park

Theodore Roosevelt National Park (TRNP) manages a herd of feral horses (Equus caballus) which was present on the landscape prior to the establishment of the park. The population presents a unique scenario in that it has experienced fairly intensive and well-documented management since the park's establishment, including herd size reductions, intentional introduction of diversity, and subsequent attempts to remove introduced lineages. This provides an interesting case study on the genetic effects of diverse evolutionary forces on an isolated feral population. To explore the effects of these forces and clarify the relationship of this feral herd with other horses, we used genome-wide markers to examine the population structure of a combined dataset containing common established breeds. Using the Illumina Equine 70k BeadChip, we sampled SNPs across the genome for 118 TRNP horses and evaluated the inbreeding coefficient f and runs of homozygosity (RoH). To identify breed relationships, we compared 23 representative TRNP samples with 792 horses from 35 different breeds using genomic population structure analyses. Mean f of TRNP horses was 0.180, while the mean f for all other breeds in the dataset was 0.116 (SD 0.079). RoH analysis indicates that the TRNP population has experienced recent inbreeding in a timeframe consistent with their management. With Bayesian clustering, PCA, and maximum likelihood phylogeny, TRNP horses show genetic differentiation from other breeds, likely due to isolation, historical population bottlenecks, and genetic drift. However, maximum likelihood phylogeny places them with moderate confidence (76.8%) among draft breeds, which is consistent with the known history of breeds used on early North Dakota ranches and stallions subsequently introduced to the park herd. These findings will help resolve speculation about the origins of the herd and inform management decisions for the TRNP herd.

Machine learning on alignment features for parent-of-origin classification of simulated hybrid RNA-seq

BACKGROUND

Parent-of-origin allele-specific gene expression (ASE) can be detected in interspecies hybrids by virtue of RNA sequence variants between the parental haplotypes. ASE is detectable by differential expression analysis (DEA) applied to the counts of RNA-seq read pairs aligned to parental references, but aligners do not always choose the correct parental reference.

RESULTS

We used public data for species that are known to hybridize. We measured our ability to assign RNA-seq read pairs to their proper transcriptome or genome references. We tested software packages that assign each read pair to a reference position and found that they often favored the incorrect species reference. To address this problem, we introduce a post process that extracts alignment features and trains a random forest classifier to choose the better alignment. On each simulated hybrid dataset tested, our machine-learning post-processor achieved higher accuracy than the aligner by itself at choosing the correct parent-of-origin per RNA-seq read pair.

CONCLUSIONS

For the parent-of-origin classification of RNA-seq, machine learning can improve the accuracy of alignment-based methods. This approach could be useful for enhancing ASE detection in interspecies hybrids, though RNA-seq from real hybrids may present challenges not captured by our simulations. We believe this is the first application of machine learning to this problem domain.

Relevance of pharmacogenetics and pharmacogenomics in veterinary clinical practice: A review

The recent advances in high-throughput next-generation sequencing technologies have heralded the arrival of the Big Data era. As a result, the use of pharmacogenetics in drug discovery and individualized drug therapy has transformed the field of precision medicine. This paradigm shift in drug development programs has effectively reshaped the old drug development practices, which were primarily concerned with the physiological status of patients for drug development. Pharmacogenomics bridges the gap between pharmacodynamics and pharmacokinetics, advancing current diagnostic and treatment strategies and enabling personalized and targeted drug therapy. The primary goals of pharmacogenetic studies are to improve drug efficacy and minimize toxicities, to identify novel drug targets, to estimate drug dosage for personalized medicine, and to incorporate it as a routine diagnostic for disease susceptibility. Although pharmacogenetics has numerous applications in individualized drug therapy and drug development, it is in its infancy in veterinary medicine. The objective of this review is to present an overview of historical landmarks, current developments in various animal species, challenges and future perspectives of genomics in drug development and dosage optimization for individualized medicine in veterinary subjects.

Geometric morphometrics of face profile across horse breeds and within Arabian horses

Horse traits under selection are largely quantitative and affected by multiple genes. Horse face shape is an example of a continuous trait, which due to the reliance on observational assessments, is classified into; "dished", "straight", and "roman-nosed". This categorization is often inadequate to convey the full spectrum of the face shape variation especially for genetic studies. The first objective of the current study was to use geometric morphometric methods to quantitatively phenotype face shapes and examine its variation across horse breeds. The second objective was to analyze the face shape variation within Arabian horses since face shape is (1) favored, valued, and genetically selected in certain lineages (e.g. Egyptian), (2) is evaluated by registries and scored in shows, and (3) in its extreme forms pose health concerns. We digitized landmarks on lateral profile photos, particularly on the dorsal curvature of the rostrum, and subjected these landmarks to Generalized Procrustes Analysis to generate independent shape and size variables which were statistically compared across breeds and within Arabians. Horse breeds varied in nasal curvature, ranging from extremely concave to extremely convex, with over 70 % of horse breeds exhibiting intermediate concavity (i.e., straight profile). Interestingly, Arabian horses possessed the highest diversity in face profile and individuals clustered into three distinct shape sub-groups (one dished and two straight profile clusters). Our quantitative phenotyping method can be the basis of future genetic studies of facial profile within Arabian lineages as a favored traits and potentially manage its extreme forms as a likely genetic disease.

Mechanisms of Rapid Karyotype Evolution in Mammals

Chromosome reshuffling events are often a foundational mechanism by which speciation can occur, giving rise to highly derivative karyotypes even amongst closely related species. Yet, the features that distinguish lineages prone to such rapid chromosome evolution from those that maintain stable karyotypes across evolutionary time are still to be defined. In this review, we summarize lineages prone to rapid karyotypic evolution in the context of Simpson's rates of evolution-tachytelic, horotelic, and bradytelic-and outline the mechanisms proposed to contribute to chromosome rearrangements, their fixation, and their potential impact on speciation events. Furthermore, we discuss relevant genomic features that underpin chromosome variation, including patterns of fusions/fissions, centromere positioning, and epigenetic marks such as DNA methylation. Finally, in the era of telomere-to-telomere genomics, we discuss the value of gapless genome resources to the future of research focused on the plasticity of highly rearranged karyotypes.

Interdisciplinary evidence for early domestic horse exploitation in southern Patagonia

The introduction of domestic horses transformed Indigenous societies across the grasslands of Argentina, leading to the emergence of specialized horse cultures across the Southern Cone. However, the dynamics of this introduction are poorly chronicled by historic records. Here, we apply archaeozoological and biomolecular techniques to horse remains from the site of Chorrillo Grande 1 in southern Argentina. Osteological and taphonomic analyses suggest that horses were pastorally managed and used for food by Aónikenk/Tehuelche hunter-gatherers before the onset of permanent European settlement, as early as the mid-17th century. DNA-based sex identifications suggest consumption of both male and female horses, while ceramic residue also shows use of guanaco products. Sequential isotope analyses on horse dentition reveal an origin in southern Patagonia and movement of these animals between the Río Coig and Río Gallegos basins. These results reinforce emerging evidence for rapid Indigenous-mediated dispersal of horses in the Americas and demonstrate that horses catalyzed rapid economic and social transformations.

Genetic characterization of Japanese native horse breeds by genotyping variants that are associated with phenotypic traits

Concerns have been raised about the loss of genetic diversity in Japanese native horses because of their declining populations. In this study, we investigated the genetic variation of four genes, myostatin (MSTN), ligand-dependent nuclear receptor corepressor like (LCORL), doublesex and mab-3 related transcription factor 3 (DMRT3), and 5-hydroxytryptamine receptor 1A (HTR1A), which are associated with horse phenotypic traits, in six Japanese horse breeds (Hokkaido, Kiso, Noma, Misaki, Tokara, and Yonaguni). MSTN, LCORL, DMRT3, and HTR1A showed polymorphisms in the Kiso; Hokkaido and Noma; Hokkaido; and Kiso, Tokara, and Yonaguni breeds, respectively. The Misaki did not show polymorphisms in any of the genes. This study may serve as a basis for developing future breeding strategies focusing on traits in Japanese native horses.

Trio-binning of a hinny refines the comparative organization of the horse and donkey X chromosomes and reveals novel species-specific features

We generated single haplotype assemblies from a hinny hybrid which significantly improved the gapless contiguity for horse and donkey autosomal genomes and the X chromosomes. We added over 15 Mb of missing sequence to both X chromosomes, 60 Mb to donkey autosomes and corrected numerous errors in donkey and some in horse reference genomes. We resolved functionally important X-linked repeats: the DXZ4 macrosatellite and ampliconic Equine Testis Specific Transcript Y7 (ETSTY7). We pinpointed the location of the pseudoautosomal boundaries (PAB) and determined the size of the horse (1.8 Mb) and donkey (1.88 Mb) pseudoautosomal regions (PARs). We discovered distinct differences in horse and donkey PABs: a testis-expressed gene, XKR3Y, spans horse PAB with exons1-2 located in Y and exon3 in the X-Y PAR, whereas the donkey XKR3Y is Y-specific. DXZ4 had a similar ~ 8 kb monomer in both species with 10 copies in horse and 20 in donkey. We assigned hundreds of copies of ETSTY7, a sequence horizontally transferred from Parascaris and massively amplified in equids, to horse and donkey X chromosomes and three autosomes. The findings and products contribute to molecular studies of equid biology and advance research on X-linked conditions, sex chromosome regulation and evolution in equids.

Centromere innovations within a mouse species

Mammalian centromeres direct faithful genetic inheritance and are typically characterized by regions of highly repetitive and rapidly evolving DNA. We focused on a mouse species, Mus pahari, that we found has evolved to house centromere-specifying centromere protein-A (CENP-A) nucleosomes at the nexus of a satellite repeat that we identified and termed π-satellite (π-sat), a small number of recruitment sites for CENP-B, and short stretches of perfect telomere repeats. One M. pahari chromosome, however, houses a radically divergent centromere harboring ~6 mega-base pairs of a homogenized π-sat-related repeat, π-satB, that contains >20,000 functional CENP-B boxes. There, CENP-B abundance promotes accumulation of microtubule-binding components of the kinetochore and a microtubule-destabilizing kinesin of the inner centromere. We propose that the balance of pro- and anti-microtubule binding by the new centromere is what permits it to segregate during cell division with high fidelity alongside the older ones whose sequence creates a markedly different molecular composition.

Single nucleotide polymorphisms (SNPs) and indels identified from whole-genome re-sequencing of four Chinese donkey breeds

This paper represents the fundamental report of the survey of genome-wide changes of four Chinese indigenous donkey breeds, Dezhou (DZ), Guangling (GL), North China (NC), and Shandong Little donkey (SDL), and the findings will prove usefully for identification of biomarkers that perhaps predict or characterize the growth and coat color patterns. Three genomic regions in CYP3A12, TUBGCP5, and GSTA1 genes, were identified as putative selective sweeps in all researched donkey populations. The loci of candidate genes that may have contributed to the phenotypes in body size (ACSL4, MSI2, ADRA1B, and CDKL5) and coat color patterns (KITLG and TBX3) in donkey populations would be found in underlying strong selection signatures when compared between large and small donkey types, and between different coat colors. The results of the phylogenetic analysis, FST, and principal component analysis (PCA) supported that each population cannot clearly deviate from each other, showing no obvious population structure. We can conclude from the population history that the formation processes between DZS and NC, GL, and SDL are completely different. The genetic variants discovered here provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Chinese donkey breeding programs.

Progress and future perspectives of livestock genomics in India: a mini review

The field of genetics has evolved a lot after the emergence of molecular and advanced genomic technologies. The advent of Next Generation Sequencing, SNP genotyping platforms and simultaneous reduction in the cost of sequencing had opened the door to genomic research in farm animals. There are various applications of genomics in livestock, such as the use of genomic data: (i) to investigate genetic diversity and breed composition/population structure (ii) to identify genetic variants and QTLs related to economically important and ecological traits, genome-wide association studies (GWAS) and genomic signatures of selection; (iii) to enhance breeding programs by genomic selection. Compared to traditional methods, genomic selection is expected to improve selection response by increasing selection accuracy and reducing the generation interval due to early selection. Genomic selection (GS) in developed countries has led to rapid genetic gains, especially in dairy cattle, due to a well-established genetic evaluation system. Indian livestock system is still lagging behind developed nations in adopting these technologies. This review discusses the current status, challenges, and future perspectives of livestock genomics in India.

Genome-Wide Detection of Copy Number Variations Associated with Miniature Features in Horses

Copy number variations (CNVs) are crucial structural genomic variants affecting complex traits in humans and livestock animals. The current study was designed to conduct a comprehensive comparative copy number variation analysis among three breeds, Debao (DB), Baise (BS), and Warmblood (WB), with a specific focus on identifying genomic regions associated with miniature features in horses. Using whole-genome next-generation resequencing data, we identified 18,974 CNVs across 31 autosomes. Among the breeds, we found 4279 breed-specific CNV regions (CNVRs). Baise, Debao, and Warmblood displayed 2978, 986, and 895 distinct CNVRs, respectively, with 202 CNVRs shared across all three breeds. After removing duplicates, we obtained 1545 CNVRs from 26 horse genomes. Functional annotation reveals enrichment in biological functions, including antigen processing, cell metabolism, olfactory conduction, and nervous system development. Debao horses have 970 genes overlapping with CNVRs, possibly causing their small size and mountainous adaptations. We also found that the genes GHR, SOX9, and SOX11 may be responsible for the miniature features of the Debao horse by analyzing their overlapping CNVRs. Overall, this study offers valuable insights into the widespread presence of CNVs in the horse genome. The findings contribute to mapping horse CNVs and advance research on unique miniature traits observed in the Debao horse.

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.

The localization of centromere protein A is conserved among tissues

Centromeres are epigenetically specified by the histone H3 variant CENP-A. Although mammalian centromeres are typically associated with satellite DNA, we previously demonstrated that the centromere of horse chromosome 11 (ECA11) is completely devoid of satellite DNA. We also showed that the localization of its CENP-A binding domain is not fixed but slides within an about 500 kb region in different individuals, giving rise to positional alleles. These epialleles are inherited as Mendelian traits but their position can move in one generation. It is still unknown whether centromere sliding occurs during meiosis or during development. Here, we first improve the sequence of the ECA11 centromeric region in the EquCab3.0 assembly. Then, to test whether centromere sliding may occur during development, we map the CENP-A binding domains of ECA11 using ChIP-seq in five tissues of different embryonic origin from the four horses of the equine FAANG (Functional Annotation of ANimal Genomes) consortium. Our results demonstrate that the centromere is localized in the same region in all tissues, suggesting that the position of the centromeric domain is maintained during development.

Transcriptome Signature of Immature and In Vitro-Matured Equine Cumulus-Oocytes Complex

Maturation is a critical step in the development of an oocyte, and it is during this time that the oocyte advances to metaphase II (MII) of the meiotic cycle and acquires developmental competence to be fertilized and become an embryo. However, in vitro maturation (IVM) remains one of the limiting steps in the in vitro production of embryos (IVP), with a variable percentage of oocytes reaching the MII stage and unpredictable levels of developmental competence. Understanding the dynamics of oocyte maturation is essential for the optimization of IVM culture conditions and subsequent IVP outcomes. Thus, the aim of this study was to elucidate the transcriptome dynamics of oocyte maturation by comparing transcriptomic changes during in vitro maturation in both oocytes and their surrounding cumulus cells. Cumulus-oocyte complexes were obtained from antral follicles and divided into two groups: immature and in vitro-matured (MII). RNA was extracted separately from oocytes (OC) and cumulus cells (CC), followed by library preparation and RNA sequencing. A total of 13,918 gene transcripts were identified in OC, with 538 differentially expressed genes (DEG) between immature OC and in vitro-matured OC. In CC, 13,104 genes were expressed with 871 DEG. Gene ontology (GO) analysis showed an association between the DEGs and pathways relating to nuclear maturation in OC and GTPase activity, extracellular matrix organization, and collagen trimers in CC. Additionally, the follicle-stimulating hormone receptor gene (FSHR) and luteinizing hormone/choriogonadotropin receptor gene (LHCGR) showed differential expressions between CC-MII and immature CC samples. Overall, these results serve as a foundation to further investigate the biological pathways relevant to oocyte maturation in horses and pave the road to improve the IVP outcomes and the overall clinical management of equine assisted reproductive technologies (ART).

Two Novel Variants in MITF and PAX3 Associated With Splashed White Phenotypes in Horses

Mutations causing depigmentation are relatively common in Equus caballus (horse). Over 40 alleles in multiple genes are associated with increased white spotting (as of February 2023). The splashed white phenotype, a coat spotting pattern described as appearing like the horse has been splashed with white paint, was previously associated with variants in the PAX3 and MITF genes. Both genes encode transcription factors known to control melanocyte migration and pigmentation. We report two novel mutations, a stop-gain mutation in PAX3 (XM_005610643.3:c.927C>T, ECA6:11,196,181, EquCab3.0) and a missense mutation in a binding domain of MITF (NM_001163874.1:c.993A>T, ECA16:21,559,940, EquCab3.0), each with a strong association with increased depigmentation in Pura Raza Española horses (P = 1.144E-11, N = 30, P = 4.441E-16, N = 39 respectively). Using a quantitative method to score depigmentation, the PAX3 and MITF mutations were found to have average white scores of 25.50 and 24.45, respectively, compared to the average white coat spotting score of 1.89 in the control set. The functional impact for each mutation was predicted to be moderate to extreme (I-TASSER, SMART, Variant Effect Predictor, SIFT). We propose to designate the MITF mutant allele as Splashed White 9 and the PAX3 mutant allele as Splashed White 10 per convention.

Construction of an individual identification panel for horses using insertion and deletion markers

Individual identification and paternity testing are important for avoiding inbreeding in the management of small populations of wild and domestic animals. In horse racing industries, they are extremely important for identifying and registering individuals and doping control to ensure fair competition. In this study, we constructed an individual identification panel for horses by using insertion and deletion (INDEL) markers. The panel included 39 INDEL markers selected from a whole-genome INDEL database. Genotyping of 89 Thoroughbreds showed polymorphisms with minor allele frequencies (MAFs) of 0.180-0.489 in all markers. The total probability of exclusion for paternity testing, power of discrimination, and probability of identity were 0.9994271269, >0.9999999999, and 0.9999999987, respectively. The panel was applied to 13 trios (sires, dams, and foals), and no contradictions were observed in genetic inheritance among the trios. When this panel was applied to the trios (52 trios) containing false fathers, an average of 7.3 markers excluded parentage relationships. In addition, genomic DNA extracted from the urine of six horses was partially genotyped for 39 markers, and 6-28 markers were successfully genotyped. The newly constructed panel has two advantages: a low marker mutation rate compared with short tandem repeats and a genotyping procedure that is as simple as short tandem repeat typing compared with single nucleotide variant typing. This panel can be applied for individual identification, paternity determination, and urine-sample identification in Thoroughbred horses.

A randomized, triple-blinded controlled clinical study with a novel disease-modifying drug combination in equine lameness-associated osteoarthritis

Objective

This study aimed to test a novel treatment combination (TC) (equivalent to sildenafil, mepivacaine, and glucose) with disease-modifying properties compared to Celestone® bifas® (CB) in a randomized triple-blinded phase III clinical study in horses with mild osteoarthritis (OA). Joint biomarkers (reflecting the articular cartilage and subchondral bone remodelling) and clinical lameness were used as readouts to evaluate the treatment efficacy.

Methods

Twenty horses with OA-associated lameness in the carpal joint were included in the study and received either TC (n = 10) or CB (n = 10) drug intra-articularly-twice in the middle carpal joint with an interval of 2 weeks (visit 1 & 2). Clinical lameness was assessed both objectively (Lameness locator) and subjectively (visually). Synovial fluid and serum were sampled for quantification of the extracellular matrix (ECM) neo-epitope joint biomarkers represented by biglycan (BGN²⁶²) and cartilage oligomeric matrix protein (COMP¹⁵⁶). Another two weeks later clinical lameness was recorded, and serum was collected for biomarkers analysis. The overall health status was compared pre and post-intervention by interviewing the trainer.

Results

Post-intervention, SF BGN²⁶² levels significantly declined in TC (P = 0.002) and COMP¹⁵⁶ levels significantly increased in CB (P = 0.002). The flexion test scores improved in the TC compared to CB (P =0.033) and also had an improved trotting gait quality (P =0.044). No adverse events were reported.

Conclusion

This is the first clinical study presenting companion diagnostics assisting in identifying OA phenotype and evaluating the efficacy and safety of a novel disease-modifying osteoarthritic drug.

Genetic characterisation of the Connemara pony and the Warmblood horse using a within-breed clustering approach

BACKGROUND

The Connemara pony (CP) is an Irish breed that has experienced varied selection by breeders over the last fifty years, with objectives ranging from the traditional hardy pony to an agile athlete. We compared these ponies with well-studied Warmblood (WB) horses, which are also selectively bred for athletic performance but with a much larger census population. Using genome-wide single nucleotide polymorphism (SNP) and whole-genome sequencing data from 116 WB (94 UK WB and 22 European WB) and 36 CP (33 UK CP and 3 US CP), we studied the genomic diversity, inbreeding and population structure of these breeds.

RESULTS

The k-means clustering approach divided both the CP and WB populations into four genetic groups, among which the CP genetic group 1 (C1) associated with non-registered CP, C4 with US CP, WB genetic group 1 (W1) with Holsteiners, and W3 with Anglo European and British WB. Maximum and mean linkage disequilibrium (LD) varied significantly between the two breeds (mean from 0.077 to 0.130 for CP and from 0.016 to 0.370 for WB), but the rate of LD decay was generally slower in CP than WB. The LD block size distribution peaked at 225 kb for all genetic groups, with most of the LD blocks not exceeding 1 Mb. The top 0.5% harmonic mean pairwise fixation index (FST) values identified ontology terms related to cancer risk when the four CP genetic groups were compared. The four CP genetic groups were less inbred than the WB genetic groups, but C2, C3 and C4 had a lower proportion of shorter runs of homozygosity (ROH) (74 to 76% < 4 Mb) than the four WB genetic groups (80 to 85% < 4 Mb), indicating more recent inbreeding. The CP and WB genetic groups had a similar ratio of effective number of breeders (Neb) to effective population size (Ne).

CONCLUSIONS

Distinct genetic groups of individuals were revealed within each breed, and in WB these genetic groups reflected population substructure better than studbook or country of origin. Ontology terms associated with immune and inflammatory responses were identified from the signatures of selection between CP genetic groups, and while CP were less inbred than WB, the evidence pointed to a greater degree of recent inbreeding. The ratio of Neb to Ne was similar in CP and WB, indicating the influence of popular sires is similar in CP and WB.

Genome-Wide Single-Nucleotide Polymorphism-Based Genomic Diversity and Runs of Homozygosity for Selection Signatures in Equine Breeds

The horse, one of the most domesticated animals, has been used for several purposes, like transportation, hunting, in sport, or for agriculture-related works. Kathiawari, Marwari, Manipuri, Zanskari, Bhutia, Spiti, and Thoroughbred are the main breeds of horses, particularly due to their agroclimatic adaptation and role in any kind of strong physical activity, and these characteristics are majorly governed by genetic factors. The genetic diversity and phylogenetic relationship of these Indian equine breeds using microsatellite markers have been reported, but further studies exploring the SNP diversity and runs of homozygosity revealing the selection signature of breeds are still warranted. In our study, the identification of genes that play a vital role in muscle development is performed through SNP detection via the whole-genome sequencing approach. A total of 96 samples, categorized under seven breeds, and 620,721 SNPs were considered to ascertain the ROH patterns amongst all the seven breeds. Over 5444 ROH islands were mined, and the maximum number of ROHs was found to be present in Zanskari, while Thoroughbred was confined to the lowest number of ROHs. Gene enrichment of these ROH islands produced 6757 functional genes, with AGPAT1, CLEC4, and CFAP20 as important gene families. However, QTL annotation revealed that the maximum QTLs were associated with Wither's height trait ontology that falls under the growth trait in all seven breeds. An Equine SNP marker database (EqSNPDb) was developed to catalogue ROHs for all these equine breeds for the flexible and easy chromosome-wise retrieval of ROH along with the genotype details of all the SNPs. Such a study can reveal breed divergence in different climatic and ecological conditions.

High-Resolution Genotyping of Expressed Equine MHC Reveals a Highly Complex MHC Structure

The Major Histocompatibility Complex (MHC) genes play a key role in a number of biological processes, most notably in immunological responses. The MHCI and MHCII genes incorporate a complex set of highly polymorphic and polygenic series of genes, which, due to the technical limitations of previously available technologies, have only been partially characterized in non-model but economically important species such as the horse. The advent of high-throughput sequencing platforms has provided new opportunities to develop methods to generate high-resolution sequencing data on a large scale and apply them to the analysis of complex gene sets such as the MHC. In this study, we developed and applied a MiSeq-based approach for the combined analysis of the expressed MHCI and MHCII repertoires in cohorts of Thoroughbred, Icelandic, and Norwegian Fjord Horses. The approach enabled us to generate comprehensive MHCI/II data for all of the individuals (n = 168) included in the study, identifying 152 and 117 novel MHCI and MHCII sequences, respectively. There was limited overlap in MHCI and MHCII haplotypes between the Thoroughbred and the Icelandic/Norwegian Fjord horses, showcasing the variation in MHC repertoire between genetically divergent breeds, and it can be inferred that there is much more MHC diversity in the global horse population. This study provided novel insights into the structure of the expressed equine MHC repertoire and highlighted unique features of the MHC in horses.

Genome-Wide Assessment of Runs of Homozygosity by Whole-Genome Sequencing in Diverse Horse Breeds Worldwide

In the genomes of diploid organisms, runs of homozygosity (ROH), consecutive segments of homozygosity, are extended. ROH can be applied to evaluate the inbreeding situation of individuals without pedigree data and to detect selective signatures via ROH islands. We sequenced and analyzed data derived from the whole-genome sequencing of 97 horses, investigated the distribution of genome-wide ROH patterns, and calculated ROH-based inbreeding coefficients for 16 representative horse varieties from around the world. Our findings indicated that both ancient and recent inbreeding occurrences had varying degrees of impact on various horse breeds. However, recent inbreeding events were uncommon, particularly among indigenous horse breeds. Consequently, the ROH-based genomic inbreeding coefficient could aid in monitoring the level of inbreeding. Using the Thoroughbred population as a case study, we discovered 24 ROH islands containing 72 candidate genes associated with artificial selection traits. We found that the candidate genes in Thoroughbreds were involved in neurotransmission (CHRNA6, PRKN, and GRM1), muscle development (ADAMTS15 and QKI), positive regulation of heart rate and heart contraction (HEY2 and TRDN), regulation of insulin secretion (CACNA1S, KCNMB2, and KCNMB3), and spermatogenesis (JAM3, PACRG, and SPATA6L). Our findings provide insight into horse breed characteristics and future breeding strategies.

Centromere Innovations Within a Mouse Species

Mammalian centromeres direct faithful genetic inheritance and are typically characterized by regions of highly repetitive and rapidly evolving DNA. We focused on a mouse species, Mus pahari, that we found has evolved to house centromere-specifying CENP-A nucleosomes at the nexus of a satellite repeat that we identified and term π-satellite (π-sat), a small number of recruitment sites for CENP-B, and short stretches of perfect telomere repeats. One M. pahari chromosome, however, houses a radically divergent centromere harboring ~6 Mbp of a homogenized π-sat-related repeat, π-satB, that contains >20,000 functional CENP-B boxes. There, CENP-B abundance drives accumulation of microtubule-binding components of the kinetochore, as well as a microtubule-destabilizing kinesin of the inner centromere. The balance of pro- and anti-microtubule-binding by the new centromere permits it to segregate during cell division with high fidelity alongside the older ones whose sequence creates a markedly different molecular composition.

5'UTR Variant in KIT Associated with White Spotting in Horses

Mutations in KIT, a gene that influences melanoblast migration and pigmentation, often result in mammalian white spotting. As of February 2023, over 30 KIT variants associated with white spotting were documented in Equus caballus (horse). Here we report an association of increased white spotting on the skin and coat with a variant in the 5'UTR of KIT (rs1149701677: g.79,618,649A>C). Horses possessing at least one alternate allele demonstrate phenotypic characteristics similar to other KIT mutations: clear borders around unpigmented regions on the body, face, and limbs. Using a quantitative measure of depigmentation, we observed an average white score of 10.70 among individuals with rs1149701677, while the average score of the control, homozygous reference sample was 2.23 (p=1.892e-11, n=109, t-test). The rs1149701677 site has a cross-species conservation score of 3.4, one of the highest scores across the KIT 5'UTR, implying regulatory importance for this site. Ensembl also predicted a "moderately impactful" functional effect for the rs1149701677 variant. We propose that this single nucleotide variant likely alters the regulation of KIT, which in turn may disrupt melanoblast migration causing an increase in white spotting on the coat. Alternatively, the rs1149701677 variant may be in linkage with another nearby variant with an as-yet-undiscovered functional impact. We propose to term this new allele "Holiday White" or W35 based on conventional nomenclature.

The Past, Present, and Future of Equine Science

A core group of 27 equine nutritionists and physiologists joined together in the late 1960s to formally address and enhance the direction of equine research, creating the Equine Nutrition and Physiology Society. In 2003, that growing society transformed into the Equine Science Society, which now serves as the preeminent, internationally recognized scientific equine organization. In recent years, it has been appreciated that equine science encompasses a wide range of focus areas, including exercise science, nutrition, genetics, reproductive physiology, teaching and extension, production and management, and mix of other specialties, qualified as biosciences. Additionally, trainees are highly valued in the society, with the clear understanding that young people are the future of equine science. Amongst tightening budgets, equine researchers must focus on timely dissemination of high-quality research studies and development of strong, interdisciplinary, cross-species, and multi-institutional collaborations to ensure sustainability of academic research programs. With a little creativity, equine science will continue to thrive for the betterment of the horse and all involved in the equine industry.

Whole-genome identification of transposable elements reveals the equine repetitive element insertion polymorphism in Chinese horses

Transposable elements (TEs) are diverse, abundant, and complicated in genomes. They not only can drive the genome evolution process but can also act as special resources for adaptation. However, little is known about the evolutionary processes that shaped horses. In this work, 126 horse assemblages involved in most horse breeds in China were used to investigate the patterns of TE variation for the first time. By using RepeatMasker and melt software, we found that the horse-specific short interspersed repetitive elements family, equine repetitive elements (ERE1), exhibited polymorphisms in horse genomes. Phylogenetic analysis based on these ERE1 loci (minor allele frequency ≥0.05) revealed three major horse groups, namely, those in northern China, southern China, and Qinghai-Tibetan, which mirrors the result determined by SNPs to some extent. The present ERE1 family emerged ~0.26 to 1.77 Mya ago, with an activity peak at ~0.49 Mya, which matches the early stage of the horse lineage and decreases after the divergence of Equus caballus and Equus ferus przewalskii. To detect the functional ERE1(s) associated with adaptation, locus-specific branch length, genome-wide association study, and absolute allele frequency difference analyses were conducted and resulted in two common protein-coding genes annotated by candidate ERE1s. They were clustered into the vascular smooth muscle contraction (p = 0.01, EDNRA) and apelin signalling pathways (p = 0.02, NRF1). Notably, ERE1 insertion into the EDNRA gene showed a higher association with adaptation among southern China horses and other horses in 15 populations and 451 individuals (p = 4.55 e-8). Our results provide a comprehensive understanding of TE variations to analyse the phylogenetic relationships and traits relevant to adaptive evolution in horses.

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.

Functional annotation of the animal genomes: An integrated annotation resource for the horse

The genomic sequence of the horse has been available since 2009, providing critical resources for discovering important genomic variants regarding both animal health and population structures. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Due to the limited availability of functional data for the equine genome, as well as the technical limitations of short-read RNA-seq, existing annotation of the equine genome contains limited information about important aspects of gene regulation, such as alternate isoforms and regulatory elements, which are either not transcribed or transcribed at a very low level. To solve above problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. Here we detail the first comprehensive overview of gene expression and regulation in the horse, presenting 39,625 novel transcripts, 84,613 candidate cis-regulatory elements (CRE) and their target genes, 332,115 open chromatin regions genome wide across a diverse set of tissues. We showed substantial concordance between chromatin accessibility, chromatin states in different genic features and gene expression. This comprehensive and expanded set of genomics resources will provide the equine research community ample opportunities for studies of complex traits in the horse.

Inbreeding depression and durability in the North American Thoroughbred horse

The proportion of the genome containing runs of homozygosity (ROH) affects production traits in livestock populations. In European and Australasian Thoroughbreds inbreeding, quantified using ROH (F_ROH ), is associated with the probability of ever racing. Here, we measured F_ROH using 333 K SNP genotypes from 768 Thoroughbred horses born in North America to evaluate the effect of inbreeding on racing traits in that region. Among North American horses, F_ROH was not associated (p = 0.518) with the probability of ever racing but was significantly associated with the number of race starts (p = 0.002). Among raced horses, those with a 10% higher F_ROH than the mean inbreeding coefficient were predicted to have 3.5 fewer race starts compared to horses with a mean inbreeding coefficient. Considering the trend of increasing inbreeding and a decline in the average number of race starts per runner in North America, mitigating inbreeding in the population could positively influence racing durability.

A Satellite-Free Centromere in Equus przewalskii Chromosome 10

In mammals, centromeres are epigenetically specified by the histone H3 variant CENP-A and are typically associated with satellite DNA. We previously described the first example of a natural satellite-free centromere on Equus caballus chromosome 11 (ECA11) and, subsequently, on several chromosomes in other species of the genus Equus. We discovered that these satellite-free neocentromeres arose recently during evolution through centromere repositioning and/or chromosomal fusion, after inactivation of the ancestral centromere, where, in many cases, blocks of satellite sequences were maintained. Here, we investigated by FISH the chromosomal distribution of satellite DNA families in Equus przewalskii (EPR), demonstrating a good degree of conservation of the localization of the major horse satellite families 37cen and 2PI with the domestic horse. Moreover, we demonstrated, by ChIP-seq, that 37cen is the satellite bound by CENP-A and that the centromere of EPR10, the ortholog of ECA11, is devoid of satellite sequences. Our results confirm that these two species are closely related and that the event of centromere repositioning which gave rise to EPR10/ECA11 centromeres occurred in the common ancestor, before the separation of the two horse lineages.

Nerve growth factor receptors in equine synovial membranes vary with osteoarthritic disease severity

Nerve growth factor (NGF) is a neurotrophin that has been implicated in pain signaling, apoptosis, inflammation and proliferation. The resultant effects depend on interaction with two different receptors; tyrosine kinase A (TrkA) and p75^NTR . NGF increases in synovial fluid from osteoarthritic joints, and monoclonal antibody therapy is trialed to treat osteoarthritis (OA)-related pain. Investigation of the complex and somewhat contradictory signaling pathways of NGF is conducted in neural research, but has not followed through to orthopaedic studies. The objectives of this study were to compare the expression of NGF receptors and the downstream regulator BAX in synovial membranes from joints in various stages of OA. The horse was used as a model. Synovial membranes were harvested from five healthy horses postmortem and from clinical cases with spontaneous OA undergoing arthroscopic surgery for lameness. Four horses with synovitis without gross cartilage changes, four horses with synovitis and cartilage damage, and four horses with synovitis and intracarpal fractures were included. Samples were investigated by immunohistochemistry and results showed that nuclear staining of TrkA, p75^NTR and BAX increases in OA-associated synovitis. TrkA expression increased in early disease stages whereas increases in p75^NTR were most prominent in later disease stages with cartilage damage and fibrosis. Clinical significance: Suppression of NGF may result in varied effects depending on different stages of the osteoarthritic disease process.