Applied equine genetics

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.

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.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

Sequence Variant in the TRIM39-RPP21 Gene Readthrough is Shared Across a Cohort of Arabian Foals Diagnosed with Juvenile Idiopathic Epilepsy

Juvenile idiopathic epilepsy (JIE) is a self-limiting neurological disorder with a suspected genetic predisposition affecting young Arabian foals of the Egyptian lineage. The condition is characterized by tonic-clonic seizures with intermittent post-ictal blindness, in which most incidents are sporadic and unrecognized. This study aimed to identify genetic components shared across a local cohort of Arabian foals diagnosed with JIE via a combined whole genome and targeted resequencing approach: Initial whole genome comparisons between a small cohort of nine diagnosed foals (cases) and 27 controls from other horse breeds identified variants uniquely shared amongst the case cohort. Further validation via targeted resequencing of these variants, that pertain to non-intergenic regions, on additional eleven case individuals revealed a single 19bp deletion coupled with a triple-C insertion (Δ19InsCCC) within the TRIM39-RPP21 gene readthrough that was uniquely shared across all case individuals, and absent from three additional Arabian controls. Furthermore, we have confirmed recent findings refuting potential linkage between JIE and other inherited diseases in the Arabian lineage, and refuted the potential linkage between JIE and genes predisposing a similar disorder in human newborns. This is the first study to report a genetic variant to be shared in a sub-population cohort of Arabian foals diagnosed with JIE. Further evaluation of the sensitivity and specificity of the Δ19InsCCC allele within additional cohorts of the Arabian horse is warranted in order to validate its credibility as a marker for JIE, and to ascertain whether it has been introduced into other horse breeds by Arabian ancestry.

Successful ATAC-Seq From Snap-Frozen Equine Tissues

An assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) has become an increasingly popular method to assess genome-wide chromatin accessibility in isolated nuclei from fresh tissues. However, many biobanks contain only snap-frozen tissue samples. While ATAC-seq has been applied to frozen brain tissues in human, its applicability in a wide variety of tissues in horse remains unclear. The Functional Annotation of Animal Genome (FAANG) project is an international collaboration aimed to provide high quality functional annotation of animal genomes. The equine FAANG initiative has generated a biobank of over 80 tissues from two reference female animals and experiments to begin to characterize tissue specificity of genome function for prioritized tissues have been performed. Due to the logistics of tissue collection and storage, extracting nuclei from a large number of tissues for ATAC-seq at the time of collection is not always practical. To assess the feasibility of using stored frozen tissues for ATAC-seq and to provide a guideline for the equine FAANG project, we compared ATAC-seq results from nuclei isolated from frozen tissue to cryopreserved nuclei (CN) isolated at the time of tissue harvest in liver, a highly cellular homogenous tissue, and lamina, a relatively acellular tissue unique to the horse. We identified 20,000-33,000 accessible chromatin regions in lamina and 22-61,000 in liver, with consistently more peaks identified using CN isolated at time of tissue collection. Our results suggest that frozen tissues are an acceptable substitute when CN are not available. For more challenging tissues such as lamina, nuclei extraction at the time of tissue collection is still preferred for optimal results. Therefore, tissue type and accessibility to intact nuclei should be considered when designing ATAC-seq experiments.

Genomic Correlations Between the Gaits of Young Horses Measured by Accelerometry and Functional Longevity in Jumping Competition

Functional longevity is essential for the well-being of horses and the satisfaction of riders. Conventional selection using longevity breeding values calculated from competition results is not efficient because it takes too long to obtain reliable information. Therefore, the objective was to identify early criteria for selection. We assessed two types of early criteria: gait traits of young horses and QTLs. Thus, our aim was to estimate the genetic correlation between gait traits and longevity and to perform a genome-wide association study (GWAS) for longevity. Measurements of gaits by accelerometry were recorded on 1,477 show jumping horses that were 4 to 5 years old. Gait analysis provided 9 principal components describing trot, canter, and walk. Longevity estimated breeding values (EBVs) for stallions were calculated using a survival analysis of more than 900,000 years of performances by 179,448 show jumping horses born from 1981 onwards. Longevity was measured as the number of years spent in competition. Model included region and month of birth, age at first competition, year, and performance level. Longevity EBVs were deregressed to obtain weighted pseudo-performances for 1,968 stallions. Genomic data were available for 3,658 jumping horses. Seventy-eight percent of the horses measured for gaits and twenty-five percent of those measured for longevity were genotyped. A GWAS of longevity revealed no significant QTLs. Genetic parameters between each of the 9 principal components of the gait variables and longevity were evaluated with a bi-trait animal linear mixed model using single-step GBLUP analysis with the relationship matrix constructed from genomic data and genealogy (24,448 ancestors over four generations). The heritability of the gait traits varied from 0.11 to 0.44. The third principal component for trot (high lateral activity) and the first principal component for canter (high dorsoventral activity and low stride frequency) were moderately genetically correlated with higher longevity: r_g = 0.38 (0.15) and 0.28 (0.13), respectively. Our study revealed that functional longevity is a polygenic trait with no major genes. We found new correlations between longevity and gait traits. Before using gait characteristics in a selection plan, these correlations need to be understood better at the biomechanical level.

