Selection signatures for local and regional adaptation in Chinese Mongolian horse breeds reveal candidate genes for hoof health.

BACKGROUND: Thousands of years of natural and artificial selection since the domestication of the horse has shaped the distinctive genomes of Chinese Mongolian horse populations. Consequently, genomic signatures of selection can provide insights into the human-mediated selection history of specific traits and evolutionary adaptation to diverse environments. Here, we used genome-wide SNPs from five distinct Chinese Mongolian horse populations to identify genomic regions under selection for the population-specific traits, gait, black coat colour, and hoof quality. Other global breeds were used to identify regional-specific signatures of selection. RESULTS: We first identified the most significant selection peak for the Wushen horse in the region on ECA23 harbouring DMRT3, the major gene for gait. We detected selection signatures encompassing several genes in the Baicha Iron Hoof horse that represent good biological candidates for hoof health, including the CSPG4, PEAK1, EXPH5, WWP2 and HAS3 genes. In addition, an analysis of regional subgroups (Asian compared to European) identified a single locus on ECA3 containing the ZFPM1 gene that is a marker of selection for the major domestication event leading to the DOM2 horse clade. CONCLUSIONS: Genomic variation at these loci in the Baicha Iron Hoof may be leveraged in other horse populations to identify animals with superior hoof health or those at risk of hoof-related pathologies. The overlap between the selection signature in Asian horses with the DOM2 selection peak raises questions about the nature of horse domestication events, which may have involved a prehistoric clade other than DOM2 that has not yet been identified.

Integrative genomics analysis highlights functionally relevant genes for equine behaviour.

Behavioural plasticity enables horses entering an exercise training programme to adapt with reduced stress. We characterised SNPs associated with behaviour in yearling Thoroughbred horses using genomics analyses for two phenotypes: (1) handler-assessed coping with early training events [coping] (n = 96); and (2) variation in salivary cortisol concentration at the first backing event [cortisol] (n = 34). Using RNA-seq derived gene expression data for amygdala and hippocampus tissues from n = 2 Thoroughbred stallions, we refined the SNPs to those with functional relevance to behaviour by cross-referencing to the 500 most highly expressed genes in each tissue. The SNPs of high significance (q < 0.01) were in proximity to genes (coping - GABARAP, NDM, OAZ1, RPS15A, SPARCL1, VAMP2; cortisol - CEBPA, COA3, DUSP1, HNRNPH1, RACK1) with biological functions in social behaviour, autism spectrum disorder, suicide, stress-induced anxiety and depression, Alzheimer's disease, neurodevelopmental disorders, neuroinflammatory disease, fear-induced behaviours and alcohol and cocaine addiction. The strongest association (q = 0.0002) was with NDN, a gene previously associated with temperament in cattle. This approach highlights functionally relevant genes in the behavioural adaptation of Thoroughbred horses that will contribute to the development of genetic markers to improve racehorse welfare.

Hair follicle regional specificity in different parts of bay Mongolian horse by histology and transcriptional profiling.

BACKGROUND: Different morphological structures of hairs having properties like defense and camouflage help animals survive in the wild environment. Horse is one of the rare kinds of animals with complex hair phenotypes in one individual; however, knowledge of horse hair follicle is limited in literature and their molecular basis remains unclear. Therefore, the investigation of horse hair follicle morphogenesis and pigmentogenesis attracts considerable interest. RESULT: Histological studies revealed the morphology and pigment synthesis of hair follicles are different in between four different parts (mane, dorsal part, tail, and fetlock) of the bay Mongolian horse. Hair follicle size, density, and cycle are strongly associated with the activity of alkaline phosphatase (ALP). We observed a great difference in gene expression between the mane, tail, and fetlock, which had a greater different gene expression pattern compared with the dorsal part through transcriptomics. The development of the hair follicle in all four parts was related to angiogenesis, stem cells, Wnt, and IGF signaling pathways. Pigmentogenesis-related pathways were involved in their hair follicle pigment synthesis. CONCLUSIONS: Hair follicle morphology and the activity of ALP differ among four body parts in bay Mongolian horse. Hair follicles of the different body parts of the are not synchronized in their cycle stages. GO terms show a regional specificity pattern between different skin parts of the bay Mongolian horse. These results provide an insight into the understanding of the biological mechanism of the hair follicle in other mammals.

Refinement of Global Domestic Horse Biogeography Using Historic Landrace Chinese Mongolian Populations.

The Mongolian horse is one of the oldest extant horse populations and although domesticated, most animals are free-ranging and experience minimal human intervention. As an ancient population originating in one of the key domestication centers, the Mongolian horse may play a key role in understanding the origins and recent evolutionary history of horses. Here we describe an analysis of high-density genome-wide single-nucleotide polymorphism (SNP) data in 40 globally dispersed horse populations (n = 895). In particular, we have focused on new results from Chinese Mongolian horses (n = 100) that represent 5 distinct populations. These animals were genotyped for 670K SNPs and the data were analyzed in conjunction with 35K SNP data for 35 distinct breeds. Analyses of these integrated SNP data sets demonstrated that the Chinese Mongolian populations were genetically distinct from other modern horse populations. In addition, compared to other domestic horse breeds, the Chinese Mongolian horse populations exhibited relatively high genomic diversity. These results suggest that, in genetic terms, extant Chinese Mongolian horses may be the most similar modern populations to the animals originally domesticated in this region of Asia. Chinese Mongolian horse populations may therefore retain ancestral genetic variants from the earliest domesticates. Further genomic characterization of these populations in conjunction with archaeogenetic sequence data should be prioritized for understanding recent horse evolution and the domestication process that has led to the wealth of diversity observed in modern global horse breeds.

Refining the evolutionary tree of the horse Y chromosome.

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski's horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity.

Fast and slow myofiber-specific expression profiles are affected by noncoding RNAs in Mongolian horses.

The heterogeneity and plasticity of muscle fibers are essential for the athletic performance of horses, mainly at the adaption of exercises and the effect on muscle diseases. Skeletal muscle fibers can be generally distinguished by their characteristics of contraction as slow and fast type myofibers. The diversity of contractile properties and metabolism enable skeletal muscles to respond to the variable functional requirements. We investigated the muscle fiber composition and metabolic enzyme activities of splenius muscle and gluteus medius muscle from Mongolian horses. The deep RNA-seq analysis of detecting differentially expressed mRNAs, lncRNAs, circRNAs and their correlation analysis from two muscles were performed. Splenius muscle and gluteus medius muscle from Mongolian horses showed a high divergence of myofiber compositions and metabolic enzyme activities. Corresponding to their phenotypic characteristics, 57 differentially expressed long noncoding RNAs and 12 differentially expressed circle RNAs were found between two muscles. The analysis results indicate multiple binding sites were detected in lncRNAs and circRNAs with myofiber-specific expressed miRNAs. Among which we found significant correlations between the above noncoding RNAs, miRNAs, their target genes, myofiber-specific developmental transcript factors, and sarcomere genes. We suggest that the ceRNA mechanism of differentially expressed noncoding RNAs by acting as miRNA sponges could be fine tuners in regulating skeletal muscle fiber composition and transition in horses, which will operate new protective measures of muscle disease and locomotor adaption for racehorses.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds.

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV) induces IL-6 production through TAK-1/JNK/AP-1 and TAK-1/NF-κB signaling pathways.

Porcine reproductive and respiratory syndrome virus (PRRSV) mainly infects monocyte/macrophage lineage and regulates the production of cytokines to influence host immune responses. Interleukin-6 (IL-6) is originally identified as a B-cell stimulatory factor and has important functions in regulating immune response, hemopoiesis, and inflammation. In this study, we verified that highly pathogenic PRRSV (HP-PRRSV) infection up-regulated IL-6 production in vivo and in vitro. Subsequently, we demonstrated that HP-PRRSV infection activated JNK and NF-κB signaling pathways to enhance IL-6 expression. We further showed that TAK-1 was important in the activation of JNK and NF-κB pathways following HP-PRRSV infection. Moreover, AP-1 and NF-κB binding motifs were found in the cloned porcine IL-6 (pIL-6) promoter, and deletion of these motifs abrogated the activation of pIL-6 promoter by HP-PRRSV, suggesting that IL-6 expression is dependent on AP-1 and NF-κB activation. These findings imply that IL-6 induced by HP-PRRSV infection is dependent on the activation of TAK-1/JNK/AP-1 and TAK-1/NF-κB signaling pathways.

Selection in Australian Thoroughbred horses acts on a locus associated with early two-year old speed.

Thoroughbred horse racing is a global sport with major hubs in Europe, North America, Australasia and Japan. Regional preferences for certain traits have resulted in phenotypic variation that may result from adaptation to the local racing ecosystem. Here, we test the hypothesis that genes selected for regional phenotypic variation may be identified by analysis of selection signatures in pan-genomic SNP genotype data. Comparing Australian to non-Australian Thoroughbred horses (n = 99), the most highly differentiated loci in a composite selection signals (CSS) analysis were on ECA6 (34.75-34.85 Mb), ECA14 (33.2-33.52 Mb and 35.52-36.94 Mb) and ECA16 (24.28-26.52 Mb) in regions containing candidate genes for exercise adaptations including cardiac function (ARHGAP26, HBEGF, SRA1), synapse development and locomotion (APBB3, ATXN7, CLSTN3), stress response (NR3C1) and the skeletal muscle response to exercise (ARHGAP26, NDUFA2). In a genome-wide association study for field-measured speed in two-year-olds (n = 179) SNPs contained within the single association peak (33.2-35.6 Mb) overlapped with the ECA14 CSS signals and spanned a protocadherin gene cluster. Association tests using higher density SNP genotypes across the ECA14 locus identified a SNP within the PCDHGC5 gene associated with elite racing performance (n = 922). These results indicate that there may be differential selection for racing performance under racing and management conditions that are specific to certain geographic racing regions. In Australia breeders have principally selected horses for favourable genetic variants at loci containing genes that modulate behaviour, locomotion and skeletal muscle physiology that together appear to be contributing to early two-year-old speed.

Genomic inbreeding trends, influential sire lines and selection in the global Thoroughbred horse population.

The Thoroughbred horse is a highly valued domestic animal population under strong selection for athletic phenotypes. Here we present a high resolution genomics-based analysis of inbreeding in the population that may form the basis for evidence-based discussion amid concerns in the breeding industry over the increasing use of small numbers of popular sire lines, which may accelerate a loss of genetic diversity. In the most comprehensive globally representative sample of Thoroughbreds to-date (n = 10,118), including prominent stallions (n = 305) from the major bloodstock regions of the world, we show using pan-genomic SNP genotypes that there has been a highly significant decline in global genetic diversity during the last five decades (FIS R2 = 0.942, P = 2.19 × 10-13; FROH R2 = 0.88, P = 1.81 × 10-10) that has likely been influenced by the use of popular sire lines. Estimates of effective population size in the global and regional populations indicate that there is some level of regional variation that may be exploited to improve global genetic diversity. Inbreeding is often a consequence of selection, which in managed animal populations tends to be driven by preferences for cultural, aesthetic or economically advantageous phenotypes. Using a composite selection signals approach, we show that centuries of selection for favourable athletic traits among Thoroughbreds acts on genes with functions in behaviour, musculoskeletal conformation and metabolism. As well as classical selective sweeps at core loci, polygenic adaptation for functional modalities in cardiovascular signalling, organismal growth and development, cellular stress and injury, metabolic pathways and neurotransmitters and other nervous system signalling has shaped the Thoroughbred athletic phenotype. Our results demonstrate that genomics-based approaches to identify genetic outcrosses will add valuable objectivity to augment traditional methods of stallion selection and that genomics-based methods will be beneficial to actively monitor the population to address the marked inbreeding trend.

The distinct transcriptomes of fast-twitch and slow-twitch muscles in Mongolian horses.

Skeletal muscle is the largest organ system in the mammalian body and plays a key role in locomotion of horses. Fast and slow muscle fibers have different abilities and functions to adapt to exercises. To investigate the RNA and miRNA expression profiles in the muscles with different muscle fiber compositions on Mongolian horses. We examined the muscle fiber type population and produced deep RNA sequencing for different parts of skeletal muscles. And chose two of them with the highest difference in fast and slow muscle fiber population (splenius and gluteus medius) for comparing the gene expression profile of slow and fast muscle fiber types. We identified a total of 275 differentially expressed genes (DEGs), and 11 differentially expressed miRNAs (DEmiRs). In addition, target gene prediction and alternative splicing analysis were also performed. Significant correlations were found between the differentially expressed gene, miRNAs, and alternative splicing events. The result indicated that differentially expressed muscle-specific genes and target genes of miRNAs might co-regulating the performance of slow and fast muscle fiber types in Mongolian horses.

IL-10 knockdown with siRNA enhances the efficacy of Doxorubicin chemotherapy in EBV-positive tumors by inducing lytic cycle via PI3K/p38 MAPK/NF-kB pathway.

High levels of IL-10 expression in Epstein-Barr virus (EBV) associated tumors have been reported and it is likely to be important for maintaining EBV latency and EBV-associated tumors. The switch from the latent form of EBV to the lytic form in tumor cells can lead to tumor cell lysis. Here, we found that knockdown of IL-10 induced EBV lytic replication. Subsequently, we demonstrated that IL-10 knockdown activated BZLF1 promoter through PI3K-p38 MAPK-NF-κB signaling pathway. Interestingly, we verified that VEGF-A was required for IL-10 knockdown to activate PI3K signaling and the accompanying EBV lytic induction. Exogenous recombinant human VEGF-A induced PI3K activation and EBV lytic infection, and inhibition of VEGF-A signaling prevented the PI3K/AKT phosphorylation and EBV reactivation responded to IL-10 knockdown. Most importantly, IL-10 knockdown synergized with chemotherapeutic agent Doxorubicin to kill EBV associated tumor cells in vitro and repress EBV-positive tumor growth in vivo. Our results suggest that inhibition of IL-10 has the potential to serve as a new supplemental strategy for the treatment of EBV-associated tumors.

Donkey genome and insight into the imprinting of fast karyotype evolution.

The donkey, like the horse, is a promising model for exploring karyotypic instability. We report the de novo whole-genome assemblies of the donkey and the Asiatic wild ass. Our results reflect the distinct characteristics of donkeys, including more effective energy metabolism and better immunity than horses. The donkey shows a steady demographic trajectory. We detected abundant satellite sequences in some inactive centromere regions but not in neocentromere regions, while ribosomal RNAs frequently emerged in neocentromere regions but not in the obsolete centromere regions. Expanded miRNA families and five newly discovered miRNA target genes involved in meiosis may be associated with fast karyotype evolution. APC/C, controlling sister chromatid segregation, cytokinesis, and the establishment of the G1 cell cycle phase were identified by analysis of miRNA targets and rapidly evolving genes.

Analysis of horse genomes provides insight into the diversification and adaptive evolution of karyotype.

Karyotypic diversification is more prominent in Equus species than in other mammals. Here, using next generation sequencing technology, we generated and de novo assembled quality genomes sequences for a male wild horse (Przewalski's horse) and a male domestic horse (Mongolian horse), with about 93-fold and 91-fold coverage, respectively. Portion of Y chromosome from wild horse assemblies (3 M bp) and Mongolian horse (2 M bp) were also sequenced and de novo assembled. We confirmed a Robertsonian translocation event through the wild horse's chromosomes 23 and 24, which contained sequences that were highly homologous with those on the domestic horse's chromosome 5. The four main types of rearrangement, insertion of unknown origin, inserted duplication, inversion, and relocation, are not evenly distributed on all the chromosomes, and some chromosomes, such as the X chromosome, contain more rearrangements than others, and the number of inversions is far less than the number of insertions and relocations in the horse genome. Furthermore, we discovered the percentages of LINE_L1 and LTR_ERV1 are significantly increased in rearrangement regions. The analysis results of the two representative Equus species genomes improved our knowledge of Equus chromosome rearrangement and karyotype evolution.