Improved reference genome for the domestic horse increases assembly contiguity and composition

Genome wide landscaping of copy number variations for horse inter-breed variability

Copy number variations (CNVs) have become widely acknowledged as a significant source of genomic variability and phenotypic variance. To understand the genetic variants in horses, CNVs from six Indian horse breeds, namely, Manipuri, Zanskari, Bhutia, Spiti, Kathiawari and Marwari were discovered using Axiom™ Equine Genotyping Array. These breeds differed in agro-climatic adaptation with distinct phenotypic characters. A total of 2668 autosomal CNVs and 381 CNV regions (CNVRs) were identified with PennCNV tool. DeepCNV was employed to re-validate to get 883 autosomal CNVs, of which 9.06% were singleton type. A total of 180 CNVRs were identified after DeepCNV filtering with the estimated length of 3.12 Kb-4.90 Mb. The functional analysis showed the majority of the CNVRs genes enriched for sensory perception and olfactory receptor activity. An Equine CNVs database, EqCNVdb (http://backlin.cabgrid.res.in/eqcnvdb/) was developed which catalogues detailed information on the horse CNVs, CNVRs and gene content within CNVRs. Also, three random CNVRs were validated with real-time polymerase chain reaction. These findings will aid in the understanding the horse genome and serve as a preliminary foundation for future CNV association research with commercially significant equine traits. The identification of CNVs and CNVRs would lead to better insights into genetic basis of important traits.

Genomic map of the functionally extinct northern white rhinoceros (Ceratotherium simum cottoni)

The northern white rhinoceros (NWR; Ceratotherium simum cottoni) is functionally extinct, with only two nonreproductive females alive. Efforts to rescue the NWR from its inevitable demise have inspired the exploration of unconventional conservation methods, including the development of induced pluripotent stem cells (iPSCs) for the in vitro generation of artificial gametes. The integrity of iPSC genomes is critical for in vitro gametogenesis to be used for assisted reproductive technologies using NWR iPSCs. We generated a chromosome-level NWR reference genome that meets or exceeds the metrics proposed by the Vertebrate Genome Project, using complementary sequencing and mapping methods. The genome represents 40 autosomes, an X and a partially resolved Y chromosome, and the mitochondrial genome. Using comparative FISH mapping, we confirmed a general gene order conservation between the NWR and horse genomes. We aligned the NWR genome with that of the southern white rhinoceros (SWR; Ceratotherium simum simum), a population that has been physically separated from the NWR for tens of thousands of years, and we found that the two subspecies are very similar on the chromosome level. Comparing long-read data from NWR iPSC lines and the fibroblast cultures used for reprogramming, we identified copy number variations that were likely to have been introduced during in vitro iPSC expansion. The NWR reference genome allows for efficient, rapid, and accurate assessment of the genomic integrity of iPSC lines to direct their differentiation. This will assist in strategies to rescue the NWR through extraordinary measures like cloning and the generation of embryos from iPSC-derived gametes.

Beyond symptomatic alignment: evaluating the integration of causal mechanisms in matching animal models with human pathotypes in osteoarthritis research

Osteoarthritis (OA) is a highly prevalent and disabling condition lacking curative treatments, with only symptomatic relief available. Recognizing OA as a heterogenous disorder with diverse aetiologies and molecular foundations underscores the need to classify patients by both phenotypes and molecular pathomechanisms (endotypes). Such stratification could enable the development of targeted therapies to surmount existing treatment barriers. From a scientific, economic, and ethical perspective, it is crucial to employ animal models that accurately represent the endotype of the target patient population, not merely their clinical symptoms. These models must also account for intrinsic and extrinsic factors, like age, sex, metabolic status, and comorbidities, which impact OA's pathogenesis and its clinical and molecular variability and can profoundly influence not only structural and symptomatic disease severity and progression but also the underlying molecular pathophysiology. The molecular definition of the OA subpopulation must also be reflected in the read-outs, as the traditional methods-macroscopic and histological scoring, along with limited gene expression profiling of established biomarkers for cartilage degradation, extracellular matrix (ECM) turnover, and synovial inflammation-are inadequate for discovering new, phenotype- and endotype-specific biomarkers or therapeutic targets. Thus, animal model characterisation should evolve to include both clinically and pathophysiologically pertinent measures of disease progression and response to treatment. This review evaluates the utility and accuracy of current animal models in OA research, focusing on their capacity to replicate the disease's pathophysiological processes.

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.

Cathelicidin antimicrobial peptides mediate immune protection in marsupial neonates

Marsupial neonates are born with immature immune systems, making them vulnerable to pathogens. While neonates receive maternal protection, they can also independently combat pathogens, although the mechanisms remain unknown. Using the sugar glider (Petaurus breviceps) as a model, we investigated immunological defense strategies of marsupial neonates. Cathelicidins-a family of antimicrobial peptides expanded in the genomes of marsupials-are highly expressed in developing neutrophils. Sugar glider cathelicidins reside in two genomic clusters, and their coordinated expression is achieved by enhancer sharing within clusters and long-range physical interactions between clusters. Functionally, cathelicidins modulate immune responses and have potent antibacterial effects, sufficient to provide protection in a mouse model of sepsis. Evolutionarily, cathelicidins have a complex history, with marsupials and monotremes uniquely retaining both clusters among tetrapods. Thus, cathelicidins are critical mediators of marsupial immunity, and their evolution may reflect the life history-specific immunological needs of these animals.

Exploring the genetic influences on equine analgesic efficacy through genome-wide association analysis of ranked pain responses

Multimodal analgesic administration is a promising strategy for mitigating side effects typically associated with analgesia; nevertheless, variation in analgesic effectiveness still poses a considerable safety concern for both horses and veterinarians. Pharmacogenomic studies have started delving into genetic influences on varying drug effectiveness and related side effects. However, current findings have narrow implications and are limited in their ability to individualize analgesic dosages in horses. Hydromorphone and detomidine were administered to a cohort of 48 horses at standardized time intervals, with dosage rates recorded. Analgesic effectiveness was scored (1-3) based on pain response to dura penetration during cerebrospinal fluid centesis. Genome-wide association (GWA) analyses identified two SNVs passing the nominal significance threshold (P < 1 ×10-5) in association with analgesic effectiveness. One SNV identified on chromosome 27 (rs1142378599) is contained within the LOC100630731 disintegrin and metalloproteinase domain-containing protein 5 gene. The second identified SNV is an intergenic variant located on chromosome 29 (rs3430772468) These SNVs accounted for 26.11 % and 31.72 % of explained variation in analgesic effectiveness respectively, with all eight of the horses with the lowest analgesic effectiveness expressing the A/C genotype at rs3430772468, with six of which also expressing the C/T genotype at rs1142872965. Whilst highlighting the multifactorial nature of analgesic efficacy, this study serves as an important step in the application of genome-wide approaches to better understand genetic factors underpinning commonly observed variation in analgesic effectiveness in horses, with the goal of tailoring analgesic dosage to minimize commonly observed side effects and improve the outcomes of equine pain management.

Extracellular Vesicle-Derived microRNA Crosstalk Between Equine Chondrocytes and Synoviocytes-An In Vitro Approach

This study describes a novel technique to analyze the extracellular vesicle (EV)-derived microRNA (miRNA) crosstalk between equine chondrocytes and synoviocytes. Donor cells (chondrocytes, n = 8; synoviocytes, n = 9) were labelled with 5-ethynyl uridine (5-EU); EVs were isolated from culture media and incubated with recipient cells (chondrocytes [n = 5] were incubated with synoviocyte-derived EVs, and synoviocytes [n = 4] were incubated with chondrocyte-derived EVs). Total RNA was extracted from recipient cells; the 5-EU-labelled RNA was recovered and sequenced. Differential expression analysis, pathway analysis, and miRNA target prediction were performed. Overall, 198 and 213 miRNAs were identified in recipient synoviocytes and chondrocytes, respectively. The top five most abundant miRNAs were similar for synoviocytes and chondrocytes (eca-miR-21, eca-miR-221, eca-miR-222, eca-miR-100, eca-miR-26a), and appeared to be linked to joint homeostasis. There were nine differentially expressed (p < 0.05) miRNAs (eca-miR-27b, eca-miR-23b, eca-miR-31, eca-miR-191a, eca-miR-199a-5p, eca-miR-143, eca-miR-21, eca-miR-181a, and eca-miR-181b) between chondrocytes and synoviocytes, which appeared to be linked to migration of cells, apoptosis, cell viability of connective tissue cell, and inflammation. In conclusion, the reported technique was effective in recovering and characterizing the EV-derived miRNA crosstalk between equine chondrocytes and synoviocytes and allowed for the identification of EV-communicated miRNA patterns potentially related to cell viability, inflammation, and joint homeostasis.

Genomic Insights into Post-Domestication Expansion and Selection of Body Size in Ponies

Horse domestication revolutionizes human civilization by transforming transportation, agriculture, and warfare patterns. Despite extensive studies on modern domestic horse origins, the intricate demographic history and genetic signatures underlying pony size remain unexplored. Here, a high-quality genome assembly of the Chinese Debao pony is presented, and 452 qualified individuals from 64 horse breeds worldwide are extensively analyzed. The authors' results reveal the conservation of ancient components in East Asian horses and close relationships between Asian horses and Western pony lineages. Genetic analyses suggest an Asian paternal origin for European pony breeds. These pony-sized horses share close genetic affinities, potentially attributed to their early expansion and adaptation to local environments. In addition, promising cis-regulatory elements influencing horse withers height by regulating genes such as RFLNA and FOXO1 are identified. Overall, this study provides insightful perspectives on the dispersal history and genetic determinants underlying body size in ponies, offering broader implications for horse population management and improvement.

Three Novel KIT Polymorphisms Found in Horses with White Coat Color Phenotypes

This paper reports three novel KIT variants likely responsible for previously unexplained white patterning phenotypes observed in three groups of horses. White spots and markings may have substantial consequences on the value and health of domesticated horses. This study aims to elucidate the genetic mechanisms underlying depigmented coat colors to aid in producing prosperous herds. Aligned whole genome sequences were manually screened to identify three polymorphisms in a family of Anglo-Arabian horses (N = 7), a family of Warmblood horses (N = 5), and a single stock-type mare with unexplained white markings. Sanger sequencing confirmed the presence of the variants, and in silico predictive programs were used to predict the functional impacts of each. We propose to term the novel variants W37, W38, and W39, respectively, per convention. The W37 polymorphism was always observed in the presence of one W35 allele, suggesting complete linkage. All three variants were predicted to alter or remove the KIT protein active domain, repressing typical protein folding and impacting pathways that upregulate pigmentation. The severe predicted impact on biological function suggests that these variants may cause increased white spotting, providing a possible explanation for the depigmentation phenotypes observed in affected individuals.

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.

Rewilded horses in European nature conservation - a genetics, ethics, and welfare perspective

In recent decades, the integration of horses (Equus ferus) in European rewilding initiatives has gained widespread popularity due to their potential for regulating vegetation and restoring natural ecosystems. However, employing horses in conservation efforts presents important challenges, which we here explore and discuss. These challenges encompass the lack of consensus on key terms inherent to conservation and rewilding, the entrenched culture and strong emotions associated with horses, low genetic diversity and high susceptibility to hereditary diseases in animals under human selection, as well as insufficient consideration for the social behaviour of horses in wild-living populations. In addition, management of wild-living horses involves intricate welfare, ethics and legislative dimensions. Anthropocentric population-control initiatives may be detrimental to horse group structures since they tend to prioritise individual welfare over the health of populations and ecosystems. To overcome these challenges, we provide comprehensive recommendations. These involve a systematic acquisition of genetic information, a focus on genetic diversity rather than breed purity and minimal veterinary intervention in wild-living populations. Further, we advise allowing for natural top-down and bottom-up control - or, if impossible, simulating this by culling or non-lethal removal of horses - instead of using fertility control for population management. We advocate for intensified collaboration between conservation biologists and practitioners and enhanced communication with the general public. Decision-making should be informed by a thorough understanding of the genetic makeup, common health issues and dynamics, and social behaviour in wild-living horse populations. Such a holistic approach is essential to reconcile human emotions associated with horses with the implementation of conservation practices that are not only effective but also sustainable for the long-term viability of functional, biodiverse ecosystems, while rehabilitating the horse as a widespread wild-living species in Europe.

A comprehensive allele specific expression resource for the equine transcriptome

BACKGROUND

Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression.

RESULTS

An equine ASE analysis was performed, using integrated Iso-seq and short-read RNA sequencing data from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues from the Functional Annotation of Animal Genomes (FAANG) project. Allele expression was quantified by haplotypes from long-read data, with 42,900 allele expression events compared. Within these events, 635 (1.48%) demonstrated ASE, with liver tissue containing the highest proportion. Genetic variants within ASE events were located in histone modified regions 64.2% of the time. Validation of allele-specific variants, using a set of 66 equine liver samples from multiple breeds, confirmed that 97% of variants demonstrated ASE.

CONCLUSIONS

This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues. Our results highlight the tissue-specific nature of allelic imbalance in the equine genome.

Genetic composition of Kazakh horses of Zhabe type evaluated by SNP genotyping

Horses are animals traditionally playing prominent role as both food source and working animals for Kazakh people. Zhabe horses are traditional type of indigenous Kazakh horses characterized by versatility and adaptation to conditions of Central Asia. The present work focuses on examination of genetic structure of Zhabe horses using SNP genotyping with addition of previously published data. Total 1038 individuals including 403 new samples of Zhabe horses and 42 sample of white horses 'Zhetysu Asyly' have been considered. DNA was extracted from hair roots using commercial DNA isolation kit and further used for analysis of SNP by Illumina iScan system with Equine80k SNP array. The analysis of population genetic parameters (expected and observed heterozygosity, linkage disequilibrium, Wright's F st ) and genetic structure (PCA, ADMIXTURE) in comparison with publicly available data on selected foreign cultivars demonstrated low between population differentiations and lack of selection factors. Genome wide association study performed for body size and weight have revealed low occurrence of SNPs with significant associations, total 57 SNPs linked to various genes with low density across all genome. The obtained results highlight difference between traditional horse breeding practices of Kazakh people and stable based breeding of foreign breeds. In contrast, the 'Zhetysu Asyly' horse breed derived from Kazakh horses demonstrate the effect of intense breeding process on the same landrace. The results provide new data on the traditional Kazakh horses of type Zhabe and will assist further studies of this original landrace.

Effect of an endothelial regulatory module on plasma proteomics in exercising horses

Elite performing exercise requires an intricate modulation of the blood pressure to support the working muscles with oxygen. We have previously identified a genomic regulatory module that associates with differences in blood pressures of importance for elite performance in racehorses. This study aimed to determine the effect of the regulatory module on the protein repertoire. We sampled plasma from 12 Coldblooded trotters divided into two endothelial regulatory module haplotype groups, a sub-elite performing haplotype (SPH) and an elite performing haplotype (EPH), each at rest and exercise. The haplotype groups and their interaction were interrogated in two analyses, i) individual paired ratio analysis for identifying differentially abundant proteins of exercise (DAPE) and interaction (DAPI) between haplotype and exercise, and ii) unpaired ratio analysis for identifying differentially abundant protein of haplotype (DAPH). The proteomics analyses revealed a widespread change in plasma protein content during exercise, with a decreased tendency in protein abundance that is mainly related to lung function, tissue fluids, metabolism, calcium ion pathway and cellular energy metabolism. Furthermore, we provide the first investigation of the proteome variation due to the interaction between exercise and related blood pressure haplotypes, which this difference was related to a faster switch to the lipoprotein and lipid metabolism during exercise for EPH. The molecular signatures identified in the present study contribute to an improved understanding of exercise-related blood pressure regulation.

Genomic insights into the genetic diversity, lateral gaits and high-altitude adaptation of Chakouyi (CKY) horses

Chakouyi (CKY) horses from the Qinghai‒Xizang Plateau are well known for their unique lateral gaits and high-altitude adaptation, but genetic mechanisms underlying these phenotypes remain unclear. This study presents a comparison of 60 newly resequenced genomes of gaited CKY horses with 139 public genomes from 19 horse breeds. Population structure analyses (admixture, PCA, and neighbor-joining tree) reveal a close genetic relationship between CKY and other highland breeds (Tibetan and Chaidamu horses). Compared with other Chinese breeds, CKY horses present reduced nucleotide diversity (θπ) and lower inbreeding (FROH coefficient), suggesting possible selective pressures. A key region on chromosome 23 (Chr23: 22.3 -22.6 Mb) is associated with the lateral gaits and harbors a highly prevalent nonsense mutation (Chr 23:22,391,254 C>A, Ser301STOP) in the DMRT3 gene, with an 88% homozygosity rate, which is strongly correlated with the distinctive gait of CKY horses. Furthermore, selection signals reveal that the EPAS1 gene is related to high-altitude adaptation, and the CAT gene contributes to altitude resilience in CKY horses. These findings suggest that preserving genetic diversity is essential for maintaining the unique gaits and high-altitude adaptations of CKY horses.

A genome-wide association study of the racing performance traits in Yili horses based on Blink and FarmCPU models

Racing performance traits are the main indicators for evaluating the performance and value of sport horses. The aim of this study was to identify the key genes for racing performance traits in Yili horses by performing a genome-wide association study (GWAS). Breeding values for racing performance traits were calculated for Yili horses (n = 827) using an animal model. Genome-wide association analysis of racing performance traits in horses (n = 236) was carried out using the Blink, and FarmCPU models in GAPIT software, and genes within the significant regions were functionally annotated. The results of GWAS showed that a total of 24 significant SNP markers (P < 6.05 × 10- 9) and 22 suggestive SNP markers (P < 1.21 × 10- 7) were identified. Among them, the Blink associated 16 significant SNP loci and FarmCPU associated 12 significant SNP loci. A total of 127 candidate genes (50 significant) were annotated. Among these, CNTN6 (motor coordination), NIPA1 (neuronal development), and DCC (dopamine pathway maturation) may be the main candidate genes affecting speed traits. SHANK2 (neuronal synaptic regulation), ISCA1 (mitochondrial protein assembly), and KCNIP4 (neuronal excitability) may be the main candidate genes affecting ranking score traits. A common locus (ECA1: 22698579) was significantly associated with racing performance traits, and the function of the genes at this locus needs to be studied in depth. These findings will provide new insights into the detection and selection of genetic variants for racing performance and will help to accelerate the genetic improvement of Yili horses.

Placental homogeneity: Characterizing transcriptional variation among equine chorioallantoic locations

The proper function of the placenta is essential for the health and growth of the fetus and the mother. The placenta relies on dynamic gene expression for its correct and timely development and function. Although numerous studies have identified genes vital for placental functions, equine placental molecular research has primarily focused on single placental locations, in sharp contrast with the broader approach in human studies. Here, we hypothesized that the molecular differences across different regions of the equine placenta are negligible because of its diffuse placental type with a macroscopic homogenous distribution of villi across the placental surface. We compared the transcriptome and stereological findings of the body, pregnant horn, and non-pregnant horn within the equine chorioallantois. Our transcriptomic analysis indicates that the variation between regions of the placenta within individuals is less than the variation observed between individuals. A low number of differentially expressed genes (DEGs) (n = 8) was identified when comparing pregnant and non-pregnant horns within the same placenta, suggesting a remarkable molecular uniformity. A higher number of DEGs was identified when comparing each horn to the body (193 DEGs comparing pregnant horn with body and 207 DEGs comparing non-pregnant horn with body). Genes with a higher expression in the body were associated with processes such as extracellular matrix synthesis and remodeling, which is relevant for placental maturation and placenta-endometrial separation at term and implies asynchrony of these processes across locations. The stereological analysis showed no differences in microcotyledonary density, and width between the locations. However, we observed a greater chorioallantoic thickness in the body and pregnant horn compared to the non-pregnant horn. Overall, our findings reveal a uniform transcriptomic profile across the placental horns, alongside a more distinct gene expression pattern between the uterine body and horns. These regional differences in gene expression suggest a different pace in the placental maturation and detachment among the placental locations.

Analyses of whole-genome sequences from 185 North American Thoroughbred horses, spanning 5 generations

Whole genome sequences (WGS) of 185 North American Thoroughbred horses were compared to quantify the number and frequency of variants, diversity of mitotypes, and autosomal runs of homozygosity (ROH). Of the samples, 82 horses were born between 1965 and 1986 (Group 1); the remaining 103, selected to maximize pedigree diversity, were born between 2000 and 2020 (Group 2). Over 14.3 million autosomal variants were identified with 4.5-5.0 million found per horse. Mitochondrial sequences associated the North American Thoroughbreds with 9 of 17 clades previously identified among diverse breeds. Individual coefficients of inbreeding, estimated from ROH, averaged 0.266 (Group 1) and 0.283 (Group 2). When SNP arrays were simulated using subsets of WGS markers, the arrays over-estimated lengths of ROH. WGS-based estimates of inbreeding were highly correlated (r > 0.98) with SNP array-based estimates, but only moderately correlated (r = 0.40) with inbreeding based on 5-generation pedigrees. On average, Group 1 horses had more heterozygous variants (P < 0.001), more total variants (P < 0.001), and lower individual inbreeding (FROH; P < 0.001) than horses in Group 2. However, the distribution of numbers of variants, allele frequency, and extent of ROH overlapped among all horses such that it was not possible to identify the group of origin of any single horse using these measures. Consequently, the Thoroughbred population would be better monitored by investigating changes in specific variants, rather than relying on broad measures of diversity. The WGS for these 185 horses is publicly available for comparison to other populations and as a foundation for modeling changes in population structure, breeding practices, or the appearance of deleterious variants.

Marsupial immune protection is shaped by enhancer sharing and gene cluster duplication of cathelicidin antimicrobial peptides

Marsupial neonates are born with immature immune systems, making them vulnerable to pathogens. While neonates receive maternal protection, they can also independently combat pathogens, though the mechanisms remain unknown. Using the sugar glider (Petaurus breviceps) as a model, we investigated immunological defense strategies of marsupial neonates. Cathelicidins, a family of antimicrobial peptides expanded in the genomes of marsupials, are highly expressed in developing neutrophils. Sugar glider cathelicidins reside in two genomic clusters and their coordinated expression is achieved by enhancer sharing within clusters and long-range physical interactions between clusters. These cathelicidins modulate immune responses and have potent antimicrobial effects, sufficient to provide protection in a mouse model of sepsis. Lastly, cathelicidins have a complex evolutionary history, where marsupials and monotremes are the only tetrapods that retained two cathelicidin clusters. Thus, cathelicidins are critical mediators of marsupial immunity, and their evolution reflects the life history-specific immunological needs of these animals.

Administration and detection of a multi-target rAAV gene doping vector in horses using multiple matrices and molecular techniques

Gene doping, which includes the non-therapeutic use of genes or genetic elements that have the capacity to enhance athletic performance, is prohibited in horseracing and equestrian sports. To provide a comprehensive assessment of matrix and detection techniques, a custom adeno-associated virus serotype 8 vector was designed to include PCR binding sites for multiple target genes and assay types. The vector was injected via an intramuscular route into two Thoroughbred horses and matrices collected at defined timepoints. DNA was analysed using 3 detection methods: qPCR, digital PCR, and NGS. Overall, there was a strong correlation across the different detection methods employed, although digital PCR was less sensitive at lower concentrations. High concentrations of vector were detected at early timepoints in plasma and whole blood, which rapidly dropped after 0.5 d to trace levels by 4 d and 9 d post-administration respectively, following a similar pattern to previous studies. Vector was detected in dried blood spots at lower levels than whole blood, but with a similar detection time. Detection in hair root bulbs in one horse was observed at over a month post-administration, which opens new avenues for future gene doping testing in humans and animals.

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.

Identification of novel cattle (Bos taurus) genes and biological insights of their function in pre-implantation embryo development

BACKGROUND

Appropriate regulation of genes expressed in oocytes and embryos is essential for acquisition of developmental competence in mammals. Here, we hypothesized that several genes expressed in oocytes and pre-implantation embryos remain unknown. Our goal was to reconstruct the transcriptome of oocytes (germinal vesicle and metaphase II) and pre-implantation cattle embryos (blastocysts) using short-read and long-read sequences to identify putative new genes.

RESULTS

We identified 274,342 transcript sequences and 3,033 of those loci do not match a gene present in official annotations and thus are potential new genes. Notably, 63.67% (1,931/3,033) of potential novel genes exhibited coding potential. Also noteworthy, 97.92% of the putative novel genes overlapped annotation with transposable elements. Comparative analysis of transcript abundance identified that 1,840 novel genes (recently added to the annotation) or potential new genes were differentially expressed between developmental stages (FDR < 0.01). We also determined that 522 novel or potential new genes (448 and 34, respectively) were upregulated at eight-cell embryos compared to oocytes (FDR < 0.01). In eight-cell embryos, 102 novel or putative new genes were co-expressed (|r|> 0.85, P < 1 × 10-8) with several genes annotated with gene ontology biological processes related to pluripotency maintenance and embryo development. CRISPR-Cas9 genome editing confirmed that the disruption of one of the novel genes highly expressed in eight-cell embryos reduced blastocyst development (ENSBTAG00000068261, P = 1.55 × 10-7).

CONCLUSIONS

Our results revealed several putative new genes that need careful annotation. Many of the putative new genes have dynamic regulation during pre-implantation development and are important components of gene regulatory networks involved in pluripotency and blastocyst formation.

Methanogenic patterns in the gut microbiome are associated with survival in a population of feral horses

Gut microbiomes are widely hypothesised to influence host fitness and have been experimentally shown to affect host health and phenotypes under laboratory conditions. However, the extent to which they do so in free-living animal populations and the proximate mechanisms involved remain open questions. In this study, using long-term, individual-based life history and shallow shotgun metagenomic sequencing data (2394 fecal samples from 794 individuals collected between 2013-2019), we quantify relationships between gut microbiome variation and survival in a feral population of horses under natural food limitation (Sable Island, Canada), and test metagenome-derived predictions using short-chain fatty acid data. We report detailed evidence that variation in the gut microbiome is associated with a host fitness proxy in nature and outline hypotheses of pathogenesis and methanogenesis as key causal mechanisms which may underlie such patterns in feral horses, and perhaps, wild herbivores more generally.

Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed

BACKGROUND

The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity.

RESULTS

To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders.

CONCLUSIONS

The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.

Widespread horse-based mobility arose around 2200 BCE in Eurasia

Horses revolutionized human history with fast mobility1. However, the timeline between their domestication and their widespread integration as a means of transport remains contentious2-4. Here we assemble a collection of 475 ancient horse genomes to assess the period when these animals were first reshaped by human agency in Eurasia. We find that reproductive control of the modern domestic lineage emerged around 2200 BCE, through close-kin mating and shortened generation times. Reproductive control emerged following a severe domestication bottleneck starting no earlier than approximately 2700 BCE, and coincided with a sudden expansion across Eurasia that ultimately resulted in the replacement of nearly every local horse lineage. This expansion marked the rise of widespread horse-based mobility in human history, which refutes the commonly held narrative of large horse herds accompanying the massive migration of steppe peoples across Europe around 3000 BCE and earlier3,5. Finally, we detect significantly shortened generation times at Botai around 3500 BCE, a settlement from central Asia associated with corrals and a subsistence economy centred on horses6,7. This supports local horse husbandry before the rise of modern domestic bloodlines.

Non-Synonymous Substitutions in Cadherin 13, Solute Carrier Family 6 Member 4, and Monoamine Oxidase A Genes are Associated with Personality Traits in Thoroughbred Horses

Retraining retired racehorses for various purposes can help correct behavioral issues. However, ensuring efficiency and preventing accidents present global challenges. Based on the hypothesis that a simple personality assessment could help address these challenges, the present study aimed to identify genetic markers associated with personality. Eight genes were selected from 18 personality-related candidate genes that are orthologs of human personality genes, and their association with personality was verified based on actual behavior. A total of 169 Thoroughbred horses were assessed for their tractability (questionnaire concerning tractability in 14 types of situations and 3 types of impressions) during the training process. Personality factors were extracted from the data using principal component analysis and analyzed for their association with single nucleotide variants as non-synonymous substitutions in the target genes. Three genes, CDH13, SLC6A4, and MAOA, demonstrated significant associations based on simple linear regression, marking the identification of these genes for the first time as contributors to temperament in Thoroughbred horses. All these genes, as well as the previously identified HTR1A, are involved in the serotonin neurotransmitter system, suggesting that the tractability of horses may be correlated with their social personality. Assessing the genotypes of these genes before retraining is expected to prevent problems in the development of a racehorse's second career and shorten the training period through individual customization of training methods, thereby improving racehorse welfare.

Ancient Mitochondrial Genomes Provide New Clues in the History of the Akhal-Teke Horse in China

This study analyzed ancient DNA from the remains of horses unearthed from the Shihuyao tombs. These were found to date from the Han and Tang Dynasties in Xinjiang (approximately 2200 to 1100 years ago). Two high-quality mitochondrial genomes were acquired and analyzed using next-generation sequencing. The genomes were split into two maternal haplogroups, B and D, according to a study that included ancient and contemporary samples from Eurasia. A close genetic affinity was observed between the horse of the Tang Dynasty and Akhal-Teke horses according to the primitive horse haplotype G1. Historical evidence suggests that the ancient Silk Road had a vital role in their dissemination. Additionally, the matrilineal history of the Akhal-Teke horse was accessed and suggested that the early domestication of the breed was for military purposes.

Population structure and genomic diversity of the Einsiedler horse

The breeding history of the Einsiedler horse is closely connected with the Benedictine cloister Einsiedeln. In the mid-nineteenth century, it was decided to use European Warmblood stallions for cross-breeding and to abandon the selection of stallions. Since that time, it has only been possible to trace back the origin of Einsiedler horses using maternal ancestry information. Here, we collected high-density genotype data for European Warmblood horses (Selle Français, Swiss Warmblood and Einsiedler) and Franches-Montagnes horses, the last native Swiss horse breed, to unravel the current population structure of the Einsiedler horse. Using commonly applied methods to ascertain fine-scale population structures, it was not possible to clearly differentiate the Einsiedler from other European Warmblood horses. However, by means of runs of homozygosity (ROH) we were able to detect breed-specific ROH islands for the Einsiedler horse, including genes involved in domestication and adaptation to high altitude. Therefore, future breeding activities should involve the screening of these breed-specific ROH segments, the revival of cryopreserved sperm and the selection of Einsiedler stallions.

slc26a12-A novel member of the slc26 family, is located in tandem with slc26a2 in coelacanths, amphibians, reptiles, and birds

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past.

A genomic hotspot of diversifying selection and structural change in the hoary bat (Lasiurus cinereus)

Background

Previous work found that numerous genes positively selected within the hoary bat (Lasiurus cinereus) lineage are physically clustered in regions of conserved synteny. Here I further validate and expand on those finding utilizing an updated L. cinereus genome assembly and additional bat species as well as other tetrapod outgroups.

Methods

A chromosome-level assembly was generated by chromatin-contact mapping and made available by DNAZoo (www.dnazoo.org). The genomic organization of orthologous genes was extracted from annotation data for multiple additional bat species as well as other tetrapod clades for which chromosome-level assemblies were available from the National Center for Biotechnology Information (NCBI). Tests of branch-specific positive selection were performed for L. cinereus using PAML as well as with the HyPhy package for comparison.

Results

Twelve genes exhibiting significant diversifying selection in the L. cinereus lineage were clustered within a 12-Mb genomic window; one of these (Trpc4) also exhibited diversifying selection in bats generally. Ten of the 12 genes are landmarks of two distinct blocks of ancient synteny that are not linked in other tetrapod clades. Bats are further distinguished by frequent structural rearrangements within these synteny blocks, which are rarely observed in other Tetrapoda. Patterns of gene order and orientation among bat taxa are incompatible with phylogeny as presently understood, implying parallel evolution or subsequent reversals. Inferences of positive selection were found to be robust to alternative phylogenetic topologies as well as a strong shift in background nucleotide composition in some taxa.

Discussion

This study confirms and further localizes a genomic hotspot of protein-coding divergence in the hoary bat, one that also exhibits an increased tempo of structural change in bats compared with other mammals. Most genes in the two synteny blocks have elevated expression in brain tissue in humans and model organisms, and genetic studies implicate the selected genes in cranial and neurological development, among other functions.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Integrated analysis of gut metabolome, microbiome, and exfoliome data in an equine model of intestinal injury

BACKGROUND

The equine gastrointestinal (GI) microbiome has been described in the context of various diseases. The observed changes, however, have not been linked to host function and therefore it remains unclear how specific changes in the microbiome alter cellular and molecular pathways within the GI tract. Further, non-invasive techniques to examine the host gene expression profile of the GI mucosa have been described in horses but not evaluated in response to interventions. Therefore, the objectives of our study were to (1) profile gene expression and metabolomic changes in an equine model of non-steroidal anti-inflammatory drug (NSAID)-induced intestinal inflammation and (2) apply computational data integration methods to examine host-microbiota interactions.

METHODS

Twenty horses were randomly assigned to 1 of 2 groups (n = 10): control (placebo paste) or NSAID (phenylbutazone 4.4 mg/kg orally once daily for 9 days). Fecal samples were collected on days 0 and 10 and analyzed with respect to microbiota (16S rDNA gene sequencing), metabolomic (untargeted metabolites), and host exfoliated cell transcriptomic (exfoliome) changes. Data were analyzed and integrated using a variety of computational techniques, and underlying regulatory mechanisms were inferred from features that were commonly identified by all computational approaches.

RESULTS

Phenylbutazone induced alterations in the microbiota, metabolome, and host transcriptome. Data integration identified correlation of specific bacterial genera with expression of several genes and metabolites that were linked to oxidative stress. Concomitant microbiota and metabolite changes resulted in the initiation of endoplasmic reticulum stress and unfolded protein response within the intestinal mucosa.

CONCLUSIONS

Results of integrative analysis identified an important role for oxidative stress, and subsequent cell signaling responses, in a large animal model of GI inflammation. The computational approaches for combining non-invasive platforms for unbiased assessment of host GI responses (e.g., exfoliomics) with metabolomic and microbiota changes have broad application for the field of gastroenterology. Video Abstract.

A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome

Background

Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression.

Results

An equine ASE analysis was performed, using integrated Iso-seq and short-read RNA sequencing data from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues from the Functional Annotation of Animal Genomes (FAANG) project. Allele expression was quantified by haplotypes from long-read data, with 42,900 allele expression events compared. Within these events, 635 (1.48%) demonstrated ASE, with liver tissue containing the highest proportion. Genetic variants within ASE events were in histone modified regions 64.2% of the time. Validation of allele-specific variants, using a set of 66 equine liver samples from multiple breeds, confirmed that 97% of variants demonstrated ASE.

Conclusions

This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues. Our results highlight the tissue-specific nature of allelic imbalance in the equine genome.

PCR-Based Equine Gene Doping Test for the Australian Horseracing Industry

The term 'gene doping' is used to describe the use of any unauthorized gene therapy techniques. We developed a test for five likely candidate genes for equine gene doping: EPO, FST, GH1, IGF1, and ILRN1. The test is based on real-time polymerase chain reaction (PCR) and includes separate screening and confirmation assays that detect different unique targets in each transgene. For doping material, we used nonviral (plasmid) and viral (recombinant adeno-associated virus) vectors carrying complementary DNA for the targeted genes; the vectors were accurately quantified by digital PCR. To reduce non-specific amplification from genomic DNA observed in some assays, a restriction digest step was introduced in the PCR protocol prior to cycling to cut the amplifiable targets within the endogenous genes. We made the screening stage of the test simpler and faster by multiplexing PCR assays for four transgenes (EPO, FST, IGF1, and ILRN1), while the GH1 assay is performed in simplex. Both stages of the test reliably detect at least 20 copies of each transgene in a background of genomic DNA equivalent to what is extracted from two milliliters of equine blood. The test protocol was documented and tested with equine blood samples provided by an official doping control authority. The developed tests will form the basis for screening official horseracing samples in Australia.

A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome

Background

Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression.

Results

In this work, we introduce and highlight the significance of an equine ASE analysis, containing integrated long- and short-read RNA sequencing data, along with insight from histone modification data, from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues.

Conclusions

This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues and foster a deeper understanding of the impact of allelic imbalance in equine health and disease at the molecular level.

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.

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.

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

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

Comparative Genomics Identifies the Evolutionarily Conserved Gene TPM3 as a Target of eca-miR-1 Involved in the Skeletal Muscle Development of Donkeys

Species within the genus Equus are valued for their draft ability. Skeletal muscle forms the foundation of the draft ability of Equus species; however, skeletal muscle development-related conserved genes and their target miRNAs are rarely reported for Equus. In this study, a comparative genomics analysis was performed among five species (horse, donkey, zebra, cattle, and goat), and the results showed that a total of 15,262 (47.43%) genes formed the core gene set of the five species. Only nine chromosomes (Chr01, Chr02, Chr03, Chr06, Chr10, Chr18, Chr22, Chr27, Chr29, and Chr30) exhibited a good collinearity relationship among Equus species. The micro-synteny analysis results showed that TPM3 was evolutionarily conserved in chromosome 1 in Equus. Furthermore, donkeys were used as the model species for Equus to investigate the genetic role of TPM3 in muscle development. Interestingly, the results of comparative transcriptomics showed that the TPM3 gene was differentially expressed in donkey skeletal muscle S1 (2 months old) and S2 (24 months old), as verified via RT-PCR. Dual-luciferase test analysis showed that the TPM3 gene was targeted by differentially expressed miRNA (eca-miR-1). Furthermore, a total of 17 TPM3 gene family members were identified in the whole genome of donkey, and a heatmap analysis showed that EaTPM3-5 was a key member of the TPM3 gene family, which is involved in skeletal muscle development. In conclusion, the TPM3 gene was conserved in Equus, and EaTPM3-5 was targeted by eca-miR-1, which is involved in skeletal muscle development in donkeys.

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.

Genetic structure and genome-wide association study of the traditional Kazakh horses

Horses are traditionally used in Kazakhstan as a source of food and as working and saddle animals as well. Here, for the first time, microarray-based medium-density single nucleotide polymorphism (SNP) genotyping of six traditionally defined types and breeds of indigenous Kazakh horses was conducted to reveal their genetic structure and find markers associated with animal size and weight. The results showed that the predefined separation between breeds and sampled populations was not supported by the molecular data. The lack of genetic variation between breeds and populations was revealed by the principal component analysis, ADMIXTURE, and distance-based analyses, as well as the general population parameters expected and observed heterozygosity (He and Ho) and between-group fixation index (Fst). The analysis revealed that the studied types and breeds should be considered as a single breed, namely the 'Kazakh horse'. The comparison with previously published data on global horse breed diversity revealed the relatively high level of individual diversity of Kazakh horses in comparison with the well-known foreign breeds. The Mongolian and Tuva breeds were identified as the closest horse landraces, demonstrating similar patterns of internal variability. The genome-wide association analysis was performed for animal size and weight as the traits directly related with the meat productivity of horses. The analysis identified a set of 60 SNPs linked with horse genes involved in the regulation of processes of development of connective tissues and the bone system, neural system, immune system regulation, and other processes. The present study is novel and introduces Kazakh horses as a promising genetic source for horse breeding and selection both on the domestic and international levels.

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.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Juvenile idiopathic epilepsy in Arabian horses is not a single-gene disorder

Valued for their temperament, beauty, athletic ability, and exhibition in the show ring, Arabian horses are an important component of the horse industry. Juvenile idiopathic epilepsy (JIE), a seizure disorder, is most often reported in Arabian foals from birth to 6 months of age. Affected foals exhibit tonic-clonic seizures lasting as long as 5 min and risking secondary complications like temporary blindness and disorientation. Some foals outgrow this condition, while others die or suffer lifelong complications if not treated. Previous work suggested a strong genetic component to JIE and proposed JIE to be a single-gene trait. In this work, we conducted a genome wide association study (GWAS) in 60 cases of JIE and 120 genetically matched controls, identifying loci suggesting JIE is not caused by a single locus. Coat color (chestnut, gray) phenotypes were used as positive control traits to assess the efficacy of GWAS in this population. Future work will attempt to future define candidate regions and explore a polygenic mode of inheritance.

Imputed genomes of historical horses provide insights into modern breeding

Historical genomes can provide important insights into recent genomic changes in horses, especially the development of modern breeds. In this study, we characterized 8.7 million genomic variants from a panel of 430 horses from 73 breeds, including newly sequenced genomes from 20 Clydesdales and 10 Shire horses. We used this modern genomic variation to impute the genomes of four historically important horses, consisting of publicly available genomes from 2 Przewalski's horses, 1 Thoroughbred, and a newly sequenced Clydesdale. Using these historical genomes, we identified modern horses with higher genetic similarity to those in the past and unveiled increased inbreeding in recent times. We genotyped variants associated with appearance and behavior to uncover previously unknown characteristics of these important historical horses. Overall, we provide insights into the history of Thoroughbred and Clydesdale breeds and highlight genomic changes in the endangered Przewalski's horse following a century of captive breeding.

Genetic diversity and signatures of selection in four indigenous horse breeds of Iran

Indigenous Iranian horse breeds were evolutionarily affected by natural and artificial selection in distinct phylogeographic clades, which shaped their genomes in several unique ways. The aims of this study were to evaluate the genetic diversity and genomewide selection signatures in four indigenous Iranian horse breeds. We evaluated 169 horses from Caspian (n = 21), Turkmen (n = 29), Kurdish (n = 67), and Persian Arabian (n = 52) populations, using genomewide genotyping data. The contemporary effective population sizes were 59, 98, 102, and 113 for Turkmen, Caspian, Persian Arabian, and Kurdish breeds, respectively. By analysis of the population genetic structure, we classified the north breeds (Caspian and Turkmen) and west/southwest breeds (Persian Arabian and Kurdish) into two phylogeographic clades reflecting their geographic origin. Using the de-correlated composite of multiple selection signal statistics based on pairwise comparisons, we detected a different number of significant SNPs under putative selection from 13 to 28 for the six pairwise comparisons (FDR < 0.05). The identified SNPs under putative selection coincided with genes previously associated with known QTLs for morphological, adaptation, and fitness traits. Our results showed HMGA2 and LLPH as strong candidate genes for height variation between Caspian horses with a small size and the other studied breeds with a medium size. Using the results of studies on human height retrieved from the GWAS catalog, we suggested 38 new putative candidate genes under selection. These results provide a genomewide map of selection signatures in the studied breeds, which represent valuable information for formulating genetic conservation and improved breeding strategies for the breeds.

Ancient segmentally duplicated LCORL retrocopies in equids

LINE-1 is an active transposable element encoding proteins capable of inserting host gene retrocopies, resulting in retro-copy number variants (retroCNVs) between individuals. Here, we performed retroCNV discovery using 86 equids and identified 437 retrocopy insertions. Only 5 retroCNVs were shared between horses and other equids, indicating that the majority of retroCNVs inserted after the species diverged. A large number (17-35 copies) of segmentally duplicated Ligand Dependent Nuclear Receptor Corepressor Like (LCORL) retrocopies were present in all equids but absent from other extant perissodactyls. The majority of LCORL transcripts in horses and donkeys originate from the retrocopies. The initial LCORL retrotransposition occurred 18 million years ago (17-19 95% CI), which is coincident with the increase in body size, reduction in digit number, and changes in dentition that characterized equid evolution. Evolutionary conservation of the LCORL retrocopy segmental amplification in the Equidae family, high expression levels and the ancient timeline for LCORL retrotransposition support a functional role for this structural variant.