Functional contexts of adipose and gluteal muscle tissue gene co-expression networks in the domestic horse

A gene's response to an environment is tightly bound to the underlying genetic variation present in an individual's genome and varies greatly depending on the tissue it is being expressed in. Gene co-expression networks provide a mechanism to understand and interpret the collective transcriptional responses of genes. Here, we use the Camoco co-expression network framework to characterize the transcriptional landscape of adipose and gluteal muscle tissue in 83 domestic horses (Equus caballus) representing 5 different breeds. In each tissue, gene expression profiles, capturing transcriptional response due to variation across individuals, were used to build two separate, tissue-focused, genotypically-diverse gene co-expression networks. The aim of our study was to identify significantly co-expressed clusters of genes in each tissue, then compare the clusters across networks to quantify the extent that clusters were found in both networks as well as to identify clusters found in a single network. The known and unknown functions for each network were quantified using complementary, supervised and unsupervised approaches. First, supervised ontological enrichment was utilized to quantify biological functions represented by each network. Curated ontologies (GO and KEGG) were used to measure the known biological functions present in each tissue. Overall, a large percentage of terms (40.3% of GO and 41% of KEGG) were co-expressed in at least one tissue. Many terms were co-expressed in both tissues, however a small proportion of terms exhibited single tissue co-expression suggesting functional differentiation based on curated, functional annotation. To complement this, an unsupervised approach not relying on ontologies was employed. Strongly co-expressed sets of genes defined by Markov clustering identified sets of unannotated genes showing similar patterns of co-expression within a tissue. We compared gene sets across tissues and identified clusters of genes the either segregate in co-expression by tissue or exhibit high levels of co-expression in both tissues. Clusters were also integrated with GO and KEGG ontologies to identify gene sets containing previously curated annotations versus unannotated gene sets indicating potentially novel biological function. Coupling together these transcriptional datasets, we mapped the transcriptional landscape of muscle and adipose setting up a generalizable framework for interpreting gene function for additional tissues in the horse and other species.

Functionally Annotating Regulatory Elements in the Equine Genome Using Histone Mark ChIP-Seq

One of the primary aims of the Functional Annotation of ANimal Genomes (FAANG) initiative is to characterize tissue-specific regulation within animal genomes. To this end, we used chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map four histone modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in eight prioritized tissues collected as part of the FAANG equine biobank from two thoroughbred mares. Data were generated according to optimized experimental parameters developed during quality control testing. To ensure that we obtained sufficient ChIP and successful peak-calling, data and peak-calls were assessed using six quality metrics, replicate comparisons, and site-specific evaluations. Tissue specificity was explored by identifying binding motifs within unique active regions, and motifs were further characterized by gene ontology (GO) and protein–protein interaction analyses. The histone marks identified in this study represent some of the first resources for tissue-specific regulation within the equine genome. As such, these publicly available annotation data can be used to advance equine studies investigating health, performance, reproduction, and other traits of economic interest in the horse.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Investigation of Known Genetic Mutations of Arabian Horses in Egyptian Arabian Foals with Juvenile Idiopathic Epilepsy

Background

The carrier status of lavender foal syndrome (LFS), cerebellar abiotrophy (CA), severe combined immunodeficiency (SCID), and occipitoatlantoaxial malformation (OAAM1) in foals with juvenile idiopathic epilepsy (JIE) is unknown.

Hypothesis/Objectives

To determine the carrier status of LFS, CA, SCID, and OAAM1 in foals with JIE.

Animals

Ten foals with JIE.

Materials and Methods

Archived DNA samples were tested for known genetic mutations causing LFS, CA, SCID, and OAAM1. The inclusion criteria consisted of having been diagnosed with JIE by ruling out other causes of seizures in foals and supported by electroencephalographic examination.

Results

Ten Egyptian Arabian horses (5 females and 5 males) were phenotyped as foals with JIE by electroencephalography (EEG). All foals were negative for the genetic mutations that cause LFS, CA, SCID, and OAAM1 except for 1 foal that was a carrier of CA.

Conclusions and Clinical Importance

Juvenile idiopathic epilepsy of Egyptian Arabian foals and LFS appear to be phenotypically and genetically distinct disorders. There was no apparent association between JIE and LFS, CA, SCID, and OAAM1.

Single nucleotide polymorphisms for DNA typing in the domestic horse

Genetic markers are important resources for individual identification and parentage assessment. Although short tandem repeats (STRs) have been the traditional DNA marker, technological advances have led to single nucleotide polymorphisms (SNPs) becoming an attractive alternative. SNPs can be highly multiplexed and automatically scored, which allows for easier standardization and sharing among laboratories. Equine parentage is currently assessed using STRs. We obtained a publicly available SNP dataset of 729 horses representing 32 diverse breeds. A proposed set of 101 SNPs was analyzed for DNA typing suitability. The overall minor allele frequency of the panel was 0.376 (range 0.304-0.419), with per breed probability of identities ranging from 5.6 × 10^-35 to 1.86 × 10^-42 . When one parent was available, exclusion probabilities ranged from 0.9998 to 0.999996, although when both parents were available, all breeds had exclusion probabilities greater than 0.9999999. A set of 388 horses from 35 breeds was genotyped to evaluate marker performance on known families. The set included 107 parent-offspring pairs and 101 full trios. No horses shared identical genotypes across all markers, indicating that the selected set was sufficient for individual identification. All pairwise comparisons were classified using ISAG rules, with one or two excluding markers considered an accepted parent-offspring pair, two or three excluding markers considered doubtful and four or more excluding markers rejecting parentage. The panel had an overall accuracy of 99.9% for identifying true parent-offspring pairs. Our developed marker set is both present on current generation SNP chips and can be highly multiplexed in standalone panels and thus is a promising resource for SNP-based DNA typing.

Single nucleotide polymorphisms associated with osteochondrosis dissecans in Warmblood horses at different stages of training

The genetic background of osteochondrosis dissecans (OCD) has been studied for years, but the compatibility of the position of markers has not been reached between results, probably because of unknown additional effects that may influence the results, such as definition of the trait, gene–environmental interactions and the dynamics of trait development. The aim of the study was to identify single nucleotide polymorphisms (SNP) associated with the occurrence of OCD in Polish Warmblood sport breed horses in two different stages of training. Warmblood horses (87 stallions and 114 mares) were phenotyped and genotyped. Horses were X-rayed twice, at the beginning and at the end of the tests (100 days for stallions and 60 days for mares). Ten images per horse were collected using digital equipment for the fetlocks, stifles and hocks. The DNA was genotyped using the Illumina Neogen Equine Array. Statistical analysis included the Cochran–Armitage test and logistic regression assuming an additive model of inheritance. The Monte Carlo Markov Chain method was also applied to determine heritability coefficients. Nineteen and twenty SNP were identified that were significantly associated with OCD using logistic regression at the first and second stage of training, respectively. Four SNP were significant for both stages of training. The estimation of the heritability of a horse’s OCD status does not achieve the same level at different stages of training. The study on the genetic background of horse OCD should include as much detailed information on their training as possible.

Genetic Susceptibility to Rhodococcus equi

Rhodococcus equi pneumonia is a major cause of morbidity and mortality in neonatal foals. Much effort has been made to identify preventative measures and new treatments for R. equi with limited success. With a growing focus in the medical community on understanding the genetic basis of disease susceptibility, investigators have begun to evaluate the interaction of the genetics of the foal with R. equi. This review describes past efforts to understand the genetic basis underlying R. equi susceptibility and tolerance. It also highlights the genetic technology available to study horses and describes the use of this technology in investigating R. equi. This review provides readers with a foundational understanding of candidate gene approaches, single nucleotide polymorphism‐based, and copy number variant‐based genome‐wide association studies, and next generation sequencing (both DNA and RNA).

Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads

The reference assembly for the domestic horse, EquCab2, published in 2009, was built using approximately 30 million Sanger reads from a Thoroughbred mare named Twilight. Contiguity in the assembly was facilitated using nearly 315 thousand BAC end sequences from Twilight's half brother Bravo. Since then, it has served as the foundation for many genome-wide analyses that include not only the modern horse, but ancient horses and other equid species as well. As data mapped to this reference has accumulated, consistent variation between mapped datasets and the reference, in terms of regions with no read coverage, single nucleotide variants, and small insertions/deletions have become apparent. In many cases, it is not clear whether these differences are the result of true sequence variation between the research subjects' and Twilight's genome or due to errors in the reference. EquCab2 is regarded as "The Twilight Assembly." The objective of this study was to identify inconsistencies between the EquCab2 assembly and the source Twilight Sanger data used to build it. To that end, the original Sanger and BAC end reads have been mapped back to this equine reference and assessed with the addition of approximately 40X coverage of new Illumina Paired-End sequence data. The resulting mapped datasets identify those regions with low Sanger read coverage, as well as variation in genomic content that is not consistent with either the original Twilight Sanger data or the new genomic sequence data generated from Twilight on the Illumina platform. As the haploid EquCab2 reference assembly was created using Sanger reads derived largely from a single individual, the vast majority of variation detected in a mapped dataset comprised of those same Sanger reads should be heterozygous. In contrast, homozygous variations would represent either errors in the reference or contributions from Bravo's BAC end sequences. Our analysis identifies 720,843 homozygous discrepancies between new, high throughput genomic sequence data generated for Twilight and the EquCab2 reference assembly. Most of these represent errors in the assembly, while approximately 10,000 are demonstrated to be contributions from another horse. Other results are presented that include the binary alignment map file of the mapped Sanger reads, a list of variants identified as discrepancies between the source data and resulting reference, and a BED annotation file that lists the regions of the genome whose consensus was likely derived from low coverage alignments